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|. INTRODUCTION 


The variational effect of wind shears, convection, and temperature 
gradients on atmospheric density distributions characterize atmospheric 


turbulence. The temperature structure parameter, Cas, and the index of 
refraction structure parameter, Oe quantify two density distributions 


that arise in acoustic and electromagnetic propagation. [Ref. 1} 


The density distributions quantified by Giz and ee induce phase 


fluctuations in propagating electromagnetic waves. The density 
distributions of atmospheric turbulence degrade the spatial and temporal 
coherences of initially coherent electromagnetic radiation. Similar, more 
pronounced effects occur in the propagation of acoustic waves. [Ref. 1] 

A number of techniques for correcting these problems have been 
proposed, particularly with respect to electromagnetic propagation. These 
corrective techniques include the use of modular mirrors [Ref. 2], optical 
phase conjugation [Ref's. 2 and 3], optical and digital signal processing 
(Ref. 4], and simple avoidance of turbulent air masses. The appropriate 


timing of transmissions or suitable location of transmission and reception 


facilities accomplishes avoidance [Ref's. 5 and 6]. These techniques can 
be applied alone or in combinations. 

A knowledge of the history of the density or refractive structure 
along a proposed transmission path is required to effectively employ a 


number of these compensatory techniques. Time-averaged er or Ce 


profiles charactize this history [Ref’s. 1 and 7]. 

Wroblewski and Weingartner developed and Moxcey refined a high 
resolution, computer controlled acoustic sounder, or echosounder, system. 
Their computer controlled echosounder system uses a transmitted pulse of 
acoustic energy as a probe of atmospheric structure. The interaction of 
the acoustic wave packet with atmospheric turbulence is used to develop a 
time-averaged Ce profile for a short spatial range [Ref’s. 5, 6 and 8]. The 
acoustic sounder may also be used to verify certain long range optical 
measurements of turbulence [Ref. 6]. 

Development of the Oe profile uses a Bi: expression developed from 
the echosonde equation summarized by Neff [Ref. 7] and the empirical 
acoustic backscatter cross section per unit volume expression of Tatarski 


[Ref. 9] [Ref’s. 5 and 6]. The accuracy of the absolute values of this C,“ 


profile depends on the calibration of the acoustic antenna parameters 


contained in the Ca expression. The antenna parameters of interest are 
oF the efficiency of conversion of electrical to acoustic power, and oe the 
efficiency of conversion of acoustic power to electrical power. 


In this thesis acoustic energy backscattered off acoustically hard 


spheres was used to determine E, andE—,. Acoustic range equations were 
software implemented to perform the necessary calculations. Data and 
calculational results provided first order verification of the calibration 


process. 


Il THEORETICAL BACKGROUND 


A. C,* REVIEW 


Wroblewski and Weingartner applied the echosonde equation, 


re Efe glee ][S sen] 29] 0 


? 


summarized by Neff [Ref. 7] and the acoustic backscattering cross 


section per unit volume expression given by Tatarski [Ref. 9] as 








Ve C 2 
6(R,f}= 0.0039 k * —t. (2) 
To 
to develop the volume averaged expression 
2 : 1 P 
0.0039 E,E, 4% cr AG i, ; 


where 


e P_ is the electrical power of the reflected signal returned to the 


acoustical antenna array, 


e Eis the efficiency of conversion of the returned acoustic power to 


p 


r? 
e P, is the electrical power supplied to the acoustical array, 


e £, 1s the efficiency of conversion of P, to acoustic power, 


e « is the average attenuation per unit distance along the transmission 


path, 


e Ris the range from the array to the target or backscattering volume, 


e c is the average speed of sound along the transmission path, 
e tis the length in time of the transmitted acoustic pulse, 

A is the antennas aperture area, 

G is the antenna’s effective aperture factor, 

kK is the wavenumber of the incident acoustic energy, and 

qT, is the average temperature (in degrees Kelvin) along the 


transmission path. 


[Ref's. 5 and 6] 
The conditions at the time and location of each application determine 


the value of the termsT,, k, c, t, P., P,, Randa in the C,* expression. 


The terms A and G depend on the antenna design. 
Adverse field effects such as dust, debris and vibration loosened 


connections cause the remaining two terms, ES and Ee, to change over the 


life of the system. This chapter develops the theoretical expressions 


supporting a determination of = and 


B. ACOUSTIC RANGE EQUATIONS AND EFFICIENCIES 

In the following two sections the acoustic range equation of Neff 
(Ref. 7] is re-developed for targets of small area relative to the 
cross-section of the ensonifying field. Neff [Ref. 7], Skolnik [Ref. 10], and 


Probert-Jones [Ref. 11] were used for guidance. 


1. Qne Way Range Equation and EtGo 
As defined previously, the electrical power supplied to the 


acoustic array is P,. It 1s converted at an efficiency E, to an emitted 


acoustic power, 


Shoei ele (4) 
P, 1s spread over a 471 solid angle. 
A far field range R is defined by the condition 
A/2 >> ILI - IRI (5) 


for the geometry in Figure |. Ais the wavelength of the emitted acoustic 


W ave. 





Radius 8 


Acoustic 
Array 


Figure 1: Calibration (Far Field) Geometry 


The acoustic intensity, 
2 2 
ier. / STIR (6) 
describes the emitted acoustic power density of an isotropic emitter at a 


far fieldrange R. |, is in units of [W/m*] 


If the emitted power is directional, the intensity I, 1s modulated 
by a geometrical gain, 
G(O) = GIF(O)I 4 = G, IF(a,8)I 2 (7) 
Note this is not the same G found in equation (1). Now, 
F(0,0) 421, (8) 
and G, is thefore the centerline gain of the acoustic array. [Ref. 11] The 
centerline direction (0,0) is taken to be normal to the face of the acoustic 


array. The acoustic intensity |, is now written 


P 
ee a (9) 
L s 
4nR 





In areal propagating medium intensity attenuation occurs along 
the transmission path. The radial term e® accounts for this 


attenuation when the atmospheric path is assumed to be reasonably 


homogeneous. The acoustic intensity is then 


i. = Cea (10) 


ATR ~ 





A target of cross sectional area Avot intercepts a power 


Prat = ba Atat - (11) 
Expanding this expression with equation (10) yields 


& A 

ae -XR ~ tgt 
le 89) | [3 ] (12) 
For a target of cross sectional area Avat with a sufficiently 


small solid angle Avot /R* little change in intensity occurs about the 


direction 0. Equation (12) is now arranged as a one way range 


equation along the centerline direction (0,0); 





1 
PE : A 5 
t Ht "XR “tgt | (13) 
- e G 
R | AT ¢ - ; 


Isolating E,G, in the above expression yields 





7 ANR* or Fat (14) 
a Ce e 
fond P 
t tgt 


Electrical power 1s expressed as 


Pa 
p - Nrms” _ Mims. (15) 
? 7a 


where 


t eae is the empirical root mean square voltage, 
2 

a Ne isV_.*, and 

eZ is the electrical impedance. 


EG, is now written 





= 2 2 OR — tgt (16) 
E,G, = 47R e gt 
a (V rns}, A tat e 


The ratio BOE eatat is defined as the acoustic intensity |. 


received by the target. Therefore, 


= REIS (17) 
E,G, = 47R” =~ e 
ee (V ms), Ir : 


2. lwo Way Range Equation and Er 
For an acoustically hard target the total power reflected is Prat: 


The power reflected from the target surface has a normalized, 
directional intensity distribution described by the differential scattering 


cross section, a. 


The direction from the target to the array is the backscattered 
direction. The differential scattering cross section in the backscattered 


direction 1s called the normalized backscattered cross section, 6,. a, 1s 


defined as the power reflected toward the source per unit solid angle, 
normalized by the incident power density over a 471 solid angle. 


The backscattered intensity from the target is 


Re t 
| sat = at OF, (18) 
ADR 


An array of aperture area A at range R intercepts an attenuated, received 
acoustic power, 


¥ -XxR 
Boa len as (19) 


Substituting equation (18) into equation (19) yields 


P 
me AGE 


ae is converted to returned electrical power, es at areturn 


efficiency, €.; 


p.=P_E.. (21) 


Since electrical power is expressed as 


2 , 
pe = Vems = V ens (15) 


ra ae 


electrical power representing the returned acoustic signal from a 


target placed along an emitter’s direction (0,0) is written as 


CV rasdy - pee te -_— tat / Exh penta (2?) 
foe éy 4) a ATS 3) 4 . 


For a passive, linear acoustic array circuit the concept of 
impedance reciprocity is invoked [Ref. 12] to yield a two way range 


equation for scattering from a small target; 


1 
,; = a 

R= [Somet : Fob Se Aik (23) 
CV rasdy ad (An) 9 : 


This is written to isolate E ; 


CY ms), E,G, A Atat CO, 


A part of the return signal's power is unwanted noise. 


subtracting noise power from the returned acoustic power yields 


_ Ones (CIR Rs 


F = ee = eee 
/ (Yims), E,&,, Ay Ata OF, 


25) 


As return voltages are quite small they are analyzed after 
electrical amplification. Calling the electrical amplification a gain 


G, allows introduction of the final modification to c— 


_ (Minslea(Vimsin (ARS 2 8B 
(Vins Gee EeGa A A age Sb 





(26) 
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C. ATTENUATION AND RANGE DETERMINATION 

Employing equations (17) and (26) in the calibration process requires 
Calculating the attenuation « and the range R from the array to the 
target. 

1. Attenuation 

The calculated attenuation x for an atmospheric propagation 

path is assumed to be the combined result of molecular and classical 
absorption in the atmosphere. The classical absorption Xo is attributed 
to heat Conduction and viscous effects while the larger molecular 
absorptiona__, is attributed to the excitation of internal energy modes of 


atmospheric gases by the propagating sound energy. [Ref. 7] 
Information from Neff [Ref. 7], Businger [Ref. 13], Neiburger 
(Ref. 14], and Fuller [Ref. 15] is applied to make a reasonable determination 


ofa. The determination is made using the average measured temperature 


ye 


T_ in degrees Celsius along the propagation path, the measured 


atmospheric pressure P in millibars, and the measured relative humidity 


Rh in percent. 


The molecular attenuation coefficient, a, 


empirical relationship 


A 
0 2 Mie 
ne (0.18 fyetio) + | ey) | 


“mol ~ 304.6 ame 


Xo) 8 in units of [-dB/m]. 
Now, 


a 0.0078 if (Tx) 49 of 771-1/T*) 


is the maximum absorption at 


= 


f_ is in units of [Hz]. 


Also, 


e_ ER 
P  P-E(1 + Rh) 


no|— 


2 
(10 + 6600 (1002) + 44,400 (100 6)" ) P* 


is expressed by the 


Kez) 


(28) 


(29) 


(30) 


is the mole ratio of water vapor in the atmosphere. The water vapor 


pressure € iS expressed in millibars. 


Ee 1Q (9-4 - 2553/T, ) (31) 


expresses the saturation vapor pressure Ee ips 1s the average absolute 


temperature along the transmission path, as defined in Section A. 


T* = (1.87, + 492)/519, (32) 
and 
px =P/ 1014. (53) 
Finally, 
' ratio calles Tm ; Wee 


where f is the operating frequency of the echosounder in [Hz]. [Ref’s. 7, 
13, 14 and15] 


The classical attenuation coefficient, «.,, 1s approximated by 
x = 1.74% 10719 (fF)? (35) 
a, is in units of [-dB/m]. [Ref. 7] 


For a total attenuation 


Ky = Xo * Knol (36) 


and the relationship 


- xR = 10 log eek (37) 
the average attenuation per meter is finally written as 


X% + & 
ao ae mal 1. 14 (38) 


Equation (37) is derived from unit definitions and equation one of Neff 
[Ref. 7]. 
2. Range 
The coordinate system for the backscattered cross section 
determination uses the center of the target spheres as the coordinate 
system origin. Therefore, the radius a of a target sphere and the time of 


return a of an acoustic signal from the acoustic echosounder are used to 


determine the range R from the target to the echosounder; 


ct 
R = Wael ; (39) 





c 1s the average speed of sound along the transmission path, and 


t.=t 


r “receive | t 


transmit - (40) 


For the range of temperatures expected in normal echosounder operation 
the speed of sound in meter per seconds in dry air is 


eymOOS dT... (41) 


For moist air c becomes 


C\=) Coa Cony C1 +0.14¢6/P), (42) 


where €/P is defined in equation (30). [Ref. 1] 


D. BACKSCATTERED CROSS SECTION DETERMINATION 
Employing equations (17) and (26) in the calibration process also 


requires calculating«,, the backscattered cross section. Recall from 


Paragraph B.2 that the backscattered cross section is the differential 
scattering cross section in the direction of the energy source. 
1. Initial Conditions and Assumptions 

In determining the differential scattering cross sectione it is 
assumed the incident wave state is unchanged after scattering. That Is, 
for a stationary, acoustically hard target it is assumed 

K = Kincident ~ Xscattered » (45) 

and that incident k has good definition, i.e. the incident wave can be 
considered monochromatic [Ref. 16]. 

Selection of acoustically hard spheres as targets further 
simplifies the determination of the differential scattering cross section 


solution. This target selection sets the boundary condition at the target 


surface as the Neumann boundary condition 


ez 0, (44) 


Y is the velocity potential defined by the acoustic propagation velocity c 
as 
ceVl. (45) 


n lies along the direction normal to the target surface. Also, 


P= —- Py =—— (46) 
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Figure 2: Incident Wave and Target Description 


The target sphere is aligned along the centroid of the main 
acoustic lobe of the acoustic arrray echosounder, that is, along the 


direction (0,0), normal to the face of the array. The target sphere is 
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selected to be sufficiently small with respect to incident beam 
dimensions so that an essentially uniform intensity is intercepted by 
the sphere [Ref. 17]. Since the target intercepts a very small proportion 
of the centroid of the transmitted sounding lobe no accommodation for 
sounding lobe geometry over the target's angular area is made. 


Acoustic gain is simply taken as G, 


Target placement is also at a far fieldrange. This allows 
spherical incident waves to be treated as planar incident waves and 
places the target outside the range associated with the ring time of the 
echosounder [Ref. 8]. Placement of the target is also close enough to 
the array to support an assumption of transmission medium homogeneity 
along the transmission patn. 

Finally, the potential describing the target falls off faster than 
l/r. The describing potential for an acoustically hard sphere with 
radius a is 

co fia 
V(r) = (47) 
= Oe k 


Thus, forr>a V(r) falls off faster than 1/r. 


Note that with the exception of the specific surface boundary 
condition for acoustic waves the above conditions and assumptions closely 
correspond to those for plane electromagnetic waves incident ona 
perfectly conducting sphere, and for perfectly elastic nuclear scattering 
in acenter of mass coordinate system. [Ref's. 16, 18, 19 and 20] 

2. Ihe incident Wave 

The asymptotic form of an incident plane wave approaching a 
scattering target but beyond the influence of the target's potential V(r) is 
the undisturbed, or free, plane wave. In the spherical polar coordinates 
(r,,6) illustrated in Figure 3 a free, viel aac plane wave of 
unit amplitude propagating in the direction Q, given by (a.,8,) is 

fir,t) = ¢ ikrlcas( 2,)ca5( 6) + sin(a,)oos(@ - @,)] g-iwt (48) 


[Ref. 16] 
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Figure 3: Spherical Polar Coordinates 


For a wave of unit amplitude incident along the z axis , 0, = 0 


and the free wave f(r,t) is presented as 
{(r,t) =e ikr cos(B) g-iwt < o ikz p-iwt | (49) 
Note this is a solution with axial symmetry for the free wave equation 
[Ref. 18]. 
The wave equation describing the behavior of the undisturbed or 


free monochromatic wave {(r,t) is 


(v2 1 fie tyso. (50) 
ieee) is 

cis the speed of propagation of the free wave. If the wave 
function {(r,t) has a simple harmonic time dependence, i.e., 

{Pi eer (51) 
then for a wave number 

kK = w/c = 2TI/A (32) 
the free wave equation becomes the Helmholtz equation 

ive koe (53) 


(Rem. 
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The stationary solution {(r) of the free wave equation can be 
written as the product 
f(r) = $,(Q,) falas) {3(Qz) (54) 
for orthogonal curvilinear coordinates q,, qn, Gz [Ref. 16]. Recall an 


undisturbed plane wave has a cylindrical symmetry about the axis of 
propagation [Ref. 22]. Using the spherical polar coordinates 
Q,=1, Qo= 8, Qz=8, (55) 
the stationary solution describing an undisturbed plane wave propagating 
along the z axis 1S 
f(r) = $,(r) jo(a) . (56) 
Substituting {(r) into the free wave equation expressed in 
spherical coordinates and separating the coordinate dependencies yields 
the separated free wave equation, 


(kr)t+ tO ,eay - . eee (sin Bf, y, (7) 


j, dr dro} f,sing de 


Using a separation constant of I(1+1) the equation for the angular 
dependency of the free wave becomes 


ee eh sina“ 4.) + iG + 1) =O (58) 
Jnsing da 


2 | 


Letting 
cos(g) = x (59) 


in the angular equation yields 


z 
(= x2), - ox etd DIZ = 0, 60) 


The previous equation is known as the Legendre differential 
equation [Ref. 21]. The most useful solutions to the Legendre differential 


equation are the Legendre polynomials, designated by P(x). They can be 


generated by the recursion relationships 


POomzall, (61a) 
P s(x) = x (61b) 
and 
P(x) = [2-1/1] x Py_ (x) - [1-1 71} Py_oQ) (61¢) 
fRef. 21]. Therefore, 
fz (4) = P\(cos 8) . (62) 


From equation (57), the equation describing the radial 


dependency of the free wave 1s 





2 
eine moe IMBALS I A  G (63) 
gre or ar io ré IN 


LZ 


Changing scale to the variable s = kr yields the differential equation 





as oat! - BER T= 0 (64) 


[Ref. 17]. 
This differential equation has the spherical Bessel function j,(s) 
with amplitude A, as a solution; 
ty j(s) = A j,(9) (65) 


(Ref. 18]. 
The stationary solution of the free wave equation in spherical 
coordinates is now written as the partial wave 


f(r) = P(cos a) [A, j,(kr)) (66) 


The most general form of this free wave equation solution is a Fourier- 
Bessel series, 
oo e 

{(r,@) = Pa j,(kr) P(cos g) (67) 

[Ref's. 17 and 21]. 
A plane, free wave of unit amplitude propagating in the direction 
5 Y 1s then 
kK oo 
Se eels) | (GOS 4) (68) 


1=0 


we) 


The norm of the Legendre polynomials, the asymptotic expression 


a i) 
sin (kt 5 1) 
kr 


jkr) ~ (69) 


and the process detailed in Elton [Ref. 18] are used to determine Ay; 
A= il {2.et) . (70) 
Therefore, 
MZ Fil [21 +1] jy(kr) Py(cos 8) (71) 
[Ref's. 18 and 21]. 
5. Detection Conditions and Assumptions 
As the dimensions of the incident acoustic wave packet are 
greater than the size of the target it is expected that a portion of the 
incident wave packet continues to propagate past the target. The total or 
disturbed wave description is then a superposition of the incident and 
scattered wave descriptions. [Ref. 17] 
The acoustic array serves as both transmitter and receiver, or 
detector. This keeps the detector outside the path of continued 
propagation of the incident wave packet. [Ref. 17] Therefore, only the 


scattered wave description is of interest in the calibration of the array. 
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The placement of the target at a far field range simultaneously 
places the detector in the region of asymptotic behavior of the scattered 
wave. 

4. Ihe Disturbed Wave 

Once an incident wave described by the velocity potential V(r) 

enters the region of influence of a scattering target's spherical potential 


Vir) the wave equation becomes the disturbed wave equation 


{ y24[k2-vin]} vir) =0 (72) 
[Ref. 17]. Comparison of the form of the disturbed wave equation to the 
free wave equation gives an expectation of a solution of the form 


V(r) = V(r) Ve) . (73) 


Substituting /(r) into the disturbed wave equation and 
separating variables yields 


1 d d Capes 1 d-_éd 
~~ —— —— . fsin g —)jV = rt [kKo- Ving) + — =r =v, (74) 
V Sin B as 7 Ae a (r)) V, dr dr | 


The angular equation for the disturbed wave is identical in form 
to the angular equation for the free wave and independent of the Neumann 
boundary condition d War =0. Therefore, 


V5 \(8) = P,(cos 9). (75) 


Za 


Using the previous substitution of s = kr yields the now modified 


Bessel equation as the radial equation for the disturbed wave: 





+20, 7)- KERIO y= 0, (76) 


Following the solution development in Elton [Ref. 18] and 


applying the Neumann boundary condition at r = a leads to a solution 


Vy yr) = i f2t+ 1] Ej kr) - ay hie). 77) 


The relationship 


J , (ka) 


pipe 
a , (ka) 


a= ie" sing = (78) 


is defined from Morse [Ref. 17] and Bowman [Ref. 16]. Additionally, 


the j,(kr) are the spherical Bessel functions of the first kind, 


argument kr, and order ], 
hY(kr) are the Hankel functions of argument kr and order 1, 


the é, are the phase shifts of the usual partial wave description for 
spherical scattering solutions [Ref’s. 17 and 20], 

the j, (ka) are the first derivatives of the j,(kr), taken with respect 
to the argument kr and evaluated at r = a, 

n() (ka) = j.'(ka) + i y,'(ka), and (79) 

the y, (ka) are the first derivatives of the y,(kr), the spherical Bessel 


Functions of the second kind with argument kr and order |, taken with 
respect to the argument kr and evaluated atr=a. 


JS. 


The general series solution for the total disturbed wave (r) is then 

Mir)= & i [21 + 1] Py(cos @) [jy(kr) - a,’ my M¢kr)] . (80) 

©). Ihe Scattered Wave and its Amplitude 

The general asymptotic form for a disturbed wave /[(r,Q) 
resulting from a plane wave /’(r,Q) incident along 0, = 0 and scattered by 
a spherical potential V(r) is 

Vir.Q) = elk2 + 5(Q) elk fp (81) 

[Ref. 18]. This description of the asymptotic disturbed wave is a 
superposition of the incident plane wave e'? and a spherical, scattered 


wave s(Q)e""/r. The term s(Q) is called the scattering amplitude 


[Ref's. 18, 19 and 22] 
Now, recalling the series equivalence for the plane wave incident 


along 0, in equation (71), and comparing the general series solution for 


the disturbed wave in equation (80) to its asymptotic form in equation 


(81) gives 


s(Q)elk /p = - = Wives P,(cos 8) a, hi (kr) Bia) 


foil 


This is a series expression for the scattered wave. Therefore, 


\! 4 = oo a 
SC) — il is (2p P(cos a) a (83) 
Kk 1=0 
[Ref's. 16 and 18]. 





A nae 


The scattering amplitude s(Q) is related to the differential 


scattering cross section o(Q) =A, by 


6(Q) = Is(Qylé (84) 
[Ref's. 18, 19 and 22]. 


Applying the series expression for the scattering amplitude in 


equation (83) to the previous equation yields 


cA, = 6(Q) = mils 1S i [21+ 1] a,” P,(cos a)I¢. (85) 
g k 120 


For the backscattered direction @ = TI this gives a backscattered cross 
section of 





4yy ,@# 
= ee 12 
[Ref's. 16 and 17]. For the spherical target Avot = 11a‘, SO 
ay t 
6,2 12 (a+ ia, /?. (86D) 
(ca) as 
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lt should be noted a plane wave of unit amplitude incident from 
the positive z direction and backscattered along g = O would introduce a 


factor of (-1)! into the previous sum [Ref. 16]. 


E. SPHERICAL BESSEL FUNCTIONS AND THEIR DERIVATIVES 


Implementation of the previous equation necessitates calculation of 


the a,” The recursion relations 


j, (ka) = [-j,(ka) + cos(ka)]/ka , (87a) 


jy(ka) = {-L1+ 1Ej(kad)/ka} + jj(ka), = (870) 
y, (ka) = [-y,(ka) + sin(ka)]}/ka (87c) 
and 
y, (ka) = {-[1 + IIfy (ka)]/ka J + y \_ (ka) (87d) 
determine the j, (Ka) and y, (ka) needed for the calculation of a, in 


equation (78) [Ref's. 21 and 23]. 

Determination of the regular and irregular spherical Bessel functions 
at r= asupports the use of the recursion relations in equations (87). As 
Ka is real the irregular spherical Bessel functions are stable enough to use 


the recursion relations 


Je) 


y (ka) = -[cos(ka)]ka , (88a) 


y ,(ka) = [y,(ka) - sin(ka)]/ka (88b) 


and 


y(ka) = {[21 - tly ,_,(ka)]/ka} - y -36ka) (88c) 


in their determination [Ref's. 23 and 24]. 

Because the size of the argument of the regular spherical Bessel 
functions may be very small the method of negative continued fractions is 
appropriate to the fast and accurate determination of a converging 
solution to these functions [Ref.24]. The regular spherical Bessel function 


reciprocals are therefore determined by the negative continued fraction 





calculation 
be 1 Elfb, ,} fb, 2b, ,). tb) ,} fb) 9.) 4) (gg) 
jy(ka) — j,(ka) (yc ee [b; a] [b, 5b) a)... - 
Now, 
Dim= 20+m- 1)/ka (90) 
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and 


| on | 





nee...) Dyes De =O. - oT 
Lm a ee eee ee Oe 
| 
=f a ra 
Do _ 
bim-3z 
ial... 
or 
[Ref. 25]. 


While the regular spherical Bessel function reciprocal calculation 
of equation (89) would appear to involve an infinite number of terms, the 
calculation can in fact be terminated when the m" numerator term equals 
the m'” denominator term to the desired precision for the spherical Bessel 


function's solution. [Ref's. 24 and 25] 
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Hl. CALIBRATION PROCESS AND SOFTWARE 


Small quantities of data were collected in the anechoic chamber at 
the Naval Postgraduate School. This data aided first order verification of 
the calibration process and its supporting software. A brief description 
of the anechoic chamber is provided in Appendix A. 

The calibration process and data development occured in two stages. 
The first stage used equation (26) and two way propagation path 


measurements to determinea,E E,G,. The second stage determined E,G, 


using equation (17) and one way propagation path measurements. 


E_ was then determined froma,E_E,G,, E,G, and the calculateds, of 
r bert OO” ted b 


equation (86b). 


Because reciprocity was assumed for E, and E, 


Eee Eee (92) 


provided the desired calibration information. 
Supporting software was written to perform calculations required for 


calibration. The next section describes the development of this software. 
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A. CALIBRATION SOFTWARE DEVELOPMENT 

Calibration software was written in BASIC. It was used with a 
BASIC 4.0 compiler on a Hewlett Packard 510 computer. Computer 
performance was enhanced by an Infotek floating point processor. 

Notation in the programs sometimes differs from that in the thesis 
text as aresult of HP. keyboard limitations. Variable definitions are 
located at the beginning of subprograms or modules. These definitions 
are provided throughout a program as an aid in determining program to text 
variable equivalence. 

The main effort of software development was calculation of 
backscattered cross sections. This was done with equation (86b). 
Information available for summation loop verification was differential 
scattering cross section information. Differential scattering section 
software was developed first. 

1. Differentia attering Cro ection software Developmen 

The differential scattering cross sections « were calculated 
using equation (85) and equation (78). 
To support equation (78) employment the subprograms 


Calculating the values of the Legendre polynomials and the values of the 


Ul 
(AN 


regular and irregular spherical Bessel functions were developed first. 
Each subprogram was developed and tested separately in test programs. 
Each test program was designed to provide a tabular output for 
verification against values found in Abramowitz and Stegun [Ref. 23]. 

The Legendre values were calculated using equations (61). The 
irregular spherical Bessel functions were calculated using equations (88). 
The regular spherical Bessel functions were calculated using equations 
(89), (90) and (91). 

Following verification these subprograms were assembled with 
subprograms calculating the derivatives of the regular and irreqular 
spherical Bessel functions, the normalized differential scattering cross 
section «, the scattering modulus of Bowman [Ref. 16], and with three 
graphics subprograms. This assembly created the differential scattering 
cross section program found in Appendix B. 

The values of the regular and irregular spherical Bessel function 
derivatives at r=a were calculated using equations (87). 

Differential scattering cross sections were calculated using 
equation (85). The number of orders 1 required for the summation 


calculation in equation (85) was determined by generating tabular outputs 
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and determining the number of orders required for numerical stability to 
occur to 12 decimal places. This stability process was based on the 


concept of the phase shift 6, of each outgoing partial scattered wave 
becoming progressively smaller with increasing | until a, (equation (78)) 


became negligible for the precision desired [Ref's. 18, 19, 20 and 22]. 

Tabular outputs covered a ka value range of one to thirty at 
varying scattering angles, including @ = 11. This ka value range selection 
was based on the expected maximum target size to be encountered during 
actual echosounder employment. 

Twenty-one orders greater than the integer value of the 
argument ka were determined to be required for equation (85)'s 
summation. The number of orders required was greater than the usually 
cited relationship 

kax JI + 1) (93) 
gives. It should be noted that the precision criteria of « was greater than 
that normally encountered. [Ref's. 18, 19, 20 and 22] 
The graphics subprograms of the differential scattering 


cross section program provided outputs for comparisons to Figure 80 in 


Morse [Ref. 17] and Figure 10.10 of Bowman [Ref. 16]. These comparisons 
were used to verify the functioning and output of the differential 
scattering cross section program. 

The graphics subprograms used to generate Figure 4 provide 
linear polar plots of the differential scattering cross section «(()) 


normalized by the cross sectional area Avot of the target sphere. These 


linear polar plots show the normalized differential scattering cross 
section of acoustically hard spheres ensonified by plane acoustic waves 
incident from the left, i.e. from o = 7. 

The first view of each sheet in Figure 4 shows the entire 
normalized differential scattering cross section pattern, including the 
forward lobe, or “shadow zone’, alongg=0. The second view of each 
sheet is a magnified plot showing the scattering detail at a scale 
comparable to the normalized backscattered cross section value. Note 
the two views are identical on Sheet 1. The need for magnification with 


increasing ka becomes apparent on subsequent sheets. 
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Figure 4, sheet | of 3: Differential scattering Cross Section 
Program Output for Comparison to Figure 80 of Morse [Ref. 17]. 
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Figure 4, Sheet 2 of 3: Differential scattering Cross Section 
Program Output for Comparison to Figure 80 of Morse [Ref. 17]. 
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Figure 4, Sheet 3 of 3: Differential Scattering Cross Section 
Program Output for Cornparison to Figure 80 of Morse [Ref. 17] 
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Figure 4 was used for comparison to Figure 80 of Morse [Ref. 17] 
The radial scale of Morse was not graduated so verification consisted of 
a check for the correct number and approximate angular location of 
minimums and maximums in the differential scattering cross section 
patterns. 

The graphics subprograms used to generate Figure 5 provide 
semi-logarithmic plots of the scattering modulus S of Bowman [Ref. 16]. 

5S =k s(a)/ 447 (94) 

is the scattering modulus 5 of Bowman [Ref. 16]. s(a) is the scattering 
amplitude s(Q), as defined by equation (83), with a factor of (-1) 


introduced into the sum, that is, 


s(8) = vam. i = (-1)' (21+ 1) Pi(cos a) a’. (95) 
s(Q) is normally used to calculate «(0), per equation (84). The scattering 
modulus S of Bowman [Ref. 16] could also be used to calculate «(Q). Note 
the introduction of (-1)' into the summation corresponds to the case of an 
acoustically hard sphere ensonified by a plane acoustic wave incident from 


g=0. 
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Comparison of Figure 5 to Figure 10.10 of Bowman [Ref. 16] 
shows a match in both @ and S. This graphical match verifies the 
functioning and output of the differential scattering cross section 


program of Appendix B. 
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Figure 5, Sheet | of 3: Differential Scattering Cross Section 
Program Output for Comparison to Figure 10.10 of Bowman [Ref. 16]. 
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Figure 5, sheet 2 of 3: Differential Scattering Cross Section 
Program Output for Comparison to Figure 10.10 of Bowman [Ref. 16}. 
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Figure 5, Sheet 3 of 3: Differential Scattering Cross Section 
Program Output for Comparison to Figure 10.10 of Bowman [Ref. 16]. 
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2. Backscattered Cross Section Software Development 

The differential scattering cross section program of Appendix B 
was modified to yield the backscattered cross section program of 
Appendix C. With the modification the normalized differential scattering 
cross section is calculated for the backscattered direction only. Equation 
(86b) is used for this calculation. 

Additionally, the graphics subprograms of the differential 
scattering cross section program were replaced by one graphics 
subprogram. The new graphics subprogram was designed to provide a 


rectangular graph of the normalized backscattered cross section, ,. 
c,, is plotted as a function of ka for a specifiable range of values of Ka. 


Figure 6 was generated by the backscattered cross section 
program. It was used for comparison to Figure 10.11 of Bowman [Ref. 16]. 
This comparison served as a verification of the calculated, normalized 
backscattered cross section values. Note the Ka axis is labeled K*A. 

Figure 7 was generated to verify the normalized backscattered 
cross section profile was approaching a value of one as Ka left the Mie 
scattering regime and entered the “geometric optics’ region. Note the ka 


axis iS labeled K*A. 
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Figure 6: Backscattered Cross Section Program Output 
for Verification of Calculated Normalized Backscattered 
Cross Section Values by Comparison to Figure 10.11 of 
Bowman [Ref. 16]. 
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Figure 7: Backscattered Cross Section Program Output 
for Verification of the Normalized Backscattered Cross 
section Asymptotic Behavior with Increasing ka. 
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5. Acoustic Array Efficiency Software Development 

The acoustic array efficiency program in Appendix D determined 
FE, In general, each major section of Chapter II of this thesis and each 
set of analyzed data is represented by a subprogram of the acoustic array 
efficiency program. 

The acoustic array efficiency program of Appendix D includes the 
non-graphics subprograms present in the backscattered cross section 
program of Appendix C. It also includes subprograms designed to furnish 


Stored data, calculate attenuation and range using equations (27) through 


(42), and use developed information to calculate E.G, with equation (17). 
Equation (86b) was used to calculate «,. 


The appropriate subprogram and main program calculational 


results were then used in equation (26) to determine E.. E_ was used in 


equation (92) to determine E_ E,. 


B. DATA ACQUISITION AND SUPPORTING EQUIPMENT 
Data collection occured 28 November 1987, 30 November 1987, and 


7 December 1987. 
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Data collected 28 November and 30 November provided information 


needed for calculation of gE E,G5. Data collected 7 December allowed 
calculation of E,G, and the calculation of Enact 


Collected data is presented in Appendix E. 
1. [Two Way Propagation Path Measurements 
Two way propagation path data were collected 28 and 30 
November. Measurements were used to calculate the product o,E E,G,, 


Figure 8 provides a connection schematic for the data acquisiton 


equipment. 
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Figure 8: Data Acquisition Equipment Set-up 
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Four target spheres were available and measured for size. They 
included one hollow aluminum sphere (sphere number one) and three solid 
aluminum and brass spheres. Table (1) lists the data and the calculated 


values of Art = Tla*, ka, and ¢,. 6, was Calculated from equation (86b). 


Table 1: Target Sphere Dimensional Data and Normalized 
Backscattered Cross Section Calculational Results for 
Frequency of 5000 [Hz]. 


Target Sphere Sphere Ay [m7] ka c, 
No. Diameter {m] Radius {m] 
2046 +.0006 ae O509 1 Hh tay 9142 
Ms W013 4.0001 O35065 008060 4.619 1.000 
6) .07634+.00003 .03817 004577 3.479 1.051 
4 .06240+.00008 .035120 WOSVS6 2.843 6836 


Figure 9 shows the target spheres’ theoretical normalized 
backscattered cross section distribution plotted on a graph generated by 


the backscattered cross section program of Appendix C. 
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Figure 9: Target spheres Normalized Backscattered Cross 

Section Distribution as Calculated from Measured Diameters, 

Using an Assumption of Perfect Acoustic Hardness. 

The target spheres were suspended in turn in the anechoic 
chamber. Placement was made using the criteria in Section D.1t of 
Chapter I]. The range R from the array to a target was determined to be 
5.30 [m]. 


Temperature i and relative humidity Rh along the propagation 


path were determined with a Weathermeasure Temperature and Relative 


0 


Humidity Meter. T_ was usually about 21 [°c] with a relative humidity of 
approximately 50%. Atmospheric pressure P was read from a barogram 


supplied by the Naval Postgraduate School meteorological station. 
Specific i Rh, and P data are furnished in Appendix E. 
The unshrouded, close packed, hexagonal acoustic array of 
Moxcey [Ref. 8] was used to ensonify each target. 
A=N. 17.4 (96) 
determined the aperture area A of the array. N,= 19 was the number of 


speakers in the array and es .0381 [m] was the average speaker horn 


radius [Ref. 8]. 

Moxcey [Ref. 8] determined the half width of the array’s main 
lobe to be 12°. The angular half widths of targets two, three and four 
were all less than 5% of the main lobe's half width. Target one's angular 
half width was 11.5% of the main lobes half width. 

A sinusoidal input from a Hewlett Packard 3314 A function 
generator drove the array via apre-amplifier. Signal input from the 
HP3314A function generator to the pre-amplifier was set for 20 cycles at 


9000 [Hz]. The voltage input to the pre-amplifier’s power amplifier 


=) | 


terminal (see Figure 8) was 3.5 [V] for spheres two, three and four and 3.0 
[V] for sphere one. 

Twenty cycles is considerably less than the number of cycles 
expected in actual echosounder application. It was necessary to keep the 
pulse packet this small because the anechoic chamber’s corner was only 
about 1.5 [m] from the target mounting location. Voltage input for sphere 
one was reduced to 3.0 [V] to avoid clipping of its greater return signal 
voltage. 


Actual transmission voltage inputs V, to the array from the 


pre-amplifier were measured with a Nicolet 5091 oscilloscope at the 
pre-amplifier’s array terminal. The array was disconnected from the 


pre-amplifier for Ve measurements. V values are provided in Appendix E. 


They are also stored in the Trans30 and Trans35 subprograms of the 
acoustic array efficiency program in Appendix D. 
The Nicolet oscilloscope was connected to a Hewlett Packard 
7090 A measurement plotting system. This connection provided a hard 
copy plot of voltage vs. time traces. 
The returned signal was amplified and filtered. 


Pre-amplifier gain G, was determined as1 1,094 by Moxcey [Ref. 8]. 
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Filtering was done with a Wavetek Brickwall filter model 753 A to remove 
the low frequency normal modes of oscillation of the speaker drivers from 
the return signal. Filtering was accomplished at 0 [dB] gain and 5000 [Hz] 
for both low and high pass settings. 


The processed return signal voltages Vv. were displayed and 
measured on the Nicolet 3091 oscilloscope. V. measurements were made 


at the Out 2 terminal of the Wavetek Brickwall filter (see Figure 8). They 
were taken at 10 [1S] intervals for two full cycles along the asymptotic 
region of the return signal trace. 

The Nicolet 3091 oscilloscope was connected to the Hewlett 
Packard 7090 A measurement plotting system. Figure 10 shows the 
Hewlett Packard plotters traces for target sphere two's return signal 
voltages as printed and subsequently labeled. The figure is representative 
of the oscilloscope traces for each target return. Data for all targets are 
provided in Appendix E. 


The time of return a was determined as the time difference 


between the first positive maxima of the transmitted pulse packet and the 


first positive maxima of the target's return pulse packet. 
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Figure 10, Sheet | of 2: Target sphere 
Two's 5000 [Hz] Return Signal Trace 





Figure 10, Sheet 2 of 2 angel spermine. 
5000 [Hz] Return Signal Trace, Scaled Down 
For Oscilloscope Measurement. 
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r peak . 





Noise voltage V_ measurements were made with the same 
measuring equipment setup used for ve measurements. The target spheres 
were removed. V_ measurements were taken at a time of return t, 
comparable to the time of measurement of Me: 

‘, measurements were actually made on occasions well 
separated chronologically from the V. measurements. Therefore, they 


were proportionally corrected in line 2530 of the acoustic array 
efficiency program. The corrections account for the slight variations in 
input voltages and propagation path attenuations that occurred between 
the two measurement occasions. 

Figure |] shows arepresentative noise trace. Data from the 
noise measurements is presented in Appendix E and stored in the Noise30 
and Noise35 subprograms of the acoustic array efficiency program in 
Appendix D. 

Unless otherwise specified above the remaining data from the 
two way propagation path measurements are stored in the subprograms 
titled spherel, sphere2z, Sphere3 and Sphere4. All remaining data are also 


presented in Appendix E. 
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Figure 11: Representative 5000 [Hz] Noise Trace. 
3.0 [V] Supplied by the HP3314A Function Generator. 


2. Qne Way Propagation Path Measurements 
One way propagation path data were collected 7 December. 
Measurements were used to calculate E,G, with equation (1 7). 
Data from the one way propagation path measurements are 
provided in Appendix E. Data are also stored in the Gain30 and Gain35 


subprograms of the acoustic array efficiency program of Appendix D. 


AO 


For one way path measurements, 1.e. for the determination of 


E.G), an Ivie Electronics Inc. IE 30 A Audio Analyzer (Serial *805 A 426) 


was used to measure received target intensity |_ in [dB]. Intensity 

measurements were referenced to 10°'* [W/m?]. The acoustic signal 

supplied for this measurement was a continuous wave at 5000 [Hz]. 

Measurements were taken at the suspension site of the calibrating targets. 
A Fluke 8060 A True RMS Multimeter was used to measure the 


root-mean-square voltage (V ), supplied to the acoustic array at its 


rms 


input. (V___), was squared to determine (V_),. 


rms 
The Weathermeasure temperature and relative humidity meter 
was unavailable for E,G, data collection. Hygrothermographs from the 
Naval Postgraduate School's meteorological station were used to estimate 


the anechoic chamber’s relative humidity Rh and temperature Ton 


7 December. Estimates were based on comparisons to the 


hygrothermographs and measured Rh and T_ of 28 and 30 November. 
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Tee SSCILARS: 

Data analysis was performed with the aid of the acoustic array 
efficiency program of Appendix D. Data used in calculations is furnished 
In Appendix E. 


The product «,£ E.G, was generated from two way popagation path 
data and application of equation (26). The average of calculated a, S 
compared to the average of the experimental producta,F FG, gives a 
normalizing product EE,G,. The experimental product gE EG and the 


normalizing product E.E,G, are used to calculate experimental a's. 


These results are supplied in Table 2. 


Table 2: Target Sphere Experimental, Normalized 
Backscattered Cross Section Calculational Results. 


Target kalm?] Calc’d cE —G. bee Geax 
No. c, o;, 
11.62 9142 55.90 947 
2 4619 1,000 59.88 1.01 
5 3.479 1.051 61.62 1.04 
4 2.843 6836 37.98 644 
AVG. She 93.85 59.0 
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The experimental a, are plotted on a graph generated by the 


backscattered cross section program of Appendix C. This plot is Figure 


2. It is similar to Figure 9 and gives a, as a function of ka. Note ka is 


labeled K*A. 
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Figure 12: Target Spheres’ Experimental, Normalized 
Backscattered Cross Section Distribution Compared to 
a Calculated Backscattered Cross Section Distribution 
Curve Generated by the Backscattered Cross Section 
Program of Appendix C. 
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As ka increases the experimental a, distribution rises from just 
below the calculated, curve to just above it. 

E.G, was determined using one way propagation path data and equation 
che): E,Go, or ae was assumed to be the same for spheres two, three, and 


four because the supply voltage and target locations were the same. The 
angular half widths of target spheres two, three and four were between 


46% and 2.8% of the echosounder's main lobes half width. 


Table 3: Acoustic Array Product E,6, 


Determined from One Way Propagation 
Path Measurements and Equation (17). 


Target HP3314A E.G, 


No. Func n Genr 
Supply Volt. [V] 


a6. 94.70 
pe 55 Sia), I 
5 Gs: ees! | 
a j.2 Sa.) 


Calculated Sie E.G, and two way propagation path data were used in 


equation (26) to determine E,. «, values were calculated from the target 
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spheres’ measured diameters, 22. Avot was calculated as Tla*. Calculation 
results are furnished in Table 4. 
Table 4 E andE,£. &, Determined From Two Way Propagation 


Path Data and Equation (26). E,E, Determined From Assumption 
of Efficiency Reciprocity. 


Target Ne [m7] Calcd Ty Atat [m7] E Se 
No. c, 
05091 9142 04654 6456 A169 
2 008060 1.000 008060 6296 3964 
5 004577 1.051] 004810 6167 58035 
+ 003058 6836 002090 084] SAT 


The E.'s determined for Moxcey’s unshrouded hexagonal array are in 


the range of .58 to .65. This is above the value range of 0.5 + 0.1 
determined by Weingartner for his square acoustic array [Ref. 6]. 
Table 4 results were used to construct Figures 13 and 14. Figure 13 


gives EF asa function of the calculated, normalized backscattered cross 
section, s,. Figure 14 gives Ee as a function of the apparent target size, 
oA 
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Figure 13: Acoustical to Electrical Power Conversion Efficiencies, 
ee as a Function of the Normalized Backscattered Cross Sections, 








o,, Calculated from Measured Sphere Sizes. 
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Figure 14: Acoustical to Electrical Power Conversion Efficiencies, 
E_ as aFunction of the Apparent Target Sizes, Oot: 
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Recalling equation (20), 


5 -CCR Att 
: = oe 3 sium 20) 
aah i= |le 80) | |—$ | , ( 
and equation (12), 
ra AnrR* ’ 


serves as areminder that the received acoustic power, P_,, depends on the 
apparent target size, oH Atat 
Figure 14's trend and the dependency of P., ona, A,_, were used as 


b° tgt 


suggestions for the construction of Figure 15. Figure 15 is a plot of E, 


Target Number 





Figure 15: Acoustical to Electrical Power Conversion Efficiencies, 
E. as a Function of the Received Acoustic Powers, oe 
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versus P.,. P_. was calculated using the definition (from equation (21)) 
Pie= P/E, (97) 
and the relationship (frorn equation (15)) 
P= (Va (98) 
Note Figure 1's trend is for E, to increase with an increase in P_, for 


the range of target sizes investigated. 
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IV. CONC USIONS AND RECOMMENDATIONS 


A. CONCLUSIONS 


The larger magnitudes of E, determined for Moxcey’s hexagonal array 


Suggest it may be more efficient than Weingartner's square array. 
However, the drivers for each array are identical. 

The 0.08 to 0.15 differences in efficiencies are more likely the result 
of the difference in calibration methods. Weingartner’s measurements 
[Ref. 6] used steady state, continuous waves. The use of steady state, 
continuous waves allows the development of standing waves. This 
calibration investigation used pulsed waves. The use of pulsed waves did 
not provide sufficient time for standing waves to completely develop in 
the backscattered energy. {[Ref. 26] This calibration method should provide 
more accurate efficiency values since it more nearly approximates the 
echosounder system's actual field operating configuration. 

Additionally, steel braid and nylon lines were used to suspend the 
target spheres. These supports probably reflected enough acoustic energy 


to have a small influence on results. Other fixtures in the anechoic 
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chamber, such as support clamps and lights, may have also contributed 
small amounts of reflected energy to the returned signal. 

Calibration of the computer controlled echosounder required 
determination of the product EE, The calibration process assumed 
Eee 


: The calibration results in Tables 3 and 4 for spheres two, three, 


and four indicate E, varied with constant E.G), or constant re Such a 
variation negates the assumption of reciprocity for E, and E,. 


Figure 15 suggests the acoustic array’s efficiency E, increased in 


) 


response to an increase inP_and(V___). 
The acoustic gain for spheres two, three and four was probably not 


constant, as assumed. However, Go is the maximum for G(Q), and Eas 
inversely proportional to Go, per equation (26). An averaging process 


accounting for the slight differences in gain encountered by spheres two, 
three and four would have increased the data points spread observed in 
Figures 14 and 15. 

Additionally, the pre-amplifier's gain may not have been constant for 
the range of received voltage values encountered. This would have 


introduced either an apparent increase in E_ or dampened an increase In E 
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with an increase in(V___).. Such a gain variation could have effectively 


= 
negated the assumption of efficiency reciprocity. However, this 


pre-amplifier is approximately 0.1% linear for less than 10 [V] peak 


[Ref. 26]. Received voltages in this investigation were all less than 10 [V]. 


B. RECOMMENDATIONS 
Because efficiency reciprocity 1s now in question it is necessary to 


determine G, by calculation or measurement. Such a determination would 


allow the desired calibration quantity EE, to be written as 





EE, =E, oo | (99) 


One method of calculating G, would involve integrating the acoustic 


arrays calculated or measured intensity pattern (see Moxcey [Ref. 8]) for 


determination of an equivalent isotropic emitter. To determine Go the 


centroid intensity of the array’s emitted intensity pattern would be 
compared to that of the equivalent isotropic emitter. This would be 


repeated for each transmission power level of interest. 
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Calculation of G, would allow a determination of E, using—E,G,. Such 
a determination of 5 would allow a reciprocity comparison to Ge and 
support a more accurate evaluation of ERE 

The accuracy of E, and E,6, could be further refined as follows. 


Use a calibrating site that allows the use of a pulse train closer in 
length to actual field use. 

Make target return and noise voltages and site intensity 
measurements as chronologically close as possible. 

Measure the barometric pressure in the propagation path directly 
along with the relative humidity and temperature. 


The received and transmitted electrical powers for E_ and the 
transmitted electrical power for Eon were all measured at the 


pre-amplifier end of the pre-amplifier/array transmission line. This 
included transmission line power losses due to impedance as efficiency 
decreases. The present calibration process includes the transmission line 
impedance as part of the array impedance. 


The amount this impedance power loss actually affects E_ could be 


evaluated by determining the impedances seen in each direction at the 


68 


pre-amplifier's array terminal. This would be an opportune time to 
verify the validity of transmission line impedance reciprocity. 

Use the Nicolet oscilloscope to record a triggered readout from a 
calibrated microphone suspended at the target site at 45° to the centroid’s 


axis. This would allow E.G to be determined for a pulse packet rather 


than for a continuous wave. Coupled with a proper calibrating site 
selection this technique should minimize interference problems such as 
anechoic chamber corner reflections. 

A larger number of apparent target sizes spanning the desired range 


of apparent target sizes should be used to redetermine E_ for Moxcey's 
array. The E, distribution should be examined for a smooth and asymptotic 
behavior. Comparative E, values for Moxcey's array could be determined 


using the method outlined by Weingartner [Ref. 6]. 


Weingartner’s Square array’s E, could also be determined with this 
thesis’ calibration process and compared to Weingartner’s E, value [Ref. 6]. 

The refined E, curve in the now calibrated echosounder system should 
be used to develop an absolute a plot of an air mass. An established 


method such as a tower mounted vertical array of thermocouples would be 
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used simultaneously to generate a comparative absolute Gre plot of the 


Same air mass. 


C. SUMMARY 

Realizing the potential of the computer controlled echosounder for 
analyzing lower atmospheric turbulence required a calibration of the 
echosounder. This thesis described theory and software for performing 
calculations crucial to a backscattered cross section calibration of the 
echosounder. The backscattered cross section calibration was evaluated 
in this thesis. 


eae values indicate the evaluated calibration process requires some 


refinement. They also indicate the calibration method possesses 
sufficient merit to warrant further development and that the software 


performs as intended. 
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APPENDIX A 


ANECHOIC CHAMBER DESCRIPTION 


The anechoic chamber is located inroom 019 of building 232 at the 
Naval Postgraduate School. Acoustic research on sound sources, sound 
receivers and sound scatterers may be conducted in the chamber with a 
minimum of interference from wall reflection and external noise. 

Wall reflection is minimized by 102 centimeter deep wedges made of 
PF.612 fiberglass. The wedges are attached to the walls, ceiling and 
floor of the chamber. Approximately 142 cubic meters of fiberglass Is 
used to absorb 99% of incident sound eneray for frequencies greater than 
100 Hertz. External noise isolation is accomplished by separating the 
chamber's inner concrete block sides from the outer 12 inch concrete 
walls with a two inch thick lining of fiberglass and cork. 

A usable region of approximately 8.2 meters by 4.3 meters by 3.4 
meters is available. This region is floored by a grid of 225 wire cables, 


with each cable placed in a tension of 150 to 200 pounds per square inch. * 


*This Appendix was compiled from the description posted by the 
entrance to the anechoic chamber instrumentation and control room. 
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APPENDIX B 


DIFFERENTIAL SCATTERING CROSS SECTION PROGRAM 


Examples of the output of this program are Figures 4 and 5, located in 


ISSSSSSSFSSESSSSSSHFSSSSSSSSSSSSSHSSFSESSSSSHSSPLESSSSHSSESSSSSSHESSSSSEHSSESSSSF 


FEFFEE EEF HEEFE EF FFHE FEF EF EFEFTEEF TE FEEF FEET EE FEEF EE EF EF 4E +4 44444 et ets 
ACOUSTICALLY HARD SPHERE SCATTERING OF PLANE ACOUSTIC saw 


THIS PROGRAM CALCULATES AND GRAPHS THE DIFFERENTIAL SCATTERING 
CROSS SECTION AND/OR SCATTERING MODULUS FOR AN ACOUSTICALLY HARD 
SPHERE OF DIAMETER D ENSONIFIEO BY AN INCIDENT PEANE WAVERGE 


EReEQUENC var. 


FHEHHHEEFHEFEELHE EEL E FEF EAF ETE EE EE FEF EEFEAL EEE FE FE EE FEF EF tet tee et+etset 


[++++ VARIABLE DECLARATIONS AND DEFINITIONS +++++++++++t+++++4++ 4444444444 


C 


Freq 


Lmax 


Ka 


Chapter III. 
12 
20 | 
30 REM 
4Q REM 
SQ REM 
6Q REM 
70 REM 
BQ REM 
90 REM 
100 REM 
138 REM 
140} 
1SQ 
162 
170 
180 REAL 
192 
20@  =REAL 
210 
220 REAL 
230 
248 REQ 
250 
260 
270 
280 
290 
30@ = REAL 
310 
320 
330 
340 
550 REAL 
360 
370 REAL 
380 
390 
420 
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4 
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° 
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« 
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d2Z 


SPEED OF SOUND. 
FREQUENCY OF THE EMITTER. 
RADIUS OF THE SPHERE. 


ORDER OF THE FUNCTION CALCULATED 
IN A SUBROUTINE. It 1SeQeie 

AS A LOOP INDEX AND TO DESIGNATS 
THE ARRAY ELEMENT FOR THE PUNCG isis 
OF CORRESPONDING ORDER. 


MAX. SIGNIFICANT ORDER OF THESSSig 
DETERMINING THE DIFFERENT ae 
SCATTERING CROSS SECTION (DSCSe 
AND SCATTERING MODULUS =GSripe 


THE WAVE NUMGER, 2*PIl*PReO7e 
K*A, THE ARGUMENT OF THE REGUEAE 


AND IRREGULAR SPHERICAL GESsaE 
FUNCTIONS AND THEIR DERIVATIVES. 


41@ 
420 
430 
440 
45@ 
46@ 
47 
480 
490 
528@ 
51@ 
52 
53 
540 
550 
560 
570 
582 
598 
620 
612 
620 
632 
640 
55@ 
66@ 
672 
682 
690 
7Q0 
71 
72@ 
730 
748 
75 
760 
1a 
78 
730 
BOQ 
810 
820 
83@ 
842 
850 
86@ 
87@ 
88@ 
890 
90@ 


PGi \ot) 


REMe Otol? 


Benl YCS!|) 


REAL Oy(5St) 


REAL Sm(36@) 


| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
REAL P(51) | 
| 
| 
! 
| 
REAL 0s(360) 
| 

| 

{ 

| 

| 

| 

| 

| 

| 


ARRAY FOR THE REGULAR SPHERICAL 
Becsee FUNCTIONS <RSB) OF 
ARGUMENT KA AND ORDER @ THROUGH 
LMAX. 


faker On Ne OERIVATIVES OF THe 
ROE PGIHE DRS). 


fRRnte ORs HE IRREGULAR SPHERICAL 
BE coe mBUNCTTONS (158) OF 
ARGUMENT KA AND ORDER @ THROUGH 
LMAX. 


BRRAYeGOR He “CERIVAIIVES OF THE 
WSS lige SOUS ere 


ARRaY FOR THE LEGENDRE 
POSYNOMLAES OF ARGUMENT COSCPHI >. 


JHE ARRoeY REPRESENTING THE 
SCC IIERING MEBUEUS (SM. 


THE ARRAY REPRESENTING THE 
OIEPERENTI AL SCATTERING CROSS 
SECTION NGRMALI ZED BY THE CROSS 
SECILONGE Gneme OF THE SPHERE 
(DSCS >). 

NOTE GSCPHT)=4*SMCPH] )°27KA 2. 


varcoue (NOS Sen eadee 

THE ARRAY OIMENSIONS (30+21) NEED 
TO BE INCREASED TO ACCOMODATE 
KA>3@. 


PEEETHEt HEH tHe tHe tet ttt Eee ttt ett tt ett t et eet egg g teeta gett gttegeet ee ttttt 


Pettt MAIN PROGRAM Ft tet tt ett t ttt tet Ft tees tte tt tttetetetettetttet+tt+ett 


| 
fpeeeinit(C Frea ,A.Lmax Ka) | 
| 
| 


ieee Oitsetr\ max ,Ka 03+) ,J(+) Dye), 


| 
! 
Poeeedutpt(C Freq A,Ka,Sm(+) Ost#)) ! 


INPUT ANDO CALCULATE REQUIRED 
PERE NERS: 


VCs 2 PCe ) Smt) Os +). 
CAECUEATE THE SMC*) AND OS(s). 


OUTPUT STHE Snt*) ANO/OR DS *). 


PETE H HTH Ht He Het HHH tH EHH Et HH tt ee et tte ee tte tet t tate tet ee eteteetttettetetetttt 


12 


91Q (++++ MAIN PROGRAM CONCLUSION +4+4+444+44+4+4+4+4¢+4+444+44+4+4+4+4++4+4+44+4+4+14+4+4+4+44444 
92Q l 


330 PRINTER Shoe 

942 PRINT 

350 PRINT “END” ! ADVISES THE USER OF THEShigin 

362 PRINT ! PROGRAM'S CONCEUSTGHE 

370 END 

38 

3932 PHEHFEEHE EEE F AEF FH HELE E LAFF HATH HEHE EEE F HEF H THE H HH ttt ttt ttt ttt t ttt ttt etttt+ 
1000 ! 


1010 ItSSSSSSSESSSSSSSSSSHSESSSSSSSHSSSSSSSSSSSSSSSSSSSSSSSSSSSSSESSSSESSSSSSSSESS$F 
1020 ! 


74 


103@ 
1040 


105@ 
1060 
1070 
188 
109 
118@ 
111 
112@ 
1130 
114 
115@ 
1160 
117@ 
1182 
119@ 
120@ 
121@ 
1220 
123@ 
124@ 
1250 
1260 
1270 
128@ 
1290 
138 
310 
1320 
133@ 
134@ 
135@ 
1360 
137@ 
138@ 
1390 
140@ 
1412 
142@ 
143@ 
1440 
1450 
1460 
1478 
148@ 
1490 
150@ 
151@ 
152@ 
Hao 8 
1540 
oo 
1S6@ 
3/0 


| SSSSSFSSSS$SSSSSSSSSSSSSSSSSSSSSHSSSF$SSSFSSSSHSSFHSSSSSFSSSSSSSSSSSSSSSSSS 


SubwiInit<¢€ Freq Al max Ka?) 
FHFEFHE EEE TE ELEFHEE TE FE FEF EEE EEE EE EEF EF HEE EE HE EF EEE ttt ett etteggetete+est 


Tice UUme ReMUES TS THE APPROPRIATE INPUT YARIABLES TO DETERMINE 


REM 
REM 

REM 
men 
| 


C, FREQ, A, LMAX, AND KA. 


FEFEEFEEHFE FEF HFEF HEF EFHFE FE EEFEFEF HEHEHE EEE F HEE HEHE EF Ett ttt tte geggtetete tt 


!++++ MODULE VARIABLE DECLARATIONS AND DEFINITIONS ++++++++4+4+4++4++4+4++44++4 


2 


REAL Temp 


REAL Diam 


| 
ReMmL. Circ 
| 
! 


AS lENT hike tEMPERATURE etn DEG €. 


CIRCUME ERENCES OF STHE IGT. SPHERE: 


OnE ene Or mee wmRoe | SPHERE . 


PHEEEEHEHEF HEF EEF E EEE FFE EFE HEHEHE HEFT EEE HEHE EE TEE EF HEE EH Ht ttt tt gtetteeeettet 


Bie Re OURS ING NOS t RECUR ING SYSTEM PARAMETERS +++4++4+4+44+++4++44+++44+ 


* 


PRINT 


IF Choice=! THEN 
oP Senor CIN ee el 7S. 


INPUT C 


END IF 


PRINT 
PRINT 
INPUT 


IF Choice=2 THEN 


Peenoe se NT ER TEMP oOIN 
PUEGREES. CEESIUS." 
Temp 


Temp=Tempt+273.15 
C=2@.@5*SQR( Temp ) 


t 


e 
' 


i 
PRINTER [5 1 | 
| 
Choice=@ 
WHILE Choice<>1 AND Choice<>2 
PRINT 
Prin “DO YOU BESTRE-1Te INPUT” 
ERENT ~ Wl) tHE SPEED OF ."; 
PRINT “SOUND ," 
niliNie eeOrR <2.) HE AMBIENT AIR “*; 
PRINT TEMPERATURE?” 
PRINT 
INPUT Choice 
END WHILE 


END IF 

PRINT | 
meoNt) FLEASE ENTER THE ECHOSOUNDER" ; 
mei) )6 Uv FREQUENCY 6IN HZ.“ 

PRINT 

INPUT Freq 


7D 


INF UTS ARES REQUESTED ON THE] CRT. 


#** DETERMINING C «#** 


DIRECTLY ENTER C, 


CReCAeewila Le 16 PROMS TEMP . 


*** DETERMINING FREQ ««« 
OIRECTLY ENTER FREQ. 


1582 
1598 


16@2 
1612 
162@ 
163@ 
1642 
[650 
166@ 
167@ 
168@ 
ese 
170@ 
1712 
172@ 
Fie 
1748 
esd 
176@ 
177@ 
178@ 
1730 
1882 
1812 
1820 
1832 
184@ 
185@ 
1862 
187@ 
188@ 
189@ 
130@ 
191@ 
1320 
133@ 
1342 
bgae 
1962 
137@ 
1982 
est si 
2000 
2010 
2020 
2030 
2040 
2050 
2060 
2072 
2880 


| 
Chevee=0 | #**# DETERMINING KA «#*#* 


WHILE Choice<71 ANU Choice<oZ AND Chores 3 
PRINT 
PRINT “DG YOU DESIRE. To Theo 
eeu ae (THE SPER TC: 
PRIN) BESSsEee aRGUneNi ek sae: 
PRINT: C2) THE “SeHERE 2s. oc; 
PRINT CIROUNFERENGE. 
PRINT. OR (37° 0he sehese sex 
PRINT DIAnNETER (= 
PRINT 
INPUT Choice 

END WHILE 


IF Choice=) THEN | DIRECTLY ENTERS 
PRIN] “PEERSE VENTE heen 
INPUT Ka 

END IF 


IF Choice=2 THEN | CALCULATE KA FROM THE 
PRINT “PLEASE ENTER THE TARGET ~; ! TARGED C1 RGU eee) eee 
PRINE “CI REUMPERENCE. Hieeie 
INPUT Gace 
Ka=Circ#(10*°(-2))*Freq/C 

END IF 


IF Choice=3 THEN | OR GALCULATE KAV EROS Tie 
PRINT “PLEASE ENTER THE TARGE! ~~: !9TGRGE ROD aie ee 
PRINE SD PANETERSIN Ch. 
INPUT Diam 
Ka=Diame*( 1@°(-2) )*#PI*Freqg/C 
END Te 


| 

| *#* DETERMINING REMAINING #*#* 

) #68 SYSTEM PARAMETERS +++ 
K=2*PI*Freq/C | CALCULATE THE WAVE NUMBER. 
A=Ka/k | CALCULATE THE SPHERE Rabitise 
Lmax=INT(Kat+21) ! MAX. SUMMATION ORDER CALCULATED. 
| 


| 
PHAFF FEE FETE FE FEFE FEAT ELHE EEE EE HEE EE LE HFEFE FEE HEE EEE EE ee tee tee tteteeetettest 
| : 
I++4++ MODULE CONCLUSION t+ 4444 tt tt ttt tt te tte tet te te tettettetetttetttitit 
SUBEND 
LEH tHE HHH HFA TE HEF ATH HAF EAE HEATH EA FEAF EFF HATHA EAE FEE FFF EEE tet eggegteeteeeettetsettet 
| 
(SSESSSSESESSSESSSSSESHHSSHSHSSHSSHSSHSSSHSSSHSHSHSSSSSHSHSHSSSSSSHSSSSSSSHSHSHSESSSSSS$S 


76 


2090 
2100 


2110 
212@ 
2130 
214 
215@ 
2160 
2170 
2180 
219 
2200 
2210 
2220 
2230 
2240 
Z250 
2260 
2270 
2280 
2239 
2300 
2310 
Za2e 
2550 
2342 
Z500 
2360 
2370 
2380 
2530 
2400 
2412 
2420 
2430 
2449 
2450 
2460 
247 
2480 
2490 
2502 
Zo0 
2520 
2530 
2540 
2550 
2560 
25/0 
2580 
259@ 
2600 
261@ 


! 
ISSSS$SHSSSSFSHSHSHSHSSSSSHSSSSSSSSESSSSSSSSSSSHSSSSSSSSSSSSSSSSSSSSSESSSSSSSSSSSSS 


SUpmOu seur vena kha O70 a0 *) Dye) ¥(*#) PC *#) SmC*) Ds(*)) 

REM — FETT HHH HHT HHA THT ETH E FE FFF FFA FTF FEET HEEL HEHEHE tHe tt ttt ttt eettedt+ettte 
REM lle ODURE s@aueulbales THE ARRAYS SM(*) AND OS(*) REPRESENTING THE 
REM SCHecR NGS TOOCEUS tanh THE DIFFERENTIAL SCATTERING CROSS SECTION 
REM HOM eee ees On@ase oF PIONAL AREA OF THE SPHERE, RESPECTIVELY. 
REM —- FFF HFHHH HHH F HT tt HEHEHE FHH HHT HHH HS Ht Ht dette a ttgtetatgttettetettttt 
| 

er OUUC enn nee WeCkarknl LONG aNOwber INI TIONS +44+++++444++4+4+4+44444 


meme Fac FACTOR COMMON TO RSM AND ISM. 

REAL Rsm mesh ecesCOMPONENT Gr Tne 
SCATTERING MODULUS, SM. 

REAL Ism THE IMAGINARY COMPONENT OF THE 
Site 

peau Phi tie ROLGR ANGE ORES THe Axis 


| 
Reme Fact ane Ocenia omiemi >? & DSCrMI 
| 
POF TP RORGorT LON. Ge THE INCIDENT 
IPPEANE WAVE WITH ORIGIN AT THE 
ioe PER eMCeNTER ae COsvenl) 1S PHE 
DeOROUMENT OF THE LEGENDRE 
1 POLYNOMIALS P(*) OF ORDER @ 
Vee ie GO Ota ele icix 
| 
LEFF FELHEAFETAFEF HFEF FEET FE EEE FE EEL FET HEHEHE FFA FETE FE HE tee etttttetetet+etes 
| 
1++++ CALCULATING THE SM AND DSCS ARRAYS +4+++ttt+tttttt+ttttetttetetett++4te 


! 
Pale OrsbiLmax ,Ka ,.D0j(*),J(*)) SC ecCUus Mes HeeDERTUATIVE Ss OF THe 
PP RECUEpReSPHER Tene BESSEL, FUNCTIONS 
! OF ORDER @ THROUGH LMAX. 
| 
| 
| 


CALL Disb(Lmax ,Ka ,Dy(*) ,Y(*)) COECURATE IRE OFRIVATIVES OF THE 
TRRESUEOR SSPHERICAL Geoocl FUNC NS 
OF ORDER ® THROUGH LMAX. 
FOR Phi=@ TO 36@ 
PoeLyLeqgtL max Ka Phi .P(*)) CHeCUEGTE THE ARRAY PC®>. 
Rsm=Q 
Ism=2 
FOR L=@ TO Lmax 
mace CU Vem Zeta P Ce) )/CDG (Le Dg(L)tOy(L)eDy(L)) 
Rsm=Rsmt(FaceDy(l)) 
Ism=Ismt(FaceDj(L)) 
NEXT L 


CALCULATING THE REAL AND IMAGINARY 
COMPONENTS OF THe SM. 


— ——_ — —_ ——_ 


ae 


2620 
2630 


2640 
265@ 
2660 
267@ 
268@ 
2690 
2700 
2710 
2720 
2730 
2740 
275@ 
276@ 
2770 
2780 
2732 
2800 
281@ 
2820 
2850 


Fact=RsmeRsmtIsmetsm 
Sm( Phi )=SQR( Fact ) CALCULATING THE ARRAY REPRESENTING 
THE. SHC eal: 

CALCULATING THE ARRAY REPRESENTING 
THE DIFFERENTIAL SCATTERING €nuSe 
SECTION NORMALIZED BY THE GaGa. 
SECTIONAL AREA OF THE SPHERE: 


Ds( Phi )=4*( Fact )/(KatKa) 


— — cm — ee — 


NEAT etni 

| 
PHEEFEE EEF E ETE FEEL HEA LETHE HATE LA LETHE HALE E + HAE E+ te tte tte ettetetett+et+ett+t 
| 
l++++ MODULE CONCLUSION t+t++tttttttettttetttttttttettttettetttttttettttt 
SUBEND 

| 
PHFEEHEEHEE TEE LE FETE LEE HEAL HEEL HEHEHE ttt ttt ttt te ttgttttetetetettt+etet+tt+tt 
| 


ISSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSS 
| 


ve 


2840 
2850 
2862 
2870 
288 
Zo90 
2900 
2310 
acl 
23150 
2942 
ashe 
2960 
2970 
2980 
2990 
3822 
3@1@ 
3220 
3030 
3040 
3Q5 
3262 
3072 
3Q80 
303 
3122 
311 
312 
3130 
3140 
3150 


ISSSSSSSSSSSSSSSESSSHSSSESSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSS 


SUB Orsb(Lmax ,Ka ,03(*) ,J¢*)) 

REM FEEL EE EE AEE EEE EEEE EEE EEE E EEE EEE tHe ttt tte ttt tt ttetetttttte te ttt ttt++e tse 
REM THI GeMODULESCALCULATES THE DERIVATIVES OFeTHE REGULAR SPHERICAL 
REM Besse eeeUNGltUNas eos) OCF ARGUMENT KA AND ORDERS © THROUGH LMAX. 
REM THEst ae CUCe ONT OEE OERIUAN IVE INVOLVES THE RSB’S OF THE SAME 
REM AND PREV LOWS SORDER- 

REM —_- FH HHH FE FE FE FE FEF FETE EE HEE EF EHF HEH tH H ett tet tstseggtettttettttettt tt 
l 


Pee CHeCUEn hOn We eOeRiVamIVES OF THE ROB t++tttrttttt++etttt++++ttt t+ 


CALL Rsb(Lmax ,Ka,J(*)) CalCubat EiGeanse: URLUES Jt #1). 
ANGLE SUNT TS ECR DRSB CALC NS. 
CALCULATE THE INITIAL DRSB DJ(O@). 
CalCupate THES ReENAINING ORSE- 


RAD . 
Dj(Q@)=(-3(@)+COS(Ka))/Ka 
FER L=! TO Lmax 

By lo=s (Cte KadedthtJ3( be} ) 
NEXT L 


. 
. 
. 
* 


! 
! 
PELE THE AHEE HEEL AFLHE EA FEE HEE EEE EEE EEF EEE EET EFF E EE HEE EHH Ete tte tggettetetet+tetes 
{++++ MODULE CONCLUSION ++ tt 44+ ttt ttt tt ttt ttt ttt ttt tt tg ¢¢¢gttt¢teetet+ttet 
! 
SUBEND 

| 
PEELE HEFE HE EE HEEL EEE EE EEE EEE EEE HEE EEE FEE EAA HEE HE EE ete tegetetge¢eeet+etetets 
| 
ISSSSSSSSSSSSSSSSSSSSSSSSSHESSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSF 
| 


79 


3160 
S10 
318@ 
319@ 
3200 
3212 
3228 
3230 
3248 
3250 
3268 
5270 
3280 
3290 
3382 
3310 
3328 
3550 
334@ 
$350 
336 
337@ 
3380 
3330 
3400 
3410 
342Q 
3430 
344@ 
3452 
3460 
3470 
3480 
3492 
3580 
3510 
3520 
5530 
3548 
355@ 
3560 
3570 
3582 
3532 
3600 
361 
3622 
3632 
364@ 
365@ 
3662 
367@ 
3682 
369@ 
3702 


ISSSSSSSSSSSSSSSSSSSSSSSSSSSSSSHSSESSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSESSSSS 


SUB Rsb(Lmax ,Ka,J(*#)) 


BESSEL FUNCTIGNS 
THE 


REM —_ FE FFEFEE HEEFT FEF EE ESE FEF EEE EEE THEE EEE HEE E tt tet ttt et ttettettttett 
REM THIS MODULE CALCULATES THE REGUEAR SSP Reker 

REM (RSB) OF ARGUMENT KA ANO ORDERS @ THROUGH LMAX. 

REM CALCULATION USES THE CONTINUED FRACTION APPROACH. 

REM —— FFF HEEFT EFEE EEE TFET EE FETE FETE FE EE THEE HEE E EHH t tt ttt t te ttt te ttt ttttt 


l++++ MODULE VARIABLE DECLARATIONS ANDO DEFINITIONS +4+4+4+4++4+++++++++4++++4++ 


REALLY RI 


REAL 


Numer 


REAL Oenr 


REAL Scexp 


REAL Nflag 


REAL Oflag 


REAL Bnm 


REAL Binc 


l 
| 
| 
| 
| 
| 
| 
| 
| 
t 


THE EVOLVING RATIO OF J(0)/7J (te 


THE EVOLVING NUMERATOR FACTOR USED 
IN CALCULAT TIONG Shi 


THE EVOLVING DENOMINATOR FACTOR 
USED IN CALCULATION: GE Sai 


EVOLVING SCALING EXPONENT FOR VERY 
LARGE OR VERY SMALL VALUES OF Re 


@ FLAG FOR NUMR. 
@ FLAG FOR DENR. 


EVOLVING TERM USEO IN THE 
CALCULATION OF NUMR AND DENR. 


INCREMENT USED IN THE EVOLUTION 
OF SNM. ; 


PEEEEEEEEE EEE EEEE THEE THEE EEE EE FEE EEE EEE EEE EFE EEE EE FEE EHH ttt te tettttt+ttts 


l++++ CALCULATING THE RSB t+t+ttttttt tt ttt tttt ttt ttt ttt tt eet tteettetttetttt 


RAD 


J(Qd=(SIN(C Ka) )/Ka 


FOR L#] 
R1l=1.@ 
Numr=Q. 
Denr=@Q. 
Scexp=@. 

Nf lag=Q. 
Dflag=@. 
Bnm=!1.@/Ka 
Binc#2.@/Ka 


TO Lmax 


80 


* 


| 
| 
I 
| 
| 
| 
| 
| 
| 


ANGLE UNITS FOR CALC NG. 


ee® CALCULATE THE INITIAL 
eee RSB, J(Q). 


ett 


#aet 


eat 


### CALCULATING THE RATIO 
eet J(OI/I(L). 
INTTIGELZ I NGea ee 
INITIALIZING NUMR. 
INITIALIZING OENR. 
INITIALIZING SGExre 
INITIALIZING NFLAG. 
INITIALIZING OFLAG. 
INITIALIZING BNM. 
CALCULATE BINC. 


#e#7 


3712 
3720 
57350 
3742 
375 
3760 
3772 
3780 
3792 
3800 
3810 
3820 
3832 
3842 
3858 
386 
3870 
388 
383 
3902 
3910 
3920 
3930 
3948 
2150 
396@ 
3378 
3980 
3992 
4000 
4810 
4820 
4030 
4042 
4052 
4060 
4272 
4080 
4092 
4122 
4110 
4120 
4132 
4140 
4152 
4160 
417 
4180 
4190 
420 
4218 
4220 
4230 
4240 


FOR Ordent=!1 TOL 
Bnm=Bnmt+Binc 
IF Nflag=@ THEN 
Numr=Bnm-Nume 
IF Numr<>@ THEN 
Rl=R1i+*Numr 
Numr=1.8/Numr 


IF ABS(R1] )>10*%25@ THEN CALL 
poe 
Nflag=] 
END TF 
ESE 
Nflag=@ 
END 1 
EX] Ordent 


eee mecwlealE THE be =NAKED* * 
! * NUMR FACTORS OF J(@)/J(L). * 
| INCREMENT BNM. 

CHECK NUMR=0 FLAG NOT SET; 

| UPDATE NUMR IF NUMRCOO. 


mUroale Re TF UPDATED NUMR<>G. 

! PREPARE NUMR FOR NEXT EVOLUTION. 
PeCHeer SCALING OF UPDATED RL. 
Scale(Rl ,Scexp) 


PSlLESNRIEAG SO  THESEREVIGUSNSTEPS ARE 
INDEGERREO FUR@ONE INCREMENT OF GNM. 


! * CALCULATE THE REMAINING NUMR * 
Pee eoONDSDENR TERMS OF SICC)7I(L). * 
| NOTE WHEN NUMR=ODENR THEN 
! RE@J(O 7 IL). 


WHILE Numr<>Denr OR Nflag=1 OR Oflag=! 


Bnm=Bnmt+Binc 


IF Nflag=@ THEN 
Numrer=Bnm-Numr 
IF Numr<>@ THEN 

Ri=R1*Numr 
Numr=1.8/Numr 


Eee sR 7G 250 THEN CALL 
ELSE 
Nf lag=1 
END IF 
ESE 
Nf lag=@ 
END IF 


IF Of lag=@ THEN 
Denr=Bnm-Denr 
1—F Denr<>O THEN 

R1=R1/Oenr 
Denr=1 .@0/Denr 


Mere SiR b> 10° 250 THEN CALL 
BIeSe 
Dflag=1 
ENO TF 
ELSE 
Dflag=@ 
END IF 


END WHILE 


8 | 


INCREMENT BNM. 


| 
| 
| CHECK NUMR=@ FLAG NOT SET; 

[SUR OATE NUMA fr eNUMR< >). 

{ UPDATE RL IF UPDATED NUMR<>@. 

| PREPARE NUMR FOR NEXT EVOLUTION. 
CHECK SGeEING UF UPDATED RL. 
Scale(Rl ,Scexp) 


(ete SNmene<>O The PREVIOUS STEPS ARE 
! DEFERRED FOR ONE INCREMENT OF BNM. 


feCHECKSUENR=O>-EAG NOT SET; 
PUP PATE DENK STE SDENR< > . 


PVeOATE RE IF UPDATED DENRO. 
| PREPARE DENR FOR NEXT EVOLUTION. 
ISGHECKS SCALING GF WPDATED RL. 
Scale(R1] ,Scexp) 


elie eno OeT ne PREVIOUS STEPS ARE 
IHUERERREDSEOR ONG INCREMENT OF BNM. 


4250 
4260 


4270 
4280 
4290 
4302 
431@ 
4320 
4330 
4340 
4350 
4362 
4372 
4380 
439 
442 


J(LI=ICO/(R1L*#10*° (Scexp)) | «+e J(L) DETERMINED See 


NE Ct te 
| 
| 
PEEP HEEEEP EEE EEE ETE PELE EE EE EE TA PEE Et tt Pt EHH GH GR etic tt 
| 
l++++ MODULE CONCLUSION +444 44+44 4444+ 444444444444 444+ ¢¢ 44+ ¢+ t+ +++ 4+4¢4+4+4+ 444+ 
| 
SUBEND 
! 
LESH tttt ttt ttt ttt ttt tte t ttt ttt tt ttt tt ttt tt tt tt ttt tt te te 
{ 


(SSSSSSSHSSSSHSSSSHSSSSSSHSSSHSHSSSSHSHSSSSSHSHSHFSSSHSHSSSSSSHSSHSSHSSHSHSSHSHSHSSSHSSSSHSSSSS 
l 


B82 


4410 
4420 
4430 
4440 
4450 
4460 
4470 
4480 
4490 
4500 
4510 
4520 
4530 
4540 
4550 
4562 
4570 
4580 
4590 
4600 
4612 
4620 
4632 
4640 
4652 
4660 
4670 
4680 
4690 
4700 
4710 


1 
ISSSSSSSSSSSSSSSSSSSSSSSSSSESSSSSSSSESSESSSHSSSSSSSSSSSSSSSSSSSSSESSESSSSSSSS 
SUB Scale(R1l ,Scexp ) 
REM HEE EEEEEEE EEE EEE ELE EEE HEE EEE FEE EA HE FEE AEE Ett ete ttetttgttget+gttett+ttttts+ 
REM THIS MODULE SCALES RL TO MAINTAIN NUMERICAL ACCURACY IN THE 
REM ChHEGING CALCULATION OF RE. 
REM FHEEEHE EEE EEE EEE FEE HEHEHE FE HEE FEF EEE Et Ft eget tgttetgtt¢tttgsgtegtetet¢ettttettt 
{ 
Pee oOne ING tt tte e te tt te ete et tt tt tt et ttt ttt tt tte et + 
| 
MeRI>I1O° 250 THEN 
RI=R1]*1@*(-250 ) 
Scexp=Scexpt250 
eNOr iF 
teens ie (=250) THEN 
RI=R1#10*°2S50 
seexp=scexp-7250 
END IF 
| 
! 
PEE EEE EEE EE EEE EEE HEE EEE TEE E ETE EE EEE EE t tt gtetteetgt¢ettt¢eetttttttetttet+tetst 
| 
l++++ MODULE CONCLUSION +4++4+44+4+4+ +4 +44 44444444 4¢t+4¢t+4¢4¢4¢4¢4¢t+44+4¢ 444+ 444444444444 
{ 
SUBEND 
PEE EE ETE HEHEHE ETE EEE HEE EE EEE FE EEE EEE FE EE EEE EE FEE FEE ttt ttetettettttetgeetteet 
1 
ISSSSSSSSSSSSESSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSHSSSSS 
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4720 
4730 
4748 
475@ 
4760 
4770 
4780 
4792 
4800 
4810 
4820 
483 
4840 
4850 
4862 
4870 
4882 
489@ 
4900 
4912 
4920 
4930 
4942 
4952 
4968 
49370 
498@ 
4990 
5800 
5212 
5820 
5032 


| 
ISSSFSFSSSFSSHFESSSSFSSESSFSFSFESSFSSSSSSFSHSSSESSSSFSSEFSSESSSSSSSSFSSSFSSSFSSSSSSSSSSSFS 
SUB Disb(Lmax ,Ka ,Dy(*) ,YC*)) 

REM —_FHFEEEEFETHEFHEEFFTHE THEFT EEE FEE EHH T HHH + ttt te ttt tet tt tt+tt+tet tt tet ttt 
REM THIS MODULE CALCULATES THE DERIVATIVES OF THE ITRREGUEARS SP HE ries 
REM BESSEL FUNCTIONS OF ARGUMENT KA AND ORDERS @ THROUGH ERAXT ie 
REM CALCULATION OF THE DERIVATIVE USES THE [S87 S"0R pee aar eae 

REM  PREVIGUS ORDER: 

REM — FF HEHEHE HEFEH HHH FTF HTH H HHH H Ett t ttt tH sett ttt gtttt+etetetetttt+t ttt 


l++++ CALCULATING THE DERIVATIVES OF THE ISB +++++++++++++++4444++++444+4 


CALL Isb(Lmasa ,Ka ,Y(*)) CALCULATE THE TS8sUAcUE Se 

RAD 

Dy(@)=(-Y(@)+SIN(Ka))/Ka 

FOR L=1 TO Lmax 
Dy(EJ==C(( LF) Kaye ce]? 

NE Xi ae 


| 
| 
| 
! ANGLE UNITS FOR CABG Nee 
| THE INITIAL ISB DERIVATIVE 
! THE REMAINING ISB DERIVATIVE Se 

| 

| 
LHHEFET EHF E EFF HEF FETE FEE HE FETE THE FEE EEE HEHEHE tt ttt tt gt ttt ett gtteteetttt+tttt 
| 
[++++ MODULE CONCLUSION +4++4++++44+444tt4++tt+++4t++¢t¢¢t++4++4+¢t+¢+¢¢t+++4444 
] 


x 


SUBEND 

| 
PEEP tH ttt ttt tert tate ett ett e ttt tte t ett e ttt ttt tee tset ett t ttt tte tt ttt tt ttt 
I . 
ISSSSFESSESSSSESSSSSSSESSSESESSSHSSSSESSESSSHSSSFESESSHSSFSHSSHSSSSESSSSHSSESESSESSFESSS 
| 
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5040 
5@50 
5062 
5070 
5080 
5090 
5102 
S11 
el20 
5) 
514 
5152 
S160 
517 
5180 
5190 
5220 
5210 
5220 
3250 
5240 
S250 
5260 
5270 
5280 
5290 
5300 
oo!2 
5320 
S550 


| 
ISSSSSSSSSSSSSSSSSSHSSSSHSSSSSSSHSSHSSHSSSSSESSSSSSSSSSSHSSESSSSSFESFSSSSESHSSSSSHS 
SUB Isb(Lmax ,Ka,Y(*)) 

REM FEFEL EEE EE EEE EEE EEF EE EET EE FEE EAE EE ETE EE FEF Et te tgteetttteetttsetgeettt 
REM (os eOUUEE en_GuUbpteS THe TRREGULAR SPHERICAL BESSEL FUNCTIONS 
REM OF ARGUMENT KA AND ORDER L THROUGH LMAX. THE CALCUEATION 1S 


REM ACCONPEISHED Ey -ORWARD RECURSION. 
REM FHEFEFEFEEEEEEEE EEF E EEE TEE E ETE E ETH EE THE t eee gtgteteeteetetettt eet ttet 


I++++ CALCULATING THE ISB t+t+t4t4tt+tt ttt tt +t ttt ttt ttt ett ttttett+tttt+ettett 
| 


| 
RAD PSHNGEESUNT TS FOR CALC NS. 
Y(@)=-COS(Ka)/Ka Pi neetas SEGAL CULATED. 
eee = CO)SSINC Ka ))/ Kha 


pense=Z2 (0 imax I REMAINING 15675 CALC*D. 
ee Sct ty Ka) ey ela V Cle 2)) 
NEXT L 


| 
PE EEEEEEE TEE EEE EE EEE HELE FEET HEHE E EEE EEE AE HEHE EE EHH ete tg t tg ee tetetts¢ett+sett 
| 
I++++ MODULE CONCLUSION 44+ tt+4 ttt 44 ttt tt ttt tt tt Ft te ttt tg ttt gt ttt tetete+ete 
SUBEND 

| 

PEEEEEEE EEE EE FEAF AL HE EEE HEEFT FEAF EE HEE EEE EEE FEE EE Ete tts sgtgetgtttetttettses+ettet 
| 
ISSSSSSHSSSSSSSSSSSHSSSSSSSSSSSSHSSSSSSSHSHSHSSSSHSSSHSSSSSSSSSSSSSSSSSSSSSSSSSS 
! 
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5340 
550 
5360 
537 
538@ 
5530 
540@ 
541@ 
5420 
5430 
5448 
5450 
5460 
5478 
5480 
5490 
5500 
5512 
Bo2e 
55308 
5540 
5550 
5560 
5570 
558@ 
559 
56@@ 
S610 
5620 
563@ 
564@ 
5650 


ISS$SF$SSSSHSSSSSSS$SSSSS$SHSSHSSS$HSSSSSSHSSHSSHSSHSHSHSHSSHS$SHSSHSSHSSSSHSSSSSSSSSSSSSSESSS 

SUB Leg(Lmax ,Ka,Phi ,P(*)) 

REM FEFHEEE ETH EE EE EE EE EF EEE FEF EEE FEE FETE EE EE EEE HEHE EAE HEE HE het egggetttegetest 

REM THIS MOOULE CALCULATES THE LEGENDRE POLYNOMIALS OF ARGUMENT 

REM COS(PHI) ANO OROERS @ THROUGH LMAX. CALCULATION 1S ACCOMPISIsHas 

REM BY FORWARO RECURSION. 

REM FEFEFELHEHEETE HELE ETHEL EEE EE HEF E HE EE EEA EL EEE EEE HE Ett ee gtetettgg¢etetttyt 

f++++ CALCULATING THE LEGENORE POLYNOMIALS ++4+++++44+7 44 tpt 
| 


I 
DEG ! ANGLE UNITS FOR CALCRNce 
X=COS(Phi ) ! ARGUMENT OF THE LEG I PGinsse 

| 
P(Q)=1 ! INITIAL LEG. POL Seay aer nes 
PCl?=* 
FOR L=2 TO Lmax ! REMAINING LEG. POLS. CABG wae 

PUL RC 2e0 17D eX PC a ee ieee 

NEXT L 


! 
| 
PPE ttt ttt tet tt ttt te Fete tet eet ee tae tt tea EEE EE ee atte etcetera tr 
l++++ MOOULE CONCLUSION +444+44444444444¢44444444444744 074974 09 foe 
| 
SUBEND 
| 
LPP PEPEE ETH EET tet tt PEST ttt Fett ete tt et tte ete ett ee ee ee 


ISSS$FSSSSSSSSSSSSSSSSSSSSSSSSSSSSSS$SSF$SSSSSSSSSSSSS$SSSSSSSSSSSSSSSSSSHSF 
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566@ 
5670 


568 
569 
570@ 
5710 
B72 
5730 
5740 
5750 
5768 
5778 
5780 
a7 20 
580@ 
5812 
5820 
5830 
5840 
5850 
5860 
5870 
5880 
5898 
530@ 
531@ 
aaZ@ 
S50 
5942 
5100 
5360 
59978 
5980 
aaa 
B00 
5@1@ 
6020 
5030 
624@ 
605 
5262 
6270 
688@ 
5Q9@ 
618 
sic 
5120 
6138 
614@ 
6150 
616@ 
6172 
6182 
6190 


| SSS$SSSSSSSSSSSSSSSSSSSSSSSFSSSSSSSSSSSSSSSSHSSSFESFSSTS$SSSSSSSSSSSSSSSS 


SUB Outpt(C ,Freg ,A,Ka,Sm(*) Ds(*#)) 
REM FEFEEE THEE EEE LETHE EE EEE EEE HEE FE PETE HH ett tte tteetteeeteteettettetttt 


REM iS nOUUEEsCOnRULS THE ,OUIPUT OF THE CALCULATED DSCS AND SM IN 


REM GRArHI CAL PReseNTAT IONS. 


REM FHEFEFEF EE FETE FE FEF ETHEFEL EEE FETE FET HHH FH tt tt teste ttttettettttt+tetttt 


Pat eevee VaR nelCemUGClLARnalLONS AND DEFINITIONS +++4+4+++4++4+4++++++++4+++ 


REAL Desire 


REAL Rpt 


. 
. 
. 

? 
+ 


INTERACTIVE LOOP CONTROL. 


POLARPLT SCALING CONTROL 3 
REPETITION INDEX. 


A 


PEEEEET HEAT HTH FATALE FEF HEA FE FATE TAFE FETE EFT EFE TAFE FEF ttt tse ttsettetsttetts 


l++++ DETERMINING THE DESIRED OUTPUT +4+4+4+4++ttet¢tttttttttttetetttttttti+t 


Desire=@Q 


WHILE Desire<?1l AND Desire<>2 AND Desire<>3 


So UNigeGO YOU DESIRE” 


ERINT * im eon OG HE DSGs f° 
me tN Caer ot OF ine Sh. 


paliwie 9 OR) (5) BOTH? 
INPUT Desire 

END WHILE 

PRINTER IS 7@l 


WHILE Desire<4 
IF Desire=! OR Desire=3 THEN 
FOR Rat=) 710 2 
CALL Polarplt(Ds(*) Rpt) 
PRINT 
NEXT Rpt 
ENDO IF 


IF Desire=2 THEN 
CALL Reecolt<( Smt *)) 
PRINT 

END IF 


pont C =" 303 "M/S." 


Pies) Fee arreq; HZ. - 
Pine | ASE SAR M.S 
ie) ke A =" fKaa". “ 
Pen, NORMALIZED *; 
PRINT “BCKSCTR = “"3Ds( 18) 


PRINT CHRS( 12) 


| 
| 


PEO MING Vine ~OSCs:. 


PEOTIING THE SM. 


ECHO SIGNIFICANT CALCULATED VALUES 


AND INPUT PARAMETERS. 


elNGER FORM FEED. 


IF Desire=l OR Desire=2 THEN Desire=4 


IF Desire=3 THEN Desire=2 
FNOD WHILE 


8/7 


6200 
621 


6220 
6230 
6242 
6252 
6260 
6272 
6282 
6298 
63202 


| 
LEEEEEEFE FEET EEE EFFEEFEFEF HEE EEE FEE FEE HEHEHE te tte get tteeeteteteteteetettte 


l++++ MODULE CONCLUSION +4+44+444+4+444444 4444444444444 4¢44¢4+4+444¢4¢44+444444+444 
SUBEND 
PEEHEEEEE EHF TEFL HEA HEAT HE AEE FEF HE EEE HHA EEA FE FAFA FEE EFA EEF HEHE +e ttt etet+ataetetst 


ISSS$FSSSSSSSSSSSSSSSSSSSSSSSESSSSSSSSSSESSSSSSSSESSSSSESSSSSHSSSSSSSESSSSSSSSS 
| 
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6312 
6322 
6332 
6342 
6350 
6362 
6378 
638 
639 
6422 
6412 
6422 
6430 
6440 
6452 
6460 
6470 
6482 
6492 
6580 
6512 
6520 
6530 
6540 
6550 
6560 
657 
6582 
6592 
6600 
661@ 
6620 
663 
6640 
6650 
6662 
6672 
6682 
669 
6702 
671 
6720 
6730 
6742 


| 
[SSSS$SSSSSSSSSSSSSSSSHSSSSSSSSSSSSSHSSSSSHSSSESSSSSSSESSSSSSSESSSSSHSSESSSSHSSSF 
SUB Polarplt(Ds(*),Rpt) 
REM HEEFT HEHE FEET FEE EE FEE E FATE FEE EE FEF EEE SHH gt et gett ttetetet+tete+stts 
Ret Valse nOOUBesiae se Ine CALCULATED. NORMALIZED DSCS ARRAY DS(*) ON 
REM Seen Orlane een sea He FiRSh PEON INCEUMES THE ENTIRE DSCS, 
REM PAT eewiHe SeGONnb LON SHOWS HE SCATTERING DETAIL AROUND THE 
REM NORMALIZED RANGE. puOmrseeeCnloN TsrGOVERNED BY A REPETITION 
REM PNOEhe het Ne tS GONIROLLED IN THE CALLING MODULE (OUTPT?. 
REM FEFEEEEEEFEE FEF FHF FEF ELA FE FEF HE FEF EFF HEH Ft ett et gttttttte+tet+etsett 
| 


a OOULE. VARIABLE DECLARATIONS AND DEFINITIONS +++4+4+44+4+4++4+44++++44 
| 


REAL Xgumax THE MAXIMUM ABSCISSon VALUE IN 


GRAPHICS OISPLAY UNZTTS: 


THE MAXIMUM ORDINATE VALUE IN 
GRAPHICS BISPLAY UNITS: 


REAL Ygumax 


REAL Gumax THE MINIMUM OF XGUMAX AND YGUMAX; 
USEO eT OSDETERMINE PLOTTING AREA: 

REAL Dtmx Cian VALUE Oe tne ARRAY OS (+32 

REAL Endpt (HE Hake Rao tal SCALE VALUE IN 
USE BER INEDCUNITTS ; 

REAL Ring GRID RING SPACING IN USER UNITS. 

Bent lick RADIAL AXIS TICK MARK SPACINGS IN 


| 
| 
! 
| 
! 
| 
| 
| 
| 
| 
| 
PSUSeER DEFINED UNITS. 
| 
| 
PEAT HEF HEE THE EE HEE THEE EE HEHEHE H tH Ht ttt ttt ttt ge ttgtggtgteteeetsgeteegsteetes 
{ 
l++++ INITIALIZING THE PLOTTER +4+4+44++4++++¢+4+++t¢+4t4t¢+¢t¢¢¢4t¢+¢t¢+¢¢tt+esstts 
| 
GINIT 
BPeQTIER IS CRT ,“ INTERNAL” 
GRAPHICS ON 
GCLEAR 
| 
PHA E HEEFT HEHEHE HEE FFE H FF FH FHF HTH FFF Ett ttt +H ttt t+ ttt tettgstssggetssegeete+tte 
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675@ 
676@ 


6772 
6782 
6790 
6822 
6810 
6820 
6832 
6842 
685@ 
6862 
6872 
6882 
689@ 
6900 
6912 
692@ 
693@ 
63942 
6952 
6960 
69702 
698 
699@ 
7002 
721@ 
722@ 
703@ 
7042 
7852 
7060 
707 
7882 
789@ 
7180 
7110 
TWZze 
7130 
7142 
715@ 
716@ 
«172 
7130 
713 
7200 
12 Ve 
7220 
7232 
7240 
7252 
7260 
7270 
7280 
729 


ee GRAPH CONSTRUCTION AND LABELLING +4++++4+4++4++++++44++444+4444444 7477 
{ 
| 
Xgumax=1@@*MAX(1 ,RATIO) 
Yqumax=1@@#MAX(1,1/RATIO) 
Gumax=MIN( Xgumax , Ygumax ) 
| 
CST Zews f ** LABELLING MAJOR RAG TARSae 
LINE yee 
MOVE Gumax ,Gumax/2 
LORG 8 
PABELS HU. 
MOVE Gumax/2 ,Gumax 
LORG 6 
LABEL oP 2: 
MOVE @ ,Gumax/2 
LORG 2 
LAGE ly 
MOVE Gumax/2 ,@ 
LORG 4 
LABEL Ssh 2: 
** DEFINES THE TOTAL PLOTTING 
VIEWPORT S Gumax=6. 5, Gdmax=sS ** AREA IN GRAPH UNITS. ++ 
#* SCALING THE GRAPH TO ** 
** USER UNITS. #* 
DETERMINING MAX. RADIUS VALUE 
REQ"D FOR fHearlaie 


Dtmax=2 
IF Rpt=t THEN 
FORs 1=9 1G 368 
IF Dtmx<Os( 1) THEN Dtmx=Ds( I) 
NEXT I 
END IF 
Endpt=INT(Otmx+tl ) 
{F Emdpti<2 THEN. Eendpt=2 
SHOW -=Endpt -Endot. -Endpteacot ! ISOTROPICALLY SCALES THE GRAPH: 


— a camer a ee — 


Pine Teh 4 { ** MAJOR RADIALS CONSTRUCTION «+ 
Ring=INT(Endpt/Stt ) 

Tick=22* hing 

AXES VVerK AP Cis.) 0. S52 


EINE SYPE cs | ** RANGE RING CONSTRUCTION «#+# 
DEG 
LORG 7 
POR R=l710 Endat7Ring 
STNG. TPE 
MOVE @, R*Ring 
LABELOR Rima f LABELLING ThE RANGE 
EINE SIY¥RESS ; 
FOR Angle=@ TO 360 ! DRAWING THE RING. 
PLOT R*Ring*COS( Angle) ,ReRing*SIN( Angle?) 
NEXT Angle 
PENUP 
NEXT R 
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7300 
7310 
7320 
7330 
7340 
7350 
7360 
7370 
7380 
Very 
7400 
741 
7420 
7430 
744Q 
7450 
7460 
7470 
7480 
74390 
7500 
7510 
7520 
7530 
7540 
7550 
7560 
7570 
7580 
7340 
7600 
7610 
7620 
7630 
7640 
7652 
7660 


| 
FOR Angle=3@ TO 15@ STEP 30 { *#* 3Q@ DEGREE RADIALS #* 
IF Angle<>30 THEN | ** CONSTRUCTION. oe 
PLOT Endpt*COS(Angle?) ,Endpt*SIN( Angle?) ,-2 
DRAW Endpt*C0S(18@+Angle?) ,Endpt*SIN( 18@+Angle?) 
END IF 
NEXT Angle 


EENUP 


| 
CHEF EEEFE THEA HE FEL AFH EEE F EEE AFH FEF HE AEE FETA H AEE H+ He ttt ttttttttt+ttttt 
| 


eee eomhinG THe eneUUe neo Dols tttt+tttttte tt te tt ttt tt ttt tee ttt ttt te 
| 


INE TYPE 1 

DEG 

FOR I= TO 360 { PLOTTING RECTANGULAR COORDINATES 
PeOmeOsti )*C0S5( 1) 0st] o«SINCI> SINS EOEGR FORMAT. 

NEXT I 


ey ee emer eee. 4444 eee eee ta e444 edad d edad 44a ed 
| 
Pee PRINTING THE PLOT t4+4++++4+ +++ 4444444444 444444444444 44444444 4444444 
DUMP GRAPHICS #701 : 
| 
| caesar na eg erg rrr Orava WP 
| 
l++44 MODULE CONCLUSION +44 4444444444444 444444444 4+4444+4+4+4+4+44++44+4+4++4+4++44+4 
| 
Sar aiCs OFF | 
SUBEND 
eee eer rey 4st t44 tes 44444 44 deat eee deeded tet 
| 
[SSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSS$SSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSS 
| 


9| 


Tou 
7680 


7690 
770 


7710 
4720 
FU 
7742 
7750 
7760 
1770 
7780 
7730 
7800 
7810 
7820 
78350 
7842 
7850 
786 
7872 
7880 
7890 
790@ 
791@ 
7920 
1350 
7942 
7350 
7960 
7370 
7980 
VaA0 
8O82 
BAI 
8O20 
8232 
8242 
8050 
8268 
8Q72 
8Q8O 
8Q9@ 
8102 
8110 
812 
8130 
8142 
815 
8160 


I SSSSSSSSSSSFSSSSSSSSSSSFESSSHSSSSSFSSHSFSHSSFEFSESSSSSSSFSSSESSSFHFSSSSSSSSSSS 


SUB Recplt(Sm(*)) 
REM FHEFEETETE EEE TEE EEF EEE FEFE EEE FEE EEE HEE EE EF EEE EE EH 4 Ett t ttt tsttetaetttets+tt+ 


REM THIS MODULE MAPS THE CALCULATED SCATTERING MODULUS ARRAY SM(*) ON 


REM APSEMl-20G Freie 
REM FHEEFHEFEFEEEE ELE E FEF E EEE TELE EE EE EEF EEE HEE HEHE EE ttt ttt tt ttt atte tttetttt 


{++++ MODULE VARIABLE DECLARATIONS AND DEFINITIONS +++++444+7+44+47 777 ee 
| 


THE MAXIMUM ABSCISSA VALU Eee 
GRAPHICS OISPLAY Nahe 


REAL Xqumax 


THE MAXIMUM ORDINATE VALUE IN 
GRAPHICS DISPLAY UNiRSe 


REAL Yquma» 


REAL Dtmx THE MAX. VALUE OF THE ARRAY SM(*). 

REAL Otmn THE MIN. VALUE OF THE ARRAY SM(*). 

REAL Mxy THE MAX. ORDINATE SCALE VALUE IN 
USER DEFINED UNITS. 

REAL Mny THE MIN. ORDINATE SCALE VALUE IN 


| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| USER DEFINED Mis 
| 
| 
PEEHEEE HE EF HEHEHE TEFL HEHE EPEAT FELT EEE ELE EEE EEE FEE FE EEF EE EE + tte tgttgetetetttets 
| 
[++++ INITIALIZING THE PLOTTER +4+++4+++4+4t+4t44¢¢+4t+4+++t+44¢¢¢tt tt t+ t+tt+tt+t+ 
GINIT 
PEOITER, 1S CRI. INTERNALS 
GRAPHICS GN 
GCLEAR 
| 
PEE EEELEAEEEETEFA LEAF EFF EEE THEE EAE EEE E EEE THE TEE t Ett tga tg tettgetgtgesteteteteses 
[++++ GRAPH CONSTRUCTION AND LABELLING +4+4+4++4+4++4+++4++4+¢+4++4¢44+44¢4++4+4++4+4+4 
| 
Xgumax=10Q@*MAX(1 RATIO) 
Ygumax=10@*#MAX( 1,1 /RATIO) 


CSIZE 4 
| 
LORG 4 | «* LABELLING THE ABSCISSA Ax 1 Gaaee 
MOVE Xgumax/2 ,@ 
LAB EES Folin CUEGREES )- 


a2 


8170 
B18 
B19 
B200 
B21@ 
B22 
8232 
B24 
B25 
8260 
8272 
B28@ 
B239@ 
B38 
B31 
8320 
8338 
8340 
8350 
B36 
8376 
8382 
8390 
8402 
B41@ 
8420 
8450 
8440 
B45 
846@ 
8470 
8482 
8492 
8500 
8512 
8520 
8530 
8540 
8552 
BSE6@ 
857@ 
BS82 
B59 
BEG2 
BE1@ 
BE20 
8630 
B64 
B65 
BEES 
8670 
BE82 
8692 
8722 
B71 


DEG 

LOIR 90 

LORG 6 

MOVE @,Ygumax/2 

BEBeE “SCTRG MODULUS © 
LDIR @ 


VIEWPORT 15 ,Xqgumax-S ,190,Ygumax-S 


Dtmx=@ 
Dtmn=Sm( 180) 
FOR I=@ TO 180 
IF Dtmx<Sm(1I) THEN Dtmx=Sm(I) 
IF Otmn>Sm( I) THEN Dtmn=Sm(I ) 
NEXT I 
Mxy=INT(LGT(Dtmx ) +1 
Mny=INT(LGT(Dtmn) ) 


WINDOW @,18@0 ,Mny ,Mxy 


fees 50 1 eOeMoayeleliei.s 
MINES TYPE 3 
FOR I=Mny TO Mxy 
Bon J=i FO 8 
MOVES Ol tLGhie) 
PeOteisQent 61 (ioe —| 
PENUP 
NEXT J 
NEXT I 
FOR I=6@6 TO i18@ STEP 6@ 
MOVE I ,Mny 
aot boxy | 
PENUP 
NEXT I 
Sele Orr 
ine TYPE 1 
LORG 6 
FOR I=@ TO !8@ STEP 30 
MOVE I Mny 
EABEL I 
NEXT I 
LORG 8 
POR T=Niny TO Mxy 
MOVE OI 
BABEL 10°! 


HOvemOs 1 tLGT (5) 
Boece tl6 Cll) )72 
NEXT I 


oS) 


— — — ——— tae 


— — a — oo 


** LABELLING THE ORDINATE AXIS #«# 


*#* GRAPH CONSTRUCTION ** 
QERPINES THE STOMA PlOTT ING vane 
IN GRAPH UNITS. 


** ORDINATE USER UNITS ** 
** DETERMINED: *# 


Oi SOTRGETCGEE YS SCAEES THE GRAPH 
TomeeeR UNITS. 
PEOTS THESAX ES: 


Pe eeOrs THE HORT ZONT AE * 
"SECCORITHMIC GRID LINES [2A- Wo 
ever Ome yotCAlLe. * 


* 
* 


Pet omIHE VERTICAL GRID LINES 
‘eles lime WORTZONTAL SCALE: 


* 


** LABELLING THE ABSCISSA SCALE +# 


Pe nBeeeINe FHE OROINATE SCALE +**« 


POWER OF MIO ABELS. 


HGepeiMe NEAT FOWER OF 1@ LABELS. 


8720 
8730 
8742 
8750 
8760 
8770 
8780 
8790 
8800 
8810 
8820 
8830 
8840 
8850 
886 
8872 
8882 
8890 
83900 
8910 
8920 
8930 
83542 
8950 
8960 
8370 
8982 
8992 
9000 
9210 


CLIP ON 
| 
! 
PHL EHH HEHEHE FE FE FEF H HET A HEHEHE HEHE EHH HE Ht HHH tte ttt te tet eetgtetet+etetttt ttt 


[++++ PLOTTING THE CALCULATED SCATTERING MODULUS +4+4+4+4+4+4+4++4++4++4++++++4++ 


LINE irre | THIS PLOT 1S FOR THE PLANESMhes 

FOR 1-070 seg | INCIDENT FROM THE RIGHT (0 OEGHes 

PLOT. I EGTCom( tailed. i USE SMC I) IN PLACE OF SHiG@iaaie 

NEA ot | FOR A PLANE WAVE INCIOENT FROM THE 
l LEFT C1S0 DEG ses 


i se Ss Sn nin ie i in i ie in ie i i in i in in a ain i in ie a in din in ie ne dn ain ain ie is i in ie ie i ie ie i ie dc in ie ie Sn ie ie i ie i nn 
! 
b++++ PRINTING THE PLOT +4+4+4+44+ +++ 444444 t 4444 $+ 44+ 4¢¢++44+4+4¢¢4¢+4 tt tte et tt ttt 
OUMP GRAPHICS 470% 
| 
PHHHHHE HHH EHH HH HHH HEH HHH ttt tHe HH tt ttt tt ttt ttt te tetttt+et+settet++et+ttett 
! 
l++++ MOOULE CONCLUSION ++ 444444444444 4444444 44+4+4+44+¢4+444+4+44+4+¢444+4+4+4+444+ 
| 
GRAPHICS OFF 
SUBEND 
PEFFHE EEE EEF HEF EE FAFFF EAL HEATH E FEAF EFA FE HH $$$ ttt et tttsettsets¢tt+ett++t++t+t++tt 


| SSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSFSSSSSSSFESSSSSSSSSSES 
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APPENDIX C 


BACKSCATTERED CROS55 SECTION PROGRAM 


Examples of the output of this program are Figures 6 and 7, located in 


ISS$SSSESSSSSSSSSSSESSSSSSESSSSSHSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSESSESS 


FHEFEFHFFEEEFE FEET EEF EEE EAE EEE FEF FE EFEE EAE HEHEHE HE HEHEHE Het tte tet tt ttt tt 
PCOUSTICAEEN Hanes srHeRe BocCKoeqmieRet NG OF PIEANE ACOUSTIC WAVES 


THIS PROGRAM CALCULATES AND GRAPHS THE NORMALIZED, BACKSCATTERED 
CxO Som see Ne onemceuscl TCALLY SHARDS SPHERES OF A SELECTED RANGE 
CF QIANETERS O ENSONIFIED BY AN INCIDENT PLANE WAVE. 


$HHEEFHEFEAF HEHEHE FFE HELE FEAF E EEE E FEE HE FEFEE HEE HE tHe ee tte ttt tetetetttetett 


f++++ MAIN PROGRAM VARIABLE DECLARATION AND DEFINITION ++4+++4++4+++++++44++ 


Nbs(3@1 ) 


Ka 


Kamin 
Kamax 


Sep 


[ieaokeiit Ree ReseNT ING THE 
NORMAEIZEO, BACKSCAIIERED CROSS 
SEC IONSCNBES). 

fhe Sea eeeNUMOER . 2*PI*FREQ/C. 
CHORUS Ole LHe  SRHERE. 

ete sR OUMENT OF THE REGULAR 
GO SLR REGULARS SPHERICAL BESSEL 
FUNCTIONS ANDSTHEIR DERIVATIVES. 
INIECERR OR eintemiIN Value OF KA. 
Pe Sen bertmeiik. JVACUE OF KA. 


KA INCREMENT. 


PEEFFHEE EEE HE FEET HE HELE EE HEEE HEF EFE EEE EEE EF EF HEF HE EF E FEE HHH tee eetee eset teet+tes 


Chapter III. 
1@ 
2Q | 
3 REM 
4Q REM 
5@ REM 
6Q REM 
70 REM 
8 REM 
90 REM 
1@@ REM 
11@ REM 
120 REM 
130 | 
14@ 
150 
160 
170 REAL 
180 
19@ 
200 
210 REAL 
220 
230 REAL 
240 
250 REAL 
260 
270 
280 
29@ = =REAL 
30@ 
31@ Ren 
320 
33Q@ REAL 
340 
350 
360 
370 


a3) 


48@ 
49@ 
422 
41@ 
420 
430 
44@ 
45@ 
460 
47@ 
480 
493@ 
S00 
51@ 
52@ 
530 
54@ 
S20 


[++++ MAIN PROGRAM +4 ttt ttt tt tte EHH tte ete ttttttetetteteteetee¢eet+et 
| 
CALL Bcksctr(Nbs(*) Kamin ,Kamax ,Stp i Game uA eS Thee ae 
| 
CALL Recplt(Nbs(*) Kamin ,Kamax Sto) ! OUTPUT THE NBS). 
I 
PEEEEFEE EEE FEF AEE EEE ETE E EE FEE FEE EE EEE EEF EEE EE HEE HE EH egg eee ttt tteeetetete+tt 
l++++ MAIN PROGRAM CONCLUSION +4+4+4+4+++4t+444444t¢+4¢¢¢t¢t¢et¢t¢tttetteetettet+ets 
PRINTER 15 070" 
PRINT CHRS( 12) 
END 
| 
PEEEEEE EEE FEE FET HEEL FFE FE FEE FEE EEE EEE EEE EEE FE EEE EE FE FH te ttt tgttetttttett test 
| 


IS$SF$SSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSESSSSSSSSSSSSSSSSSSS 
| 


96 


56@ 
57@ 
982 
592 
622 
61 
622 
630 
640 
658 
660 
670 
68 
69 
78@ 
710 
720 
732 
74Q 
75 
76@ 
7708 
78 
790 
B22 
810 
82Q 
830 
842 
852 
862 
87@ 
882 
890 
5@Q 
912 
920 
93 
94@ 
952 
96 
9708 
982 
992 
1082 
1812 
1022 
193 
1040 
1052 
1060 
127 
182 
109 


ISSSSSSSSSSSSSSSSSSSSSSSSSSHSSSSSSSSSSSSSSHSSHSSHSSSSSSSSSSSSSSSSSSESSSSSSESS 
SUB Bcksctr(Nbs(*) Kamin ,Kamax ,Stp) 


Rem 
REM 
REM 
REM 
REM 
| 


FHF E HEHEHE FHE HEHEHE EE FHF FEE EE EH EE HEE EEE Eee tg eget geese gett tee teeteeet et 
itsosOoUBeweneCUERTES Tae magn, NES(*) REPRESENTING THE 
Bote Ten ebmtROooeseCTION NORMALIZEO BY THE CROSS SECTIONAL AREA 


OFo HE sPHERe 
FHEFHEFtHtFFHF HHH tHt HEHE t HEH Et HHH t teat geet gt tttetgttettetetettettt 


t+4+++ MOOULE VARIABLE OECLARATIONS ANO DEFINITIONS ++44++4++4+4++4++4++++7++4 


REAL 


REAL 


REAL 


REAL 


REAL 


REAL 


REAL 


REAL 


REAL 


BEAL 


Lmax 


JCS) 


CaS) > 


Code) 


DyCS|? 


toil) 


Sh) 


| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
! 
{ 
| 
| 
| 
| 
| 
| 
! 
| 
| 
| 
| 
| 
! 
! 
! 
| 
| 


OROER OF THE FUNCTION CALCULATED 
[iecesUEROUTINE.~@at 15 USED 

Ao © LOOP INDEX, (ND OsBESIGNATE 
THE ARRAY ELEMENT FOR THE FUNCTION 
OF CORRESPONDING ORDER. 


Nines ONIEICANT ORDER OF THE SUM 
DETERMINING THE OIFFERENTIAL 
SUMIER ING CROSs SECTION (BSCS ) 
BNOese nm lER ING MODULUS sit)- 


ARRAY FOR THE REGULAR SPHERICAL 
Bese ee rUNC TIONS (RSB Or 
ARGUMENT KA AND OROER @ THROUGH 
LMAX. 


Sen aOR Te wen! YVATIVEoe OF THE 
ob meltiem Oro). 


G@uRtyY POR TIRHE IRREGULAR SPHERICAL 
BESET UNCTIONS (ISB) OF 
ARGUMENT KA AND OROER @ THROUGH 
LMAX. 


CAR er On mine OERIVATIVES OF THE 
odes Ur les |O) ec) op ie 


ARRAY FOR THE LEGENORE 
POLYNOMIALS OF ARGUMENT -1. 


THES REE COMPONENT OF THE 
SCATTERING MNOBGLUS. 


THE IMAGINARY COMPONENT OF THE 
SCATTERING MOOULUS. 


FACTOR COMMON TO RSM ANO ISN. 


Tia eNO Edt tet 

THE ARRAY OIMENSIONS (38+21) NEEO 
POSSE NCRERSED hO ACCOMOOATE 
KAMAX > 3@. 


1182 
111 
Eze 
Piroad 
114 
he 
1160 
117@ 
1182 
1139 
1280 
121 
1228 
123@ 
1242 
[258 
126@ 
1278 
1280 
1290 
1300 
Voie 
1320 
1332 
1342 
1358 
1368 
1370 
138@ 
1330 
1480 
141@ 
1420 
1438 
144@ 
1452 
146@ 
1472 
1480 
149@ 
1502 
151 
ro2u 
153@ 
154@ 


| 
PHEEHEHHFETEFFE HEE EEEE EEE HEE FE FETE EEE EEE EE EEE EEF HEHEHE ttt ttt ttseteteeeeteete ttt 
{ 
1++++ CALCULATING THE NBCS ARRAY +4+tt4+4+t+4tt4tt 44 tt +t 44+ 4tttt¢¢ttetttett+t 
ee UNGER MS | 
INPUT Kamin,Kamax ! DETERMINE DESIREO RANGE Ghai 
Kamin=INT (Kamin) IP VABRUES. 
Kamax=INT(Kamax ) 
Stp=. 1 
FOR Ka=KamintStp TO Kamax STEP Stp 


Lmax=INT(Kat21) 
CALL Drsb(lmax shar Diy( ex) ! CALCULATE THE DRSB OF ARGUMENTS 
! AND ORDERS @ THROUGH LMAX. 
CALL Disb(Lmax ,Ka ,Dy(*),Y(*)) ! CALCULATE THE DISB OF ARGUMENT IS 
| AND ORDERS @ THROUGH LMAX. 
CALL Legt8@(Lmax ,Ka ,P(*)) ! CALCULATE P(*) OF ARGUMENT =i 
! ORDERS @ THROUGH LMAX. 
Rsm=@ ! CALCULATING THE REAL AND IMAGINARY 
Ism=@ | COMPONENTS OF THE SCATTERING 
FOR L=@ TO Lmax | MODULUS. 
Fac=(DjCl )*(24F |) 24P Ceci) #0) (ay Ci) iy eee, 
Rsm=Rsmt(FacetDy(L )) 
Ism=Ismt+(Fac*#Dj(L)) 
NEXT L 
| CALCULATING THE ARRAY REPRESENT 
ines NECS. 
Nbs((Ka-Kamin)/Stp )=4*(RsmeRsmtism*eiIsm)/(Ka#Ka) 
NEXT Ka 
{ 
EEE EF EEF HEE FEEL HEHEHE FETE ELTA FAL HEHEHE LAE ETH FE FE FFF FFE Ett ttt ttt tetetttett+ 
1++++ MODULE CONCLUSION +++ 4444+ 4t 44444444 tt 44 ttt ttt tt ttt ttt ttt tt+ttttett 
SUBEND 
PEEFEFEFEEFEF HEEFT HELE EFE FEF HEF EE EEE EFF E EEF FE EEE HEF Ht ttt ttt tetttttt+e++tttt 
|SSESSS$SSSSSSSSSSSSSSSSSSSSSSSSSSSSSESSSSSSSSSSSSSSSSSSSSSESSSESSSSSSSSSSS 
| 


98 


faa? 
1560 


157@ 
1580 
159@ 
1620 
161@ 
162@ 
163@ 
1640 
165@ 
1660 
167@ 
168@ 
1639 
1782 
i? 1Q 
172@ 
1730 
174 
175@ 
176@ 
177@ 
1782 
ier 38 
1800 
1812 
1820 
183 
1842 
18S 
1860 


| 
ISSSSSSSSSSSSSSSSSSSSSSSSHSSSSSSSSSSSSHSSHSSSSSSSSSSSESSSSSESSESSSSSSESESSSSS 
SUB Drsb(Lma» Ka Djl*#) JC*#)) 

Rien HEFELE EEE EEE FEE EEEEEHEEFEF EA EEE EE EELE ETE EEE TAH HEH tte tt tttttt+t++e t+ et 
REM Dio ODUM neecUlAtes THE DERIVATIVES OF THE REGULAR SPHERICAL 
REM BESSEL FUNCTIONS (ORSB) OF ARGUMENT KA AND ORDERS @ THROUGH LMAX 
REM (ieee eewiin nO Or THe DERIVATIVE INVOLVES THE RSB'S OF THE SAME 
REM AND PREVIOUS ORDER. 

REM FEAEFHEEFEFHEE EEF HEE EEE EE FE FEE EEF HEHEHE EEE EEE FEE EH Ett tte ttett+t+t+e+ttt 
l++++ CALCULATING THE DERIVATIVES OF THE RSB +++4+4+4++4++4+++++t4+t++++t+++++++ 


Gabe Rsbtimax ,Ka,J(*)) ECACCUEATESTHE RSS VALUES J(#). 
ANGee UNI Toe POR DRSe CALC NS. 
CAeeCUEnIe THE - INIT TAL IOR SBS DCO). 
ChECURATE THE REMAINING. DRSE- 


RAD 

wa) =< (0 )tCOS( Kha) )/Ka 

foRee= | TO imax 
Oei==—( CEP Kale ICL IUCR ) 

NEXT L 


Hy 

a 
a 
+ 


| 
PH EEF EEFFHEE FEE E ALE PEE FEE ALF H HEFT TAFEP HEHE THE H+ ett t tet ttt+etttt+tetett+tt++et 
l 
Dee OOUEERCONCMUGION ttttttttttttttt treet eter ttt tet eet ttt ttt ttt ttt + 
SUBEND 
LEAF TAF TELE A FATE LHF EL HEHEHE FETAL EA EEE HEE HEFT HEHE HEHEHE HEHEHE tet ttt ttt estat eettest 
| 
ISSSSSSSSSSSSSSSSSSSSSSSESSSSSSESSSSSSHSSSHSSSSSSSSSSHSSSSSSSSHSSSSSSSSSSSESS 
| 


99 


1872 
188 
1890 
1900 
191 
1920 
1932 
1949 
I dae 
1962 
1972 
198 
1992 
2000 
2010 
2020 
2032 
2040 
205 
2062 
207 
2280 
2292 
2182 
2110 
2120 
2130 
2140 
2030 
2162 
2178 
2180 
219@ 
2202 
2210 
2220 
2230 
2240 
2250 
2260 
lapels 
2280 
2292 
23080 
2312 
2320 
2330 
2348 
2350 
2368 
Zone 
2380 
2598 
2400 
2410 


[SSSSSSSSSSSSSHSSSSSSSSSSSSHSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSS$SSSSSSSSSS 


SUB Rsb(Lmax ,Ka ,J(#)) 


REM SHH HHH HHH tt HHH tt HHH HHH HHH ttt tt ttt tt ttt tte ts etttettt+tt+++ ttt t+ 
REM THIS MODULE CALCULATES THE REGULAR SPHERICAL SBESSseReR tee 

REM (RSB) OF ARGUMENT KA AND ORDERS @ THROUGH LMAX. THe 

REM CALCULATION USES THE CONTINUED FRACTION APPROACH. 

REM FHHtH HHH Ht tH HHH Ht HHH HH tH HH EHH Ht tt tH Ht ttt ttt ttt sttggtttttettt++tt+tet 


f++++ MODULE VARIABLE DECLARATIONS AND DEFINITIONS ++4+4+++4+4+++4+++4++++++4+ 


RENE 


REAL Numer 


REAL Denr 


REAL Scexp 


REAL Nflag 
REAL Dflag 


REAL Bnm 


REAL Binc 


THE EVOLVING RATIO OF J(@)7 Gis 


THE EVOLVING NUMERATOR FACTOR USee 
IN. CALCULATIONS GER as 


THE EVOLVING DENOMINATOR FACTOR 
USED IN CALCULATION Cpe ne 


EVOLVING SCALING EXPONENT FOR VERY 
LARGE OR VERY SMALL VALUES OF See 


@ FLAG FOR NUMR. 


0 FEAG FORS DENS 


EVOLVING TERM USEO SDNe tae 
CALCULATION OF NUMR AND DENR. 


INCREMENT USED IN THE EVOLUTION 
OF BNM. 


PHF HEHEHE E THEE FES HEHE HALE HEHE HEHEHE HE tHe ttt ttt ttt et eetttetettt+et+t++ttt tt 


1++++ CALCULATING THE RSB +++4ttttttt tt ttt ttt te ttt tt ttt ttt tettt+ttet+ + ert 


RAD 


J(Q)=(SIN(Ka))/Ka 


FOR E=i 10 tmax 
R1=1.@ 
Numr=@. 
Denr=@. 
Scerxp=@. 
Nflag=Q. 
Oflag=@. 
Bnm=1.@/Ka 
Binc=2.@0/Ka 


1OQ 


ANGLE UNITS FOR CALC Nae 


*** CALCULATE THE INITIAL **# 
ree RSB, JO). tae 


#*## CALCULATING THE RATIO #+* 
eee (017 3 Cie nae 
INITIALIZING Rice 

INITIALIZING NUMR. 

INITIAL ZING (OER 
INITIALIZING SSGey ee 
INITIALIZING NE EAGe 
INITIALIZING DFEaGe 
INITIAL TZ INGER 

CALCULATE BI itGe 


2420 
2430 
2448 
2450 
2460 
2470 
2480 
2490 
2580 
2518 
2220 
2530 
25408 
2552 
2560 
2570 
2580 
2592 
2600 
261@ 
2620 
2630 
2640 
2652 
2660 
2670 
2680 
2690 
2700 
2710 
2720 
2730 
274 
2750 
2760 
2770 
2780 
2792 
2802 
2810 
2820 
2830 
2842 
2850 
2860 
2870 
2880 
2890 
2920 
2912 
2920 
2930 
2940 


FOR Ordcnt=) 10" 
Bnm=BnmtBinc 
IF Nflag=@ THEN 
Numr=Bnm-Nume 
IF Numr<>@ THEN 
Rl=R1*Numer 
Numr=!.@/Numr 


IF ABS(R1 )>10*250 THEN CALL 
les 
Nflag=1 
ENG) ir 
ELSE 
Nflag=@ 
END IF 
NEXT Ordcnt 


CAC UEATe THE L “NAKED * 
TaN MC IORS OF JC G97 ICL). * 
INCREMENT BNM. 

CHECK NUMR=0 FLAG NOT SET; 
UPDATE NUMR IF NUMR<>O. 


UPDATE RL IF UPDATED NUMR<O@. 
PREPARE NUMR FOR NEXT EVOLUTION. 
CHECK SCABINGOF UPDATED RL. 


scale(R] ,Scexp) 


[Eenmeno Oe PREV TOUS STEPS “ARE 
DEFERRED FOR ONE INCREMENT OF BNM. 


* CALCULATE THE REMAINING NUMR * 

STONOVDENR TERMS GE d(GI7I(L). * 
NOTE WHEN NUMR=DENR THEN 
RL=J(O)/ICL). 


WHILE Numr<>Denr OR Nflag=! OR Dflag=! 


Bnm=Bnm+Binc 


IF Nflag=@ THEN 
Numr=Bnm-Numr 
IF Numr<>@ THEN 

Rl=R1*#Numr 
Numr=1.8/Numr 


PewnoscR ly 210°250 THEN CALL 
eles 
Nf lag=! 
END IF 
EESE 
Nf lag=Q 
END IF 


IF Dflag=@ THEN 
Denr=Bnm-Denr 
IF Denr<>@ THEN 

R1l=R1/Denr 
Denr=1.@/Denr 


fee oottl te c50 THEN CALL 
ELSE 
Oflag=I 
END IF 
ELSE 
Df lag=@ 
END IF 


END WHILE 


10] 


INCREMENT BNM. 


CHECK NUMR=@ FLAG NOT SET; 
UPDATE NUMR IF NUMR<>O. 


UPDATE RL IF UPDATED NUMR<>®@. 
PREPARE NUMR FOR NEXT EVOLUTION. 
CHECK] SCALING OF UPDATED RL: 


Scale(R1l ,Scexp) 


Caieena. Oe Hearne VIOUS STEPS ARE 
DEFERRED FOR ONE INCREMENT OF BNM. 


CHECK DENR=@ FLAG NOT SET; 
UPDATE DENR IF DENR“D®@. 


Ue OAlesRES Tr UPDATED DENR<7@. 
PREPARE DENR FOR NEXT EVOLUTION. 
CHECK SCALING OF UPDATED RL. 


Scale(R] ,Scexp) 


° 
5 


PePOGEno 7 wine PREVTOUS STEPS ARE 
DEFERRED FOR ONE INCREMENT OF BNM. 


Zda0 
2960 


29370 
2382 
2990 
3000 
3010 
3820 
30350 
3042 
3050 
5@60 
3070 
3080 
309@ 
3102 
3118 


J(LISICOI/ CRI #10 Gseexn 7 eee J(L) DETERMINER 


NEXT L 
| 
PHAEHEEHE EEE FE FAFEFAFFE FFF H EEE H FHT HATH FEE HE HHH HH ttt ttt t+ tt tt+ttet+et++tt++tttt 
| 
[++++ MODULE CONCLUSION ++4+4+t++t+t++4tt +444 ttt ttt ttt ttttettttettettt++t++++t+ 
| 
SUBEND 
| 
LHHEEHEEHHAFEFEF FALE FFA AFA FE HEATH FATEH HATH HEHEHE H+ 4+ 44+ tt tetttt++ttt+++++t¢+t++tt 
! 


| SSSS$SSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSS 
| 


Oz 


3120 
3130 
314@ 
3150 
3160 
a1 70 
3180 
3190 
3220 
3210 
3220 
3230 
324 
3250 
3260 
3270 
328 
3290 
3302 
3310 
2520 
53332 
3340 
3352 
336 
35 70 
3380 
3390 
3400 
3410 
3420 


| 
(SSSSSSSSHSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSHSSSSSSSSSSSSSSSSSSSSSSSSHSSSS 
SUB Scale(Rl ,Scexp) 
REM LEFEFFEE EEE EHEEFEF EEE EEE EEE HEE EEE HEHEHE eet te teeteeteteeeet+ettettett ttt 
REM THIS MODULE SCALES RL TO MAINTAIN NUMERICAL ACCURACY IN THE 
REM ONGOING CALCULATION OF RL. 
REM A LELEEEHEEE EEF EEE EE EE EHEFE LAE EEE EEE HEF EEE EEE He eee tte teteteteetttett 
| 
l++++ SCALING RL ttt ttt Fett ttt FFF ttt ttt t ttt ttt ett ggeetetetetttttttt++t+tt 
teeRID18° 2590 THEN 
R1=R1+*#10*(-250) 
Scexp=Scexp+250 
END IF 
IF R1<10*%(-250) THEN 
R1=R1+*10°250 
Scexp#Scexp-250 
END IF 
| 
| 
LHEEEE EEE EEE EAL HEE FE FEE ELE EEE EE EEE FEE EEE EEE FEF ETE FE +H e+ es teeeete+e+eeses 
(++++ MODULE CONCLUSION +4+++++4t+t+ ttt ttt ttt tt ett ttt tegttettttttetet+etettet 
{ 
SUBEND 
| 
PEAEEEE HEE EEE LEAF EEE FEF ELE EEA HEE FEE FEE EEE FEE HFEF HFEF +44 ett eteteteeeet++et+et 
| 
ISSSSSSSSSSSSSSHSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSHSSSSSSSSSS 
| 


lO3 


3430 
3440 
3450 
3460 
3470 
3480 
3490 
3500 
SS) 
3520 
395¥ 
3540 
Joa0 
3560 
3570 
358 


S590. 


3609 
3610 
3620 
3630 
3640 
3650 
3660 
5670 
3680 
363 
3700 
S119 
3720 
3730 
3740 


1SSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSESS 
SUB Disb(Lmax ,Ka ,Dy(+#) ,Y(*#)) 

REM FELHHEFE HEE EE EEE EEF HEE EE ALE HEF HEE EE EEE HEH tte tgt ttt tetett¢tttt++etteetetts 
REM THIS MODULE CALCULATES THE DERIVATIVES OF THE IRREGUEAR SSPE ee 
REM BESSEL FUNCTIONS OF ARGUMENT KA AND ORDERS @ THROUGH LMAX. THE 
REM CALCULATION OF THE DERIVATIVE USES THE ISB’S GFeiHes shies 

REM PREVIOUS ORBER- 

REM FHETHE HEH HH HEHEHE ttt HE HHH ett tt tt Het ttt tte tt sett tte eter tte tttt+ttttt+t+ 
| 

l+t+t+t CALCULATING THE DERIVATIVES OF THE ISB +++++4++++++++++7+7777 Pee 


CALL Isb(Lmax ,Ka,Y(*)) CALCULATE THE ISB VabWESe 
ANGLE UNITS FOR CAG shoe 

THE INITIAL ISB DERIVATIVE 
THE REMAINING ISB DERIVATIVES. 


RAD 
Dy(@)=(-Y¥(0)+SIN(Ka))/Ka 
FOR L=!] TO Lmax 
OyCE y==( CET! Kade ye +y Cle] 
NEXT L 


| 
PE EELTE FEE EEE EEF EE EEE EET EEE EEE HE ALE E EEE EEF EE EE EHH Ht te ttt get ttteeettttttt 
| 
l++++ MODULE CONCLUSION ttttt ttt tt ttt ttt tte ttt tte tte et tt ge ttte ttt etetetttt 
| 
SUBEND 
| 
PEEEEEEE EEE EE EE EEE EE EAE AEE EE EEL AL EEE HELE FEE EEE PETES HEE TET Et Het tee ettt+itt 
ISS$SSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSISSSSSS 
| 
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3750 
3760 
3770 
378 
3790 
3800 
3812 
3820 
3830 
3840 
3850 
3860 
3872 
3880 
3890 
3900 
Bait Q 
392 
3930 
3940 
6450 
3960 
52/0 
3980 
S240 
4902 
4010 
4@20 
4030 
4042 


l 
ISSSSSSSSSSSSSSSSSSSSSSSSHSSSSESSSSHSHSSSSSSSSHSSSSSSSSSSSSSSESSSHESHSSSSSSSSSS 
SUB Isb(Lmax ,Ka,Y(*)) 

REM —- FF HHH HEHEHE THF HTH t HHH tHE HEHEHE ttt ett g tt tet gett ttettetettt+ tt tttet ttt 
REM il owe bUPENCAeeCUlLAIES THE TRREGULAR SPHERICAL BESSEL FUNCTIONS 
REM OF ARGUMENT KA AND ORDER L THROUGH LMAX. THE CALCULATION IS 

REM ACCOMPLISHED BY FORWARD RECURSION. 

REM FEFHETTHEFHEFFEFE EHF EEE FLEE FEF EEF E EET ETE HEHEHE ttt Hee eegeteetttt++t++ett 


| 
'++++ CALCULATING THE ISB +++4t+++++++ 4+ t+ tt ttt t+ +++ ttt tett tt t+ttt+tt++t+t+++ 


| 
RAD l ANGEE UNITS FOR CALC'NS . 
Y¥(@)=-COS(Ka)/Ka (PIN ET TAS 1S68 S CALCULATED. 
Y(t )=CY(QI-SIN( Ka) )/Ka 


FOR t=Z2 TO Lmax Sire Matin iNnGs)Sse’ S CALC’ OD. 
gle =e, tL lo/Kka)* 1 Cle) —y Cl =2 >) 
NEAT L 


PHAFHE HEHEHE FE FH FEAF FF EAE HEHE HEHEHE HEHE E+E +H tHe ttt tt tt ttt ttt tt tt+tet+tet+ tee t+et 
l++++ MODULE CONCLUSION +4++++++++++++ +t t+ t+ t+ tt ttt +t t+t+ttt+tt+e++ttete++tt 
] 
SUBEND 

| 
PHHFE HFEF ELH FEAF EE HATA AEE L HEE FHF EEE HE tH Fee tt tette tt tte tttttt tee ett ttt ttt 
ISSSSSSSSSSSSSSSSSSSSSSSSSESSESSSSSSSSSSESSSSSSESSSSSSSSSSSSSSSSSSSESSSSSSSSF 
| 


Os, 


4050 
4062 
4270 
4082 
4092 
4100 
411@ 
4120 
4152 
414@ 
4152 
4160 
417@ 
4180 
419@ 
4200 
4210 
4220 
4230 
4242 
4252 
426@ 
4272 
4280 
429@ 
4322 
4310 
432@ 
4330 


| 
ISSSSSSSSSSSSSSSSSSSSSSESSSSSESSESHSSSSSSSSSSSSSSSSSSSSESSSSSSSSESSESSSSSSSSHS 
SUB Leg!i8@(Lmax ,Ka,P(*)) 

REM LEFF EEE E EE EEE EE EET ETE ETE EH PEP Pat tt Feet ee 
REM THIS MODULE CALCULATES THE LEGENDRE POLYNOMIALS OF ARGUMENT =! ANG 
REM ORDER @ THROUGH LMAX. CALCULATION IS ACCOMPLISHED BY FORWARD 

REM RECURSION. 

piled FHEEEEEE HFEF HAAFE EEA AEE E FLEA LEE EEE FEF E FET FEE H Ett t ttt t tt tttettt+ettt 
| 

l++++ CALCULATING THE LEGENDRE POLYNOMIALS ++4+++++4444+344 44 54 5 7 se 


2h ! ARGUMENT OF THE LEGS PGi 
j 
Foe =1 ! INITIAL LEG. POR SSeS Cale aan 
EA lee x 
FOR E=2 10 "imax | REMAINING LEG. POL’ S) CaBG Es 
PCL peC2=C1 7b be A EP UL= 1 SCS ee 
NEXT FL 


| 
LEE EEEF EEE AE HEE HEEEE ELE EE EEE EE EE EEF EEE EEF EEF HEE A FEF FF egg gtggtgttgtggtt+tet 
| 
l+4+44 MODULE CONCLUSION t+ t+tt4t +++ ttt tt ttt Ett Ht tt ttttttt¢ tt ttttet+t+++ett 
{ 
SUBEND 
| 
LEEEEEEE EAE HEHE FEE EEE EFEF EEE EF EFEEE EEE ETE EEE EEE FEE AE HEFE FFA tt eet ttettt++ test 
| : 
ISSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSESHSSSSSSSSSSSSSSSSSSESSSSSSSSSS 
| 
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4340 
4352 
4360 
437 
4380 
4390 
4400 
4410 
4420 
4430 
4440 
4450 
4462 
4470 
4480 
4490 
4502 
4510 
4520 
4530 
4540 
4550 
4560 
4570 
4582 
4532 
4602 
4612 
462@ 
4630 
4640 
4652 
4660 
4670 
4680 
4690 
4700 
4712 
4720 
4730 
4740 
4750 
4762 
4770 
4780 
4730 
4800 
481@ 
4820 
4850 
4840 
4850 
4860 
4870 
4860 


l 
ISSSSHSHSSSSSSSSSSSSSSSSSSSSSSHSSSSSSSSSSSSHSSHSHFSSSSSSHSSSFESSSSSESSSSESSSSSSESS 


SUB Recplt(Nbs(*) Kamin ,Kamax ,Stp) 
REM FEF EEFE EEA EEE EEF HEEL EEE FHF EEE FEE EEE Et tt ttt ttetettetetttt+tset+tt+tt 
REM THIS MODULE MAPS THE CALCULATED SCATTERING MODULUS ARRAY SM(*) ON 


REM fi SEMISIEOG PLOT: 
REM FHEFFFEFF HATE FF TH EEE FEE FEF ETH EET HEHEHE tte ttt ts ttt ttt tte t+ttttt+ttet 


| 
Bee TOUULE VARIABLE DECEARATIONS AND DEFINITIONS t+++++4++4444+44+444+4+4 
| 


Fae Manu apo lash VALUE “IN 
GRAPH TG omeimorlLiAy Und TS. 


REAL Xgumax 


THE MAXIMUM ORDINATE VALUE IN 
GRApihes SUPSriEny aUNT TS. 


| 
| 
| 
REAL Yguman | 
| 
| 


| 

| oh oapab bape SP Sr te o 8 ap aes s qeees page aya sme noe sae eee eae oa ee ye eaemn sneeh Sen om senna em oh Si Sem de oi Me ote oh Se ota oe or or | 

| 

See eae Gel HeePeO) eR tttt+tttet ttt ett tees ete tt eee ttt tet ttt 
| 

GINIT 

peOtTER IS CRT, INTERNAL” 

GRAPHICS ON 

GCLEAR 

ee et tet tthe tat tt ttt te eet ett ttt tt ett tt 

| 

met ec ONo IeUet LON HNO SERGELLING +4tttetttttttettt4tt+tt+tt+++t+++44+ 
| 
| 

Xgumax=100*MAX( 1 RATIO) 

Ygumax=10@0*MAX(1,1/RATIO) 


Gal7e— 4 
Ne TYPE 1 


LORG 4 te nee UNG aiGfeEenBoelSons ax ls *«* 
MOVE Xgumax/Z ,O 
BABEL “Ke” 


DEG ls * SPABELEING THE OROINATE AXIS *+# 
EBIR 9 

EOKRG 6 

MOVE @,Yqumax/2 

meee NORM DO BCKSCTRG CROSS SEC’N*” 

LDIR @ 

** GRAPH CONSTRUCTION ** 

DEPT Nessie. OTA PEOTTING AREA 
IN GRAPH UNITS. 


VIEWPORT 18,Xgqumax-S ,1@,Yqumax-S 


WINDOW Kamin ,Kamax ,@,1.2 
miv= |YPE 4 
oe lO Stp*i0,.1 ,Kamin,®,61 .10,.2 


1Q7 


4839 
49302 
4912 
4320 
4930 
4942 
4950 
49360 
4972 
49380 
49390 
5000 
5010 
5020 
5030 
5042 
SQ50 
5062 
5070 
5080 
5092 
5100 
Slee 
5122 
S138 
5142 
5150 
S160 
Sig 
5180 
5130 
5202 
5210 
S220 
5230 
5242 
5252 
526 
S270 
5280 
529 
5382 
53.10 
5520 


Clair OFF 
Gale) 
LT NE yee 


LORG 6 | *«* LABELLING THE ABSCISSA SCAlENNS 
FOR IT=Kamin TO Kamax STEP Stp*2Q 

MOVE I ,@ 

CABEL =! 
NEXT I 


LORG 6 | ** LABELLING THE ORDINATE SCAR@Eaas 
FOR] f=0/ TO Ser ss 

MOVE Kamin,!I 

Lyayste|e Il 
Nei 


Glir ON 
| 
LEE EEFEFTHEFEEFEEEE ETE EEF EEE FFF FEE EE PETE E PEPE EEE te ee ttt ee ee 
[++++ PLOTTING THE CALCULATED NBCS +++4+++++++4++44+44++44+4++4+444++47747 2 
sma 
FOR I=KamintStp TO Kamax STEP Stp 
PEOT 1 Nest (i-Kemined7scip) 
NEXT I 
| 
PEEEFEEEE EP EEER EEE FETE EEE FETE EEF EE HEE Pat at Ee ett a 
l++++ PRINTING THE PLOT t+ t+44+t4tt4t+4 +++ 44444 FF 4444+ 444+ +444 +++ +444 +4+4+444 
| 
DUMP GRAPHICS #701 
| 
LEEEEEEEELEF EEE EEE EEE EE EEE EEELE FEF ELE FEF EEE EHF Fett tttt ttt tett++4++44rrst 
! 
l++++ MODULE CONCLUSION +4+4+444+ 44 +444 +4444 44444444444 4+4+4+++4+4+4+4+4++4+4+4+4+4+4+4++4 
| 
GRArH CS OFF 
SUBEND 
| 
PEELE EEE HEEF ELE TL AFA THE FEF TEFL FEE EE FLEE ALE FTE FEE Ft tte +t ttt 4444+ FFE tee 


ISSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSESSSSSSSSSSSSSSSSSSSSSS 
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APPENDIX D 


ACOUSTIC ARRAY EFFICIENCY PROGRAM 


This program was used to analyze data and determine the calibration 


product E_£,. The output of this program was also used to construct 


Figures 13, 14, and 15. 


The power supply mentioned in the Gain30 and Gain35 subprograms 


was the HP3314A function generator. All data contained in the data 


subprograms are repeated in Appendix E. Data units are available in 


Appendix E with their respective data. 


1@ 
20 


4 
~ 


4Q 
5Q 
62 
7Q 


92 

180 
13@ 
14Q 
15@ 
162 
170 
18 
192 
222 
21 
220 
230 
240 
200 


ISSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSESSSSSSSSSSSSSSSSSSSSSSSSESSS 


| 

REM 
REM 
REM 
REM 
REM 
REM 
REM 
REM 


REM 
| 


S4t Hts gt tg ps Hs tots sss ss sss pss sss ggg stg seg est eee gssgt tte ttt+see 
HeOuatlhGeaRnkay EFF ICLENCY 


fp ROSRal CARCULAITES THE EFFICDENCY AND CENTERLINE GAIN OF 
GNeACOUSTiCeaRmny Bosco ON TRE RETURN FROM AN ACOUSTICALLY HARD 
Se ee eon Oa eiEnmO® sENSGONIFIED BY AN INCIDENT PLANE WAVE OF 
ERE QUENCY Ee 


FHH+EFFEFF HEE LELEAFHALE ELLE HEEFT TELE F LETHE FEF FEL ELF EF HHH H+ +H H ++ ttt tet 


eve hibee HeGeomalloNS AND MWEFINITIONS ++++++4+++4+4+++4+4+4+4+++44+4++++++ 


. 


REAL Vesr MEANS SOURRE VOLTAGE PECETVED EROM 
tle leecel sore ne eAND MEASURED 
Or UE Ree Re=mrr ele] CAlEON . 

REAL Vmsn MEANS SOUARE VOLTAGE RECEIVED IN 


THESABSCENCe OF A TAaRGE! AND 


l 
] 
| 
l 
| 
PVea over len FReE-oOnMr , 
l 
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262 
270 
280 
290 
300 
510 
320 
330 
34 
350 
362 
570 
582 
390 
400 
41@ 
422 
430 
440 
452 
460 
472 
482 
492 
520 
510 
520 
530 
542 
aau) 
560 
570 
580 
530 
500 
612 
620 
630 
642 
650 
662 
672 
680 
592 
700 
710 
720 
730 
740 
750 
762 
772 
780 


REAL 


REAL 


REAL 


REAL 


Rene 


REge 


REAL 


Regie 


REAL 


REALE 


REAL 


REAL 


REAL 


REAL 


| omen: j ‘el 
mm es 
tee * fhe 


Vmst 


ee 


Go 


ELEO 


At 


C.Cmorsd 


Efreaq 


Ka 


Aa 


Gpre 


Nbs 


Aitot 


Sasph 
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MEAN SQUARE VOLTAGE TRANSMITTED. 


THE RANGE FROM THE PLANE OF Sie 
ARRAY SPEAKER DIAPHRAGMS 10™Tae 
CENTER OF THE TARGET SPHERE AND 10 
THE SPECTRUM ANALYZER MICROPHONE. 


THE EFFICIENCY OF CONVERSION 
FROM TRANSMITTED ELECTRICAL POWES 
TO TRANSMITTED ACOUSTICAL POWERS 


THE CENTERLINE GAIN OF THE ACOWSiame 
ARRAY. 


ET*G0 


THE ATTENUATION COEFF ICT EIi@® 
SiGe fous anaes 


THE SPEED OF SOUND IN HUNTON aie 


THE SELECTED FREQUENCY OF Sine 
ACOUSTIC ARRAYS 


THE WAVE NUMBER, 2*P1*ERREQ 7 ee 
THE RADIUS OF THE TARGET SPHGeiae 


K*A, THE ARGUMENT OF THE REGULAR 
AND IRREGULAR SPHERICAL BESSem 
FUNCTIONS AND THEIR DERIVATIVES. 


THe APERATURE AREA OF THE AR eee 
THE PRE-AMP’S ELECTRICAL Game 


THE VARIABLE REPRESENTING THE 
NORMALIZED, BACKSCAT TERED  GRase 
SECTIONS NEGo 7] 


CROSS SECTIONAL AREA OR Sits 
TARGET] SSP Benes 


SURFACE AREA OF AN IMAG] NARs 
SPHERE CENTERED ON THE ine 
AND HAVING A RADIUS EQUAL TO 
THE RANGE TO THE CENTER OF SiGe 
TARGET, Serene e 


VED TReeal 
OS Tee 
ICAL POWERS 


THE RATIO G5ee=e = 
ELECTRICAL FPaNe aa 
TRANSMIT IED Veber. 


790 
820 
810 
820 
830 
840 
350 
860 
870 
880 
890 
920 
910 
920 
930 
940 
950 
960 
970 
980 
990 
1020 
1010 
1020 
1030 
1040 
1050 
1060 
1070 
1080 
1090 
1100 
1112 
1120 
1130 
1142 
1152 
1160 
1170 
1182 
1190 
1200 
1210 
1220 
1230 


fee cr PPE seer eC lENCY) OF S@@NVERSION 

| FROM RECEIVED ACOUSTICAL POWER 
{ 

| 


TO RECEIVED ELECTRICAL POWER. 
ee ee ee ee eee 
| 
b++++ MAIN PROGRAM ttt ttttt HHH tHF HE ttt ttt tet tt tte tttetetttt+ettetet++tt+ 
! 
y 
CALL Init(VUmsr ,Vmsn,Vmst ,R,Etgo,At ,A,Ka,Aa,Gpre) 
| INPUT ANDO CALCULATE REQUIRED 
! PARAMETERS. 
l 
CaLlL BeksetrckKea ,Nbs? Miaecurnlo. Hey NBCS. 
PRINT “NBS=";Nbs 
| 
Atgt=PI+#A#A 
PRINT “ATGT=";Atgt 
Sasph=4+PI*R+#R 
Prat=(Umsr-Umsn)/(Gpre+Gpre#Vmst ) 
Er=(Prat )*(Sasph*Sasph/Etgo )#*(EXP(Z2*At#R )/(AatAtgt ))/Nbs 
| 
aren ee fer 
PRUNdesEtGo= se1G0 
PRINT “Go=";Etgo/Er 
Ph Cee Pee see ter 
! 
| 
PHEEEEE HEE EEE EEE FFE EHF HEE LEH FEFE FE FE FHE FEE FFF EEE HEHEHE Ht tee tte tetttt+ettttset+ 
l++++ MAIN PROGRAM CONCLUSION ++t+t+tt+ttttttt +++ ++ ttt tt tttttttt+ttttttt++tt 
PRINT CHRS( 12) 


BrONIeER {3 | 

PRINT 

PRINT "END" (MAbvices The USER OF THE THAIN 
PRINT PROGRAMS CONCLUSION: 

END 


| 
PALFFLFLE EHH LHE LEFF FEF HEHE FFE FFFFFLE FEF HE LL HELE FH EEE HEHEHE HEHEHE Het tte tt tttetee 
(SS$SSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSFSSSSSSSSSSS 
! 
| 


11] 


124@ 
1250 
1260 
1270 
1282 
129 
130@ 
131@ 
1320 
133@ 
134@ 
135@ 
136@ 
1370 
1382 
pers 
140@ 
141@ 
142@ 
145@ 
144@ 
1450 
146@ 
1472 
148@ 
149@ 
150@ 
Se 
1520 
153@ 
154@ 
is)s) 
156@ 
157@ 
158@ 
159@ 
160@ 
161@ 
162@ 
163@ 
164@ 
1652 
166@ 
167@ 
168@ 
169@ 
172@ 
171@ 
1720 
1730 
174@ 
1752 
176@ 
177@ 
1782 


ISSSSSSSSSSSSSHSSSSSSSSSSSSSSSSHSSHSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSS 
SUB Init(VUmsr ,Vmsn ,VUmst ,R,Etgo ,At ,A,Ka,Aa,Gpre) 

FEFFFHE THEE FEA FEE FEE EHF EF HEF HE ELE FE FEF E THEE EEE EHH He ee eteteeeeeeetesestet 
THIS MODULE DETERMINES THE VALUES OF THE PARAMETERS REQUIRED FOR 


REM 
REM 
REM 


REM 
! 


ltee+ MODUCE VARI ABEE PECEOR MT ONorn 


REAL 


REAL 


Rete 


REAL 


REAL 


REAL 


Rea 


REAL 


REAL 


REAL 


REAL 


REAL 


REae 


REAL 


REAL 


THE GAECUE ATION - OF sen. 


FHEFFEFFHFEFEFEAFE EE FFT HEFT THEEF FEET E+E FEE + +++ ttt tetegtetetetetettt+tt 


Ch 


Vr(6Q) 


J,dJt,Jn,Jtn 


ie 


Patm 


Rh 


erin 


Dia 


Nuspkr 


Diaspkr 


Vt(60) 


Noise(6Q@) 


Tnoise(6@) 


Ten 


Rhn 


le 


ND DEFINITIONS +++++++++++++4++++++4++ 


A CONTROL VARIABLE USEGSEGE 
CONTROLLING SELECTION OF THE 
CORRECT GAIN UAlAVSETse 


THE ARRAY REPRESENTING THE 
MEASURED TARGET RETURN VOLTAGES. 


THE NUMBER OF ARRAY ELEMENTS TORR 
ONE AND A LOOP CONTROL VARIABLES 


THE AMBIENT TEMPERATURE IN OEGeaiEe 


THE ATMOSPHERIC PRESSURE MEASURES 
LN MILEIBARS. 


THE RELATIVE HUMIDITY MEASURE Dae 
PERCENT. 


THE MEASURED TIME FOR THE ACOUSIIR 
SIGNAL TO MAKE THE ROUND TRIP 
BETWEEN THE ARRAY AND THE TARGET. 


THE MEASURED DIAMETER: Ost 
TARGET Sener e 


THE NUMBER OF SPEAKERS I Niiige 
aunrlany 


THE AVERAGE DIAMETER OF AN ARRAY 
SPEAKER. 


THE ARRAY REPRESENTING Sie 
MEASURED TRANSMISSION VOETAGESs 


THE ARRAY REPRESENTING THE 
MEASURED NOISE RETURN VOLTAGES. 


THE ARRAY REPRESENTING Tie 
MEASURED TRANSMISSION VOLTAGES SGe 
NOISE DETERMINATION: 

Tc FOR NOISE DETERMINATION. 


Rh FOR NOISE DETERMINATION. 


179 
1800 
181 
1820 
1830 
1840 
1850 
1862 
1872 
188 
1890 
1902 
191@ 
1920 
193@ 
194@ 
195@ 
196@ 
197@ 
198 
1990 
2000 
2012 
2820 
2030 
2042 
2058 
2060 
2070 
2082 
2890 
2180 
21102 
2120 
2130 
2142 
Zia’ 
2162 
2170 
218 
2192 
2280 
2210 
2220 
ZL50 
2240 
2250 
2260 
2270 
2280 
ZEAQ 


fone Faun Patm FOR NOISE DETERMINATION. 


REAL Efreaqn Efreqn=Efreq. 


REAL Cn C FOR NOISE DETERMINATION. 
REAL Umstn Umst FOR NOISE DETERMINATION. 
REAL Rspkr THE AVERAGE RADIUS OF AN ARRAY 


| 

1 

| 

J 

{ 

REAL Atn | At FOR NOISE DETERMINATION. 

{ 

| 

| 

{ 

| SPEAKER. 

! 


| 
LEEEEEEEEEEFEFEFEEEE EEE FEE FEE EEL EEE EEF EEE EEE EEE EEE EF EF EEE EEE HE HH tte tte e tet 
| A 
l++++ DETERMINING SYSTEM PARAMETERS +4++++++4+4++4+ +444 +t+t+44+444+4++ 44+ 4444444 
Se INTER 1S | 
arnt WHICH TAREE! DATA SEY “; 
Pein’) OO YOU WISH ANALYZED?” 
PEeINTER 1S 701 
INPUT Sphere 
PRINT “TARGET NO. "Sphere 
PRINT 


IF Sphere=! THEN 
Ch=3.@ 
CALE Gaino(Ch ,Etgo?) 
CALL Spherel(Vr(#),J,Te,Patm,Rh,Tottm Oia ,Efreg ,Nuspkr ,Oiaspkr ,Gpre) 
GALE TranmssO( Jt Vic *)) 
CALL Noise3@(Jn,Jtn,Noise(*) ,Tnoise(*),Ton,Rhn ,Patmn ,Efreqn) 
END IF 
| 
IF Sphere=2 THEN 
Ch=3.5 
Coie cainocch,Etoo ) 
CALL Sphere2(VUr(*),J,Tc ,Patm,Rh,Tottm,Dia ,Efreg ,Nuspkr ,Diaspkr ,Gpre) 
Paria camsso( Jt. Vit *)) 
CALL Notse3S(Jn,Jtn,Noise(*) Tnoise(*),Ton,Rhn ,Patmn ,Efreqn) 
ENO IF 
i 
IF Sphere=3 THEN 
Cn=4.5 
Shee ecaimnovwch Etgo) 
CALL Sphere3(Vr( +) J,Tce,Patm Rh,Tottm Dia ,Efreq ,Nuspkr ,Diaspkr ,Gpre ) 
Paee transa5(Jt Vic *)) 
CALL Noise35(Jn Jtn,Noise(*) ,Tnoise(*),Tcon,Rhn,Patmn ,Efreqn) 
ENO IF 
i 


yo 


2300 
2310 
(ici) 
233 
2342 
2352 
2560 
2570 
2530 
2390 
2400 
2410 
2422 
24350 
2442 
2450 
2460 
2470 
2480 
2492 
2500 
251 
2520 
2532 
2540 
230 
2562 
2570 
2580 
2550 
2602 
2610 
2620 
2630 
2642 
2650 
2660 
2672 
2680 
2020 
2722 
2710 


IF Sphere=4 THEN 
Cha3.5 
CALL GainotGn Etea. 
CALL Sphere4(VUr(+),J,Toc ,Patm,Rh,Tottm,Dia ,Efreq ,Nuspkr ,OlaspkKr ope 
CALL. Irans3SG@iie vt ts) 
CALL Noise3S5(Jn,Jtn,Noise(*) ,Tnoise(#),Teon,Rhn ,Patmn Efregn) 
ENO @ie@ 
| 
CALL Voltsms(VUr(#), J ,Umsr) 
PRINT “VUMSR=";Umsr 
CALL Attent len hha Patane Ga otne et requ) 
CALL Voltsms(Tnoise(*#) Jtn,VUmstn) 
CALL Voltsms(Noise(#) Jn, ,Umsn) 
CALL Voltsms(Vt(#) Jt ,Umst) 
PRINT. “UMST= 3Umst - 


A=Dia/2 
PRINT “A=";A 
CALL Range(lattm Ale Rn Petm. © At citroen 
PRINT SCs 
PRIN A= Fat 
PRINT “ATTEN="s3EXP(At#R ) 
PRINT "“R=";R 
Umsn=VUmsn*EXP(2*#Re(Atn-At ))*Umst/VUmstn 
PRINT “VUMSN="3;Umsn 
Ka=2*PI#Efreq#A/C 
PRINT “KA=";Ka 
Rspkr=Diaspkr/2 
Aa=NuspkrtPleRsokreRspkr 
PRINT “AA=";Aa 

| 
PHAHHHHEHEHE HEHEHE HEHE HEHEHE HEHEHE HELE EHH EHH HEHE FHF E FFE EEE FEF HE ++HE HEH HH H+ HHH 
| 
l++++ MODULE CONCLUSION ++ +t ttt tte tHHHHt HHH HHHHHH HEHEHE HEHEHE HHH ee ttt e+e ttt 

| 
SUBEND 
PHRF HFEF ELE HHA FHF LFF HE HHH HEE EEE HH EFL FFE HAE FEE HEHEHE HFEF HE +4 F¢ 4444444444 
| 
ISSSSSSSSSSSSSSHSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSS 
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2720 
21350 
2740 
2750 
276@ 
2770 
2780 
2790 
2800 
2812 
2820 
2830 
2840 
2852 
2862 
2870 
2880 
2890 
2900 
2210 
2320 
2350 
2940 
2950 
2962 
2970 
2980 
2998 
3@0@ 
3@1@ 
3020 
3030 
3040 
305@ 
3060 
3078 
3@80 
303 
31@@ 
3112 
3120 
21:50 
3142 
3150 


| 
ISSSSSSSSSSSSSSSSSSSSSSSSSSSSS$SSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSS 
Sie Gainot Ch, Etgo) 
REM LHELEFEFHFHEFE HEHEHE FEL HEE FEF HEEFT EEE HEE FHEE HEF EET ETH E+ +e Hee tee teetettttt 
el lalownOPUle BeleRMUNEs HE VALUE OF Ettbo. 
REM LEFHEEEFE FEET FEEL FE FEE HFEF EEE FETE FETE AHH THF Het teetteteettetettet ttt 
| 


I++++ MODULE VARIABLE DECLARATIONS AND DEFINITIONS ++4+4++++++4++++++++++++ 
i 


REAL Vrms ' ROOT MEAN SQUARE VOLTAGE MEASURED 
! IN GAIN DETERMINATION. 
| 
REAL Ir ISMEASCURED INTENSITY Al R; 
le Oty Hot. 
| 
REAL Za | ACOUSTIC ARRAY IMPEDANCE. 


PEELE EHTHEHE HEHEHE HEE HEHE EHH EFA LEE EHEEFE HEHEHE HEF EEE HH tte ete e tte ttt tteetettett 
l++++ DETERMINATION OF EtGo t++ttttttttettetttetttt tt te tet tttttetetttett+tetet++ 
| 

[& Ch=3.@ THEN CALL Gain3@(Vrms ,fc,Rh,R ,Patm,Ir ,Efreaq,Za) 
Peeehn=3-5 THEN CALL Gain3s5(Vems ,Tc RAR ,Patm.Ir ,Efrea.Za) 
PObtaNINGevme JES FROM APPROPRIATE 
Oe Ute. 
| 
CALE Atten(Te .Rh,Patm,C At ,Efrea ) ISGReCeeo INGA I. 
Pt=VUrms+Vrms/Za 
Brago-ir*4*Ple*R*eR*EXPCAt#R )/Pt Pe CHeCULATING JEtGo-. 
| 
PEE E HEHEHE LH EEAEEF HLH EEE HEHEHE FH HEHE HLF HEHEHE EE HEHE THEFT HEHE HF tt tt te tte ttt+ 
l++++ MODULE CONCLUSION 4+4+4+++4++4+ +++ 44+ +++ + +4++tt+tetettteettttttet+ettett+ 
| 
SUBEND 
| 
PHF EE FEE HEE HEHEHE FEET HEHEHE EEE HE FEET HEE HEF EH HEHE FH HEF HEFT ET HEF H tt tt ete te tteteet+ 
ISSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSHSS 


|} 


3162 
3172 
3182 
6132 
3202 
3212 
3220 
3230 
3242 
3258 
3260 
3270 
3280 
3290 
33200 
3310 
3320 
3330 
3342 
3350 
3362 
3378 
3380 
$450 
34200 
3412 
3420 
3430 
3442 
3450 
3462 
3470 
3480 
3490 
3520 
5512 
$520 
3550 
3542 
S450 
3560 
3570 
358@ 
35:30 
5602 
3612 
5620 
3650 
3640 
6) eye)’ 
5660 


| 
ISSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSS 
SUB Voltsms(Volt(*#) J ,Ums) 


REM FHEFEF HEHE EFT H ALTE HFEF ELH FFL EE EFAS FELT HE FFE PLE FHF H+ +E Hees eee tegteete+ 
REM THIS MODULE DETERMINES THE MEAN SQUARE VALUE OF THE INPUT VOD ieaes 
REM ARRAY. 

REM FEFHEFAF HEHE FA F HEHEHE FHL ELEAF EF HFFA FFF FAT EF HEF AF HE FFE +L HEHEHE HHH Het ttt 


eee MODULE VARIABLE DECLARATIONS AND DEFINITIONS +4++44++++++++++++++4+4++ 
! 


THE ARRAY REPRESENTING THE INPUT 


(REAL VOLT(*) 
rm; VOLTAGE ARRAY. 


REAL Volts THE SUM OF VOLT(4)°S ELEMENiGas 

REAL V THE DIFFERENCE BETWEEN VOR ies 
AND VAVG. 

PREG, Vins THE MEAN SQUARE VALUE OF VOL T= 


| 
| 
| 
| 
| 
REAL Vavg ! THE AVERAGE VALUE OF VOBIEae 
| 
| 
| 
| 
| 


PHAAEHEEEHHHHE HEHEHE HEE HEE HEHEHE FEE HEF H EEE HEHEHE EEE HEE HEHEHE EHH ETH Hee te tte 
l++++ DETERMINATION OF UMS ter tttt et tHe tHe tr HHE HHH HHH HHH HHH Het tte teste t 
| 
| 

Voltis=2@ 
FOR 1-0 102d 
Volts=Voltst+Volt( I) 
NEXT I 
Vavg=Volts/(Jt+l) 
FOR [=@ 10 J 
V=eVolt(I)-Vavg 
Ums=(V#V )+Ums 
NEXT I 
Ums=VUms/( J+1) 
| 
| 
PHRF HHHE HEHEHE FHL ELE HH FHS FHF HHH HE HEHE FEET FE FH FEL ET HEP EE EEF HE TH HTH Het ree ttt+ 
b++4++ MODULE CONCLUSION +44 tttttrttte ttt ttt etete eee te tte tte tt tt+ t+ treet 
SUBEND 
PHAEHHEHEHEETHEHHEE HE FE EF HEE EEE HELE HEHEHE FE HE EHH EHH tHe He te tt teeter tett rere tttt+ 
|SSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSISSSSSSSS 
| 
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3670 
3682 
369@ 
3702 
371@ 
3728 
E750 
3740 
B750 
376@ 
3770 
378@ 
27390 
3802 
381@ 
3820 
383@ 
3842 
3850 
386@ 
387@ 
3882 
3892 
3902 
5310 
3920 


CALL Atten(Tc ,Rh,Patm,C,At ,Efreq) | 
R=Tgttm#C/2+A 1 R DETERMINED. 
| 


SUBEND 


| 
ISSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSS 


SUB Range(Tgttm,A,Te ,Rh,Patm,C ,At ,Efreg,R) 


REM FALELLLE EHF EFEEEEFEF HEF ELE FE FLEE FLEET FHA FAFHE FF H EH E+ HH tt + ttttttttttt 
REM THIS MODULE DETERMINES THE AVERAGE PATH ATTENUATION AND THE RANGE 


REM BETWEEN THE ARRAY AND THE TARGET CENTER. 
REM FEFEEEFEEHRE EF E EE FE FE EEEE FE FEELALE AFA LHEAFF FH TH+ HHH HET + HH ttt ee ttttttset 


| 

l++++ DETERMINATION OF R AND AT t+t+ttttttttttttt He HH HH ttHtH Het tte tte tttettet 
i 

Tattm=Tgqttm/1@@@ 

AlemNO CUE TERMINED. 


| 
PHAEEEFEE LEE EEE EEE FEF E EEF EEE HEE HEHE AEE EEE EAHA EFA PEE FHF EFT HE HH HH +t tte ttt 
1++++ MODULE CONCLUSION +t+tt+ttttt ttt tet tt tH HEHE HH Ht He tte tet tt ttt ttettet 
| 


LAFEFFE FELLER EL FFE LTE FEEL LFF TELE FE EEE ELE EEF FF AF EEE T HET HE t Ete ete ete 
ISSSSSSESSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSS 
! 


lige 


5930 
3940 
3950 
3960 
3972 
3980 
3930 
4202 
4010 
4220 
4230 
4040 
4050 
4060 
4072 
4280 
4090 
4122 
4110 
4120 
4130 
4148 
4150 
4162 
4170 
4182 
4130 
420 
421 
4220 
4230 
4240 
4250 
4260 
4270 
4280 
4290 
4300 
4310 
4320 
4350 


rs 


ISSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSHSSSSSSSSSSSSSSSSSSSSSSSSSS 
SUB Atten(Tc,Rh,Patm Cmoist ,At ,Efreq) 


REM 
REM 
REM 
REM 


REM 


REAL 


Remi 


REAL 


REAL 


Rene 


Rene 


REAL 


Rie mie 


REAL 


REAL 


KREME 


REAL 


FEFEEHEEFHEFLA EE EFE FFE EEE EEE HE FEE FFE LAFF FEF E FETE EHH +E HEHE HHH tet ett 
THIS MOOULE CALCULATES THE ATTENUATION USING RELATIVE aurea 
IN PERCENT, ATMOSPHERIC PRESSURE IN MILLIBARS, AND TEMPenn ie 


IN DEGREES CENTIGRADE. 
FHEFFFHEHEHEFAFHE ETE FE FHT HEF EEF HEEE EEE EEE FAFEFE FE EEF EE H+ Hee te eee eettttt 


1++++ MODULE VARIABLE DECLARATIONS AND DEFINITIONS +++++++++++++++++++++ 


1 
Tk ' THE AMBIENT TEMPERATURE IN DEG. K. 
| 
Es 1 THE SATURATION VAPOR PRESSURE. 
l 
Pratio ! THE RATIO OF WATER PRESSURE IN 
| MILLIBARS TO ATMOSPHERIC PRESSURE 
| IN MILLIBARS. 
| 
Cary | THE SPEED OF SOUND IN ORY AIR. 
| 
H | 
| 
Past | P ASTERISK. 
| 
Tast ! T ASTERISK. 
| 
Fm | MAX. FREQ. 
| 
Amax | MAX EXPECTED ATTENUATION. 
| 
Frat | EFREQ/FM 
l 
F2 | FRAT#FRAT 
| 
Acl Amol i ATTENUATION EXPONENT COMPONENTS. 
| 


LHAFEEE FLEE EFL F HFEF F EF EE FHF EF HEEL EFA F EEE FETE FE FH ETE ETE FEAF EF E+E E +E HH HEH tte tt 
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4340 
4350 
4362 
4372 
4380 
4330 
4400 
4410 
4420 
4430 
4440 
4450 
4460 
4470 
4480 
4490 
4502 
4510 
4520 
4530 
4540 
4550 
4560 
4570 
458 
4590 
4620 
4612 
4620 
4630 
4642 
4650 
4662 
4672 
4682 
4692 
4700 
4712 


(++++ CALCULATING THE ATTENUATION ++t+t+++++t+t +t ttt ttt tte ttettetetetetett 
| 


CONVERSION OF TEMP IN CELSIUS 10 
TEMP IN KELVIN. 


re Ter2 (Sols 


eo=10°(9.4-2553/ 1k} 


| 
| 
l 

Rh=Rh/ 100 1 CONVERSION OF RH IN PCT. TO PRATIO. 
| 
Pratio=Es*#Rh/(Patm-Es+*#(1+Rh)) | 
| 
| 


Cdry=20.05*SQR( Tk) CmeculAltON OF THE SPEED OF SOUND. 


Cmoist=Cdry#(1+.14#Pratio) 


H=Pratio+!@0 
Past=Patm/1014 
Tast=(1.8#7Tc+492)/S19 
Fm=(10+6600#H+44400#H#H )*#Past/(Tast*.8) 
| 
Amax=.0@078*Fm*( Tast*(-2.5))*EXP(7.77#(1-I/Tast )) 
| 
Frat=Efreq/Fm 


F2=Frat*Frat 
| 


moci—Amax*oORC. 18%. 18*F2tcZz*F2/(1+FZ2))*2)/304.8 
Acl=1.74#10°(-10 )*EfreqtE freq 
At=(Amol+tAc! )#*L06(1@)/10 

| 
PH HEHEHHHFE HEHE EEE HEEEE FE HHE HEHEHE EEE EE FHF HFEF ELEAF HEHE HEHEHE FEF HEHEHE FH tte teetttet 
| 
b++++ MODULE CONCLUSION +++ +++ tt ttt ttt ttt ttt tte FHF tHHH te ttt tt tttttttttttett 
SUBEND 

| 
PHEAE EEL HEHEHE HHEHEE HEHEHE HEHE FE HEE FEE EHH HEHEHE HEE HEE FE FET H +H +++ ttt tttet+ett+ 
| 
ISSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSHSSSSSSSSFSSSSSSSSSSSSSSSSSSSSSSSSSSS 
| 


Ls 


4720 
4730 
4742 
4752 
4762 
477 
4782 
4790 
4802 
4810 
4820 
4832 
4842 
4850 
4860 
4870 
4880 
4892 
49200 
4912 
4920 
4930 
4940 
4952 
4962 
4970 
4980 
4990 
5000 
5Q1@ 
5020 
5032 
5040 
5052 
5060 
5870 
5282 
509 
5180 
511 
5122 
5130 
5140 
5152 
5160 
5172 
5182 
5190 
5222 


| 
ISS$SSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSS 
SUB Bcksctr(Ka ,Nbs) 

REM HHAFHA HEHE FHEA FAFA EEE EA FAF HFEF HEAL E FEE EEE FFE EAE HE FF HEHE F4E EHH eH 4e4eett 
REM THIS MODULE CALCULATES THE BACKSCATTERED CROSS SECTION NORMA Aas 
REM BY THE CROSS SECTIONAL AREA OF THE SPHERE: 
REM FHEHHHE HEHEHE HEHE He Het te tet e ete te tet eet tee te te ¢e¢ete te tee te tet+4e++4+++ttes+ 

1 


I++++ MODULE VARIABLE DECLARATIONS AND DEFINITIONS +4+4+++4+4++4+4+4++4+4+44+44+4+444 
l 

REAL L ORDER OF THE FUNCTION CALCULATED 

IN A SUBROUTINE. IT IS USED 

AS A LOOP INDEX AND TO DESIGNATE 

THE ARRAY ELEMENT FOR THE FUNCTION 

OF CORRESPONDING ORDER. 

REAL Lmax MAX. SIGNIFICANT ORDER OF THE SUM 

DETERMINING THE NORMALIZED 

BACKSCATTERING CROSS SECTION 

(NBCS). 

REAL J(S1) ARRAY FOR THE REGULAR SPHERICAL 

BESSEL FUNCTIONS (RSB) OF 

ARGUMENT KA AND ORDER @ THROUGH 

LMAX. 

REAL Bo.) ARRAY FOR THE DERIVATIVES OF THE 

RSS8S wepHE DRSB): 


BESSEL FUNCT IONS (isege age 
ARGUMENT KA AND ORDER @ THROUGH 
LMAX. 
REAL Dytat > ARRAY FOR THE DERIVATIVES OF Sige 
ISB “ThE DISBo. 
REA PiSl,) ARRAY FOR THE LEGENDRE 
POLYNOMIALS OF ARGUMENT COSC IS@ae 


NEW ec FACTOR COMMON TO RSM AND ISM. 


REAL Rsm THE REAL COMPONENT OF THE 


SCATTERING MODULUS eo. e 


THE IMAGINARY COMPONENT OF THE 
rte 


l 
| 
| 
| 
| 
| 
I 
| 
| 
{ 
| 
| 
| 
| 
| 
REAL Y(S1) ! ARRAY FOR THE IRREGULAR SPHERICAL 
| 
| 
i 
| 
I 
I 
i 
| 
1 
| 
| 
| 
| 
REAL Ism 
| 
| 
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5212 
5220 
5230 
5240 
5250 
5260 
5270 
5280 
5290 
5300 
5312 
5320 
5330 
5340 
SSS 
5360 
5372 
5380 
5390 
5400 
5410 
5420 
5430 
5440 
5450 
5460 
5470 
5480 
5490 
5500 
5510 
Bol 
SSRN 
5540 
S550 
SS6Q 
=a 
5580 
5590 
5600 
5610 
5620 


SGnle Ursbitmax ,Ka ,0j{ *),J(+*)) 


| ++++ NOTE ++++ 

| THE ARRAY DIMENSIONS (30+21) NEED 
IPO se INCREASED 10 ACCOMODATE 

| KA>3@. 

| 


! 
et ee ee ee tte ee ee eee eee tte ett eT tt tt te tte te ett ttt t rt 
I++++ CALCULATING THE NBCS tt ttttttt ttt tt tt ttt ttt ttt tt ttt tttttetttttt+tt+t 
! 
Lmax=INT(Kat21 ) 


CneCULATe. THEBOERIURTIVES OF THE 
REGULAR SPHERICAL BESSEL FUNCTIONS 
OF ORDER @ THROUGH LMAX. 

CHeECUERS) PHe DERIVATIVES OF THE 
TRREGUCAR SPRERICAL BESSEL FUNC NS 
OF ORDER @ THROUGH LMAX. 


CALL Disb(Lmax ,Ka ,Dy(*#) YC *#)) 


Pole eegiovt imax Ka ,Pt )) COLCURMies THE ARRAY FC): 


CALCULATING THE REAL AND IMAGINARY 
COMEONENT Se Or wile SM. 


Rsm=Q 

Ism=0 

FOR L=0 TO Lmax 
deci tiet ee C7 FL ey LE +Dy( EL )+Dyt(L)) 
Rsm=Rsmt+(Fac*Dy(L)) 
Ism=Ismt(Fac+#Dj(l)) 

NEXT £ 


— _—— + ae — _— + —- — + — — + 


{ 
Nbs=4+(Rsm#eRsmt+tIsmeIsm)/(KatKa) i CALCULATING. THE NBCS. 
! 


| 
PEF FEE HEE HE FEL EF HEF FEF FHF FALAFEL EFL E LTE LEE EEE EL EE HEEFT HEH E+ tet tet tet sete 
! 
b++4++ MODULE CONCLUSION +4 +444 44444444 4 +t FF FHF ttt ttt tt tt tetett+ettttt++++ 
! 
SUBEND 
| 
PEHFL HEF EEFE FFF ALE EE ELE LEELA LEH E LFA TEE THEE TEEPE FH AEE HEHE Hee ttt tsetteeetettt+ 
ISSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSHSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSES 
| 


an 


5630 
5640 
5652 
5660 
567@ 
568@ 
5630 
5700 
Sie 
572@ 
S7o0 
5740 
5750 
5760 
S770 
578@ 
Sv a0 
5822 
5810 
5820 
5830 
5840 
5850 
5860 
5870 
5880 
589@ 
5302 
591@ 
5920 
593@ 
594@ 


ISSSSSSSSSSSSSSSSSSSSSSSHSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSS 
SUS Orsb(Lmax {he 0 )\ *) 2)? 

REM FEET HE HEHE HEHEHE HEHEHE EEE HEE ELE ETE HE EAE EEL FH FEF EH HEH HH tet eeetetee+eteet 
REM THIS MODULE CALCULATES THE DERIVATIVES OF THE REGULAR SPHERICAL 
REM BESSEL FUNCTIONS (ORSB) OF ARGUMENT KA AND ORDERS @ THROUGH LMAX. 
REM THE CALCULATION OF- THE DERIVATIVE INVOLVES THE RSB”S OF SRE eae 


REM AND PREVIOUS ORDER. 
REM FHEHFFTETHF FHF HEHEHE EEF E FET ELAL LEFF FE FEEL FEF EEE LEE FETE TEE F+ ++ e+ +e eet ete 


! 
l++++ CALCULATING THE DERIVATIVES OF THE RSB ++t+++ttttttt+ttteettttttttet 


CALL Rsb(Lmax ,Ka,J(*)) CALCULATE THE RSS VALUES Jitgaam 
ANGLE UNITS FOR DRSB CALC hice 
CALCULATE THE INITIAL ORSS Die 
CALCULATE THE REMAINING ORGSae 


RAD 
D3(@)=(-J3(@)+COS( Ka) )/Ka 
FOR t=! TO Lmax 
Oe ee re en) Se lt 
NEXT L 


! 
! 
| 
! 
! 


| 
PAA E HEHEHE FEEL HEHEHE HEEL HEE HEHEHE HEE HEE HEE EEF EE HEF HE +t tte te ettet+ttttet+t+ ttt 
! 
b++++ MODULE CONCLUSION +4+t+t+t+t+ +++ t+ ttt ttt ttt ttt rt tt tttt te tttttettett+ttt+ 


SUBEND 

| 
PELE HEHEHE FEF HEHEHE FE FHEFL FLEE FEET E LTE F EEF E FE FEFE HFEF E+E+ HEH HEF +e tet ttt+t++et 
| 
ISSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSS 
{ 


lize 


S550 
S960 
S370 
S38 
5990 
6200 
6Q10 
6222 
6232 
6240 
6052 
6062 
6272 
6280 
6092 
6180 
6110 
6122 
6130 
6142 
61S5@ 
616 
6172 
618 
613 
6202 
621 
6220 
6232 
6240 
6250 
6262 
6270 
6282 
6292 
6320 
6310 
6320 
63350 
6340 
E352 
63560 
6370 
6382 


ISSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSHESS 
SUB Rsb(Lmax ,Ka,J(*#)) 

REM FEEEEE EEE HEE FEE EEE FE FE FHL EEE EEF EEE EEE EE FEE EEE FEE FEE Feet tee tet ee ttettet 
REM THIS MODULE CALCULATES THE REGULAR SPHERICAL BESSEL FUNCTIONS 

REM (RSB) OF ARGUMENT KA AND ORDERS @ THROUGH LMAX. Wale 

REM CALCULATION USES THE CONTINUED FRACTION APPROACH. 

REM soa seyhgpsp ga gnse cose Se ecb Sh age Sasa E aaa dea a asian orp chen ela sa ota ease ohaah sateen soon S1G Za papa smae Se one alec 
! 


[++++ MODULE VARIABLE DECLARATIONS AND DEFINITIONS +4+4+++4+++++++++4+++++++ 
| 


OF BNM. 


| 
REGAL RI ihe sey OU TNGmr On Oe. )/ J Clee 
| 
REAL Numer ISTHE EVOLVING NUMERATOR FACTOR USED 
Pet Ne CALCURATION OF RE: 
REAL Denr | THE EVOLVING DENOMINATOR FACTOR 
[SUSE TNeCAECUEAT ION SOF RE- 
l 
REAL Scexp PeEVOEVING SCAUING EXPONENT FOR VERY 
| LARGE OR VERY SMALL VALUES OF RL. 
| 
REM Nf lag | @ FLAG FOR NUMR. 
REAL Dflag ! @ FLAG FOR DENR. 
| 
REAL Bnm PeeVOEVINGSTERMSUSED INSTHE 
| CALCULATION OF NUMR AND DENR. 
REAL Binc ! INCREMENT USED IN THE EVOLUTION 
| 


| 
Peete e tt ete tte tet tt eet ett tt ttt ttt t+ teeter ttt ett tt ttt ttt tttettttt t+ t+tttttset 
| 
l++++ CALCULATING THE RSB ++t+tttttttttt ete ttt eet tttt tte te tte tttttetetetetet 


RAD ANGLE SUNT TS FOR CALC NS. 
J omenecueate THe INITIAL #+« 


se* RSB, J(Q). as 


| 
| 
| 
| 
J(Q)=(SIN(Ka))/Ka | 
| 


eZ 


639@ 
6482 
6412 
6420 
6430 
6442 
6450 
6460 
6470 
6482 
6490 
65280 
651@ 
6520 
6530 
6540 
6550 
6560 
6572 
6580 
6530 
6600 
6610 
6620 
6630 
6640 
665 
6662 
6670 
6680 
6690 
6700 
6712 
6720 
6730 
6742 
6750 
6762 
6772 
6782 
679 
6800 
6810 
6820 
6830 
6842 
6852 
686 
6872 
6880 
6890 
6900 
6912 
6920 


FOR L=!1 TO Lma~x 
R1=1.@ 
Numr=2. 
Denr=0. 
Scexp=0. 
Nflag=@. 

Of lag=@. 
Bnm=1.@/Ka 
Binc=2.0/Ka 


FOR Ordcent=! TOL 
Bnm=BnmtBinc 
IF Nflag=@ THEN 
Numc=Bnm-Numr 
IF Numr<>@ THEN 
R1=R1*Numr 
Numr=!.0/Numr 


TF ABSCR I) 10 S256 Tren Gack 
ELSE 
Nflag=i 
ENG Sie 
ESSE 
Nf lag=@ 
ENO: ste 
NEC Cedent 


oo —— me .—— :— om et ee —— —— — — — + 


I 


##* CALCULATING THE RATIO #+#« 
##e J(QO)/I(L). fae 
INITIALIZING Sai 

INITIALIZING NUMR. 
INITIALIZING VENRE 
INITIGEDZING See 
INITIALIZING NFLAG. 
INITIALIZING DEEAG: 
INITIALIZING BNM. 

CALCULATE BING 


* CALCULATE THE L “NAKED" * 
* NUMR FACTORS OF J(@)/J(L). * 
INCREMENT BNM. 

CHECK NUMR=@ FLAG NOT SET; 
UPDATE NUMR IF NUMR<O@. 


UPDATE RL IF UPDATED NUMR Ge 
PREPARE NUMR FOR NEXT EVOLUTION. 
CHECK SCALING OF UPOATEDSRES 


Scale(R1] ,Scexp) 


| 
| 
| 


IF NFLAG<>@ THE PREVIOUS STEPSeiiae 
DEFERREO FOR ONE INCREMENT OF SSRs 


* CALCULATE THE REMAINING NUMR * 

* AND DENR TERMS OF J(O)/J(L). * 
NOTE WHEN NUMR=DENR THEN 
RL=J(@)/ ICL). 


WHILE Numr<>Denr OR Nflag=! OR Of lag=!t 


Bnm=BnmtBinc 


IF Nf lag=@ THEN 
Numc=Bnm—-Numr 
IF Numr<>@ THEN 

R1=R1 *Numr 
Numr=1.0/Numr 


IF ABStR1 O10 256" THEN. Cages 
PESE 
Nflag=l 
END IF 
ELSE 
Nflag=@ 
END IF 


IF Oflag=@ THEN 
Denr=Bnm-Denr 
IF Denr<>@ THEN 

R1=R1/Denr 
Denr=1.@/ODenr 
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Py 
. 


INCREMENT: BNM. 


CHECK NUMR=@ FLAG NOT SET; 
UPDATE NUMR IF NUMR<O@. 


UPDATE RL IF UPDATED NUMB@ Ze 
PREPARE NUMR FOR NEXT EVOLUTION. 
CHECK SCALING OF UPDATED REF 


Scale(R! ,Scexp) 


IF NFLAG<>@ THE PREVIOUS STEP Sie 
DEFERRED FOR ONE INCREMENT OF BNM. 


CHECK DENR=@ FLAG NOT SET: 
UPDATE DENR IF DENR ZO Ge 


UPDATE RL IF UPDATED DEIR 
PREPARE DENR FOR NEXT EVOLUTION. 


6930 
6940 
6952 
6962 
6970 
698@ 
6990 
7002 
7010 
7020 
7030 
7040 
7050 
7060 
707 
7082 
7090 
7100 
7110 
7120 
€130 
714 
715@ 
7160 
7170 
718 
Tpely 


! CHECK SCALING OF UPDATED RL. 
IF ABS(R1)>10*%25@ THEN CALL Scale(R1 ,Scexp) 
Bor 
Dflag=!I 


END IF 
EUSE Peimveenoa70 ihe PReVvlOUS STEPS s ARE 


Df lag=0 | DEFERRED FOR ONE INCREMENT OF BNM. 
SNORE 


END WHILE 


«s+ J(L) DETERMINED. «++ 


| 
| 
mE I=100)/CRI*10° (Scexp)) | 
| 


NEXT L 
| 


| 
PEELE EAL HEHEHE EEA F LEA ALE AL HALAL ETHEL EEE HEF H+ +e FH ett tt t+etetettt te tettt+ 
| ; 
f+4+++ MODULE CONCLUSION +4+ttt+tttttttt ttt ttt te ttt tte tttetettttetet+tett+t+++ 
| 
SUBEND . 
| 
PEELE LEELA LEE EEE LA LEA LALLA LAL HEEL LA LEAF ALA FEE EFF ++ tt teet+etette+et+++t+ 
(S$SSSSS$SSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSES 


W2> 


7200 
721@ 
7220 
aay) 
724 
7250 
7260 
7270 
7280 
7290 
730@ 
T3510 
7320 
7530 
7340 
1250 
7360 
7370 
7380 
7530 
7400 
741@ 
7420 
7430 
7442 
7450 
7460 
7470 
7482 
7490 
75200 


ISSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSES 
SUB Scale(R1 ,Scexp) 

REM FETE EEEE EEE EEE FEET FEE FEELS E HELE FETE ELE HEEFT EHH TEE F tee eee etettet+tte+ 
REM THIS MODULE SCALES RL TO MAINTAIN NUMERICAL ACCURACY IN The 


REM ONGOING CALCULATION OF RL. 
REM FEF EHEFEFEFEH EEE EE EEE FEE LE FELE LL LEAF EEE LEFF EEF EEE E TEE HE ++ +++ 4444444 


I 
1++++ SCALING RL ttt tHHEtEHHE tHE tHE HEHEHE HEHE tHE HHH HE tet tte tttettt+ettt+tttt 
| 
| 
IF R1>10*250@ THEN 
R1=R]1+10*(-250) 
Scexp=Scexpt+250 
END IF 
TP R1IiG?O (2500 eh 
R1=R1+#10°250 
Scexp=Scexp-25@ 
ENO IF 
| 
| 
PERERA EEE EE EEE LE HEHEHE HEHEHE EEL LE LAE ELE HHL LEH ELLE HHH HEHEHE HEHE HH ttt tee tete ttt 
1++++ MODULE CONCLUSION Ft tte te tHE HEHEHE EEE EEE EE EEE EEE EEE EEE EHH Ht te t+ ett 
I 
SUBEND 
| 
PERE HEHEHE HEHEHE LEE EEE HEE EE HEHEHE EE HEE EE EEE HEEL EHH EEE ETE HEHEHE FETE HH tre tetete ett 
| 
'SSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSHSSSSSSSSSSSSHSSSSSSSSSSSSSSSSSSSSS 
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7518 
7520 
(ees, 
7540 
(SM) 
7562 
7570 
7580 
7 aa0 
7680 
7612 
7620 
763@ 
7640 
7650 
766@ 
7678 
768@ 
7698 
770@ 
7710 
7720 
iroU 
774@ 
7758 
7760 
7770 
7780 
7730 
7802 
7810 
782@ 


] 
[SSS$SFESSSSSESSESSSSSSSESSSSSSSSSSSSSSSSSSSSSSSSSSSSSESSSSSSSSSESSSSSSSESFSSSSSS 
SUB Disb« imax ,Ka ,Dyt*), YC *)) 

REM FHEEEEF HEHE THEFT THLE FEF TELEF FH TE EEE FETE EE HEHE EE HEHEHE H+ Het tee tt ttt tt+tt+ 
REM fe en OOUEEeCCeCURAeomiee DERIVATIVES OF THEsIRREGULAR SPHERICAL 
REM Besse FUNCTIONS OF ARGUMENT KA AND ORDERS OeiHROUGH LMAX. THE 
REM Peeve velo AOE winemten valve wSeS THE 158"S OF THE SAME AND 


REM PREVIOUS. ORDER. 
REM FEAL ELEE FEEL EEEHEE HEHE FLEE FEE HEHE SHEET ETE HH tte etetetetettt+ett+t+tte+ 


! 
eae CwEn Nomi ewern VIVES OF THE LoS Tete ttt++++t++tte+t++++44++4 


CALL Isb(Lmax ,Ka ,Y(#)) CAREUEATE IRE sioseuMlUES. 
ANGLE UNITS FOR CALC’NS. 

Chew IN Ii TeabS i Se DERIVATIVE. 
(RoemoriuNINGo tos DERIUAIIVES. 


RAD 
Dy(@)=(-Y(Q@)+SIN(Ka))/Ka 
FOR L=1 TO Lmax 
yO = et Kade (iy (1) 
NEXT L 


. 
. 


| 
PELEFHEHETEFEFHE HEEFT TELE FEET FEF ETE FEET LE FET TELE FETE EEE EE +H H+ tt tte tettetteetttst 


I 

f++++ MODULE CONCLUSION +++ ++ 4+ +++ ttt ttt tttt ttt tte te tr ttt tet tetttttet+tt+e+ 
SUBEND 

LEEEEHEEE EEE HEHEHE FEE FE FEE EEE FEF EEE EEF FEE HFEF HEFT HE FHF HEHEHE HEH Ht tte ett ttet 
ISSSSSSSSSHSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSS 


Aad 


7830 
7840 
7852 
7862 
787@ 
7882 
7890 
7982 
791@ 
7320 
793@ 
7940 
7950 
7960 
7327 
7382 
7990 
8282 
8010 
8822 
8830 
8240 
8250 
8850 
807 
8282 
8290 
8120 
8110 
8122 


ISSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSHSSSSSSSSSSSSSSSSSSSSSSSS 


SUB Isb(Lmax ,Ke ,Y(*)) 
REM FEFFEEFHEL FEF EE EEF HEHEHE FEEL EF EF ELE EEF EE FEAF EEE FEF FEE Ht ete tet tetetee+te+ 


end THIS MODULE CALCULATES THE IRREGULAR SPHERICAL GESSEU Uae ites 
REM OF ARGUMENT KA AND ORDER L THROUGH LMAX. THE CALCULATIONS 


REM ACCOMPLISHED BY FORWARD RECURS TOT] 
REM SHEE EHH LEE E FE FEE EEF HEE FEET E EF EEE FEE FE FEET FEE FFE + tt tet tettt+tetet++ttet+ 
t++4++ CALCULATING THE ISB +++++++ ttt ttt +t FFF tt tt ttt ttt tt tt+etet+tetetettt++tt+ 


| 
RAD 1 ANGLE UNITS FOR CALC'NS. 
Y(@)=-C0S(Ka)/Ka | INITIAL I1SB°S CALCULA TSE 
Y(1)=(Y(@)-SIN( Ka) )/Ka 


FOR L=2 TO Lmax | REMAINING 156° SeCare me 
YC J=¢C24 Soy Rape Cy Cie 2) 
NEXT L 


| 
| 
LEEEEHEFHEEFHE EEE EE FE FEE FE EEF EF HE EFL EEF THEE FFL FTE ETE THEE EEE HHH tH ttetettttt 
(++4++ MODULE CONCLUSION ++t+t44t+4t+t Ft tt FFF FF +t F4 tte ttt tetttte+tt+tt+tte+te+ 
| 
SUBEND 
LEE EEEHEFE FEEL FE FEF HT E LE FEL ELLE FEF LEE EEL F EEE FF + HHH + +tt ttt ttttt¢++trtte+ 
| 
ISSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSESSSSSSSSSSS 


Pee 


8130 
8140 
8152 
8162 
8172 
8182 
8192 
8222 
8212 
8220 
8230 
8242 
8250 
8262 
8272 
8280 
8298 
8382 
8310 
8520 
8350 
8340 
8350 
836@ 
8370 
8380 
8392 
8422 
8412 


l 
1\S$SSSSSSSSSSHSSSSSSSSHSSSSSSSHSSSHSSSSHSSSSSSSSSSSSSSSSSSHSSSS$SSHS$SS$SSSSSSSSSSS 
SUB Leg!8@(Lmax ,Ka ,P(*)) 
REM FEEEFHE HEHEHE HFEF HEHEHE ETE HE FE FE FEE EF E HEF FFE EEF EEF EH HH Het t ete teeeetetttetttt 
REM THIS MODULE CALCULATES THE LEGENDRE POLYNOMIALS OF ARGUMENT -! AND 
REM ORUERTO THROSGHVENAX . ween CULATION IS ACCOMPLISHED BY FORWARD 
REM RECURSION: 
REM FEFEFEEEA THEE E FE LEE EEE EE FE EEE LEFT ELH LEE EEE EF HHH EHH ttt et ete tt+tt+ttt 
| 


eae: Eaeevea Gr HeseeceNUKe FOEYNOMNITALS tttt+t+tt+ttt+t+t+ttt+t ere ret 
| 


X=-] ARGUMENT. OF THE EEG Polk? s. 
| 
P(Q)=1 Peale ease OL See CALC D- 
P(t )=X 
FOR L=Z [0 Lmax IR ohimiiN New te ost OL Sa Gale “OF 
pile) Crate AL) Xe PCL =) )—C1=( 1/L) +P (La 
NEXT L 


! 
PEALE L HAF F FEAL EHF L FALE EHTEL EL FALE EEE E+E ETE HEHEHE Ht Ht ttt tet ttttte tsetse ttt 
! 
l++++ MODULE CONCLUSION +++4+4+4++t44t4+ 4444+ tt ++ 4+ tt ¢4¢¢4+¢+ ++ t+ ++ +++ +¢+4t++4+4 
| 


SUBEND 
| 


PEE EEEEFHEFEFEEE FFE TEE FEF ESE FEES E EEE FEES FEF TEESE FH EE HH tet settee tetetttes 
IS$SSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSHSSSSSHS$SSSSSSHSSSSSSSSSSSSSSSSSS 
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8420 
8430 
8440 
8450 
8460 
8470 
8480 
8490 
8500 
8512 
8520 
8530 
8540 
8550 
8SE0 
8570 
8580 
8530 
8600 
8612 
8620 
8630 
8640 
8650 
8662 
8672 


| 

ISSSSSSSSSSSSSSSSSSSSSSSSSSSSESSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSS 

SUB Gain3sO(VUrms ,Tco Rh. R Patm, ir El regqwea) 

REM HEF EEE EE EEF HE EEE EEE EEE TEEPE EEE EEE EEE EEF HE HE HHH Ht ttt ttt ttt tt ttttttt++et 

REM THIS MODULE CONTAINS THE DATA USED FOR DETERMINATION OF Et Gong 

REM 3.0V INPUT FROM THE POWER SUPPE = 

REM HEHEHE FHEE HFEF HEHEHE EF HEE EEF EE EEE HEE EEE EEE EEE HEHEHE HEH ttt tt tetttetttt+tet 

l++++ LOADING DATA ttttttttte tt ttt ttttt ttt tt tt ttt tttttttttettettttt+ttttt 
| 

RESIORE 650 

READ Urms,Tc,Rh,R ,Patm,iIr ,Efreq,Za 


630: DATA .426 ,21.0,50.0,5.30,1021.5 52.951 E75 SoUC eae 


LEFEFEEFETHEALE TFET FEFT FLEE E TEFL FL FEF FFE FEF HEF EE HHH tt tet te tttetettetet+tttt+tete 
| 
l++++ MODULE CONCLUSION +++4tt+tt+t ttt Ft tt tt ttt ttt tte ttttttt+tt+t+ttetttet+ 
| 
SUBEND 
| 
LEEFHE EEE HEFH FE LETHE FFF T FTTH FATE FEET FEET EHE FEF E LEFF TEE ETE HEHE ttt t+ tetetett 
ISSSSSSSSSTSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSS 
| 


Lore 


8680 
86 90 
8700 
871 
8720 
8730 
8740 
8752 
8762 
8770 
8780 
8730 
8800 
881@ 
8820 
8850 
8840 
8850 
8860 
8870 
8880 
8890 
8900 
8910 
8320 
8930 


(SSESSSSSSSSESSSSSSSSSSSSSSSSSSSSSESSSSSHSESSSSSESSSSSSSSSSSSSESSSESSESSSESSSSESSS 
SUB GaimcostVems leCern.n Patm,IiIr,Efreq,Za) 

REM FEEFEEEEEEEFEEEE EEE HEE EHE FEE EE FEE HEF EEE FHT FEE HEHEHE HHH ttt tet tttttttte tt 
REM THES =MODUEE CONTAINS THE Data USED FOR DETERMINATION OF EtGo WITH 
REM Se50 I NEUISEA Orme POWER SUPPLY. 

REM FEFEEEEFEFEE EEE FEE EEE EEE EEE EEE FETE EEE HEHEHE tHe te tet tttttttteteetttttett 
l++++ LOADING DATA ttetete ttt terete tte tt tt ttt te te ttete tee etttetetetetetetetttttt+ 
PEolORE G55 

READ Vrms ,Tc,Rh,R ,Patm,Ir ,Efreq ,Za 


Boo: Wiese o.2|.0,00.0.5,50,1021.5,4.467E-5 ,5000 ,15.8 


| 
PEEEEEEE EEE EEE EEL EEE EE FEE EEE LE FEEEE EE ELE EEE EEE EEE EEE EEF E EE +EE HEHE tt tet+ttt 
| 
{++++ MODULE CONCLUSION Ft tt ttt tHE tHE HH HHH HEHEHE HEHEHE tte ttetetettetetetetett 
| 
SUBEND 
PEE EEE HEHEHE FEE EEE EEE HEHEHE HEHEHE HEHEHE HEHE EEE FEE HEHEHE HEE FFE Fete tet tte teeeetteetes 
'SSSSSSSSSSSSSSESSSSSSSSSSSSSSSSSSSSSSSSSSSESSSSSSSSSESSSSSSSSSESSSSSSSESSSES 


om 


8942 
8950 
896@ 
8978 
8982 
8992 
9200 
9012 
9022 
9050 
9042 
905 
9060 
9072 
9080 
9090 
9122 
9110 
S20 
clo 
9142 
915@ 
916 
917@ 
9182 
sia 
9200 
g20 
9220 
9250 
9240 
d2o0 
9262 
9272 


| 
ISS$SS$SSSSSSSSSSSSSSSSSHSSSSSSSSSSHSSHSSSSSSSSSHSSSHSSHSSSSSHSSSSSSSSSSSSHSSSSSSS 
SUB Spherel(Dt(*),J,Tc,Patm,Rh,Tgttm,Dia ,Efreq ,Nuspkr ,Diaspkr ,Gpre) 

REM TFHEFFEFEE HEH FF HEHEHE EEE HEE FFE EEE FEE FEE FEL F HEEFT EET H+ Fe teste eee etetetetet 
REM THIS MODULE CONTAINS THE DATA USED FOR DETERMINATION OF Sere ee 


REM TARGEI NOS? CHE as ene 
REM FHEFEEHFF HH FEFHEL FFL FHFEEFAL HEHEHE EFF F EFF FF ETH FEE +E FE HEHEHE Hee te ee ee +ee+ +e 


b++++ LOADING DATA ttttt tt ttt ttt ttt tHe te ttt te Ht tte te tet+ tte ttt+ttt+et++etetstet 
| 
J=39 
Re SORE 9S. 
FOR I=@ TO J 
REAL DEC? 
NEXT I 
READ Tc ,Patm,Rh ,Tottm ,Dia ,Efreg ,Nuspkr ,Diaspkr ,Gpre 


Sl: DATA $9.54 ,9.4,8.32 ,6.42,3.88 ,. 35 ,—-2.07 aa 


DATA -9.64 ,-9.45 ,-8.32 ,-6.41 ,-3.85 ,—.69 2. Ge 
DATA 3.6,9.58,8.22 6.26 53.67 ,.%3952. S49 ee 

DATA -9.64,-9.39 ,-8.2 ,-6:21 ,-3.61 (—.6S 22595 Sale e oe e 
DATA 20.2,1016.8 ,52.8,239.31,.2545 (5000 [oe Oia 


| 
LEE EEHFE ELE FFF F AFH FE FEF EFF HELE FFE FEE HFEF FEE F HELE HHH +H tt rete tte tee seseestet 
| 
l++++ MODULE CONCLUSION ++++4++4+t4 Ft tt t+ +4444 44+ttt te ttteeettt+etteteteet 
| 
SUBEND 
| 
LEFF HFEF HE HELE FEL EF FEF HEHE FE EE PEE ELLE FEE EEF EPEAT FEF EFE ETE + HE HH t tee tee tse teseee est 
| 
IS$SSSSSSSSSSSSSSSSSSSSCSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSS 


| 2 


9280 
3290 
9300 
a1 0 
9320 
9330 
3340 
9350 
936 
2570 
9382 
9390 
9400 
5341@ 
9420 
9430 
944@ 
9450 
9460 
5470 
9480 
949@ 
950@ 
9510 
952@ 
953@ 
954@ 
J500 
9560 
9570 
358@ 
gag 
9600 
3610 


! 
ISSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSS 
SUB Sphere2(Dt(#),J,Tc,Patm,Rh,Tgttm,Dia,Efreq ,Nuspkr ,Diaspkr ,Gpre) 
REM FEFEEFEFHE FE FE FEEL EEE FE FEF EEE EE FE ELE FEEL E FEE FEE HEHEHE FE HEE TE FH HHH H+ Het tt 
REM ithowmogvee CONTAINS (HE DATA USED FOR DETERMINATION OF Er WITH 
REM Ramee aNOe 2 EMPEY ED: 
REM FHEEEHEFHEFE FE FE FHL EEE FEE EEF FE FEE FE EEF E EE FE FEE HEHE FEE FE HEHEHE HEHE tHe te tet 
i : 
l++++ LOADING DATA tttttttt FHF HEF HEF HHT HHHHE FE HE HHH+ tHe ttt tee ttetetttet+ ttt 
! 
| 
J=39 
ReslORE SZ 
FUR 1=0 10 J 
Remon Ot (1 ) 
NEXT I 
READ Tc ,Patm,Rh,Tgttm,Dia ,Efreq ,Nuspkr ,Diaspkr ,Gpre 


See cligete eee O64 wo o.2 0.2.0), .500,-. 940 -2.55,-3.52 ,-4.36 


Dele ao ORO oO mone S|). 9O,=,.500 ,.950,2.34,3.51,4.34 

ieee o Ob 4) Solidi 95,.500,-1.00 -2.40,-3.57 ,-4.38 

Oengtee ome bo 4 |S =o, 20 ,-tod2,-2.480,1.02 ,2.59,5.54 4.34 

Crilgee Ure od, 449.6 50755,.1013 5000 ,19,.0762 ,11094 

| 

epee et ee ett tee ee et ee ete ee et tte ttt ett tt tte tet 
b++++ MODULE CONCLUSION +++ +++ +++ +++ +4 4+ +44 $+ ++ ttt ++ ttt tttett+tt++t+t+t 
: ; 
SUBEND 
PSEHEEHEEE EEE PFE EE FE FEF FEEL EE HEHEHE EE FEET FE LHEE EEE EEE EEE THEFT ETE EHH ++ +++ +t+t+ 
ISSSSSSSSSSSSSSSSSSESSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSESSSSSSSSSESSESSHSSSHS 
| 


3622 
96 30 
9642 
ey SA, 
966@ 
carU 
968@ 
96 93@ 
37 G2 
9712 
9720 
973@ 
9742 
9750 
9760 
9772 
9780 
3732 
9800 
9810 
9820 
9832 
984@ 
3850 
9862 
987@ 
9880 
9890 
9982 
3910 
J220 
9930 
$94@ 
3350 


| 
ISSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSS$SSS$S 
SUB Sphere3(ODt(+) J ,Tc ,Patm,Rh,Tottm,Dia ,Efreq ,Nuspkr ,Diaspkr ,Gpre) 
REM LEFF FE FLA FEE FEAF H EEE HE FEL E FEL EEF EF EE HEE ETE PE HEH tee ee tee teteetee+ete 
Ren THIS MOOULE CONTAINS THE OATA USED FOR DETERMINATION OF eee 
Ren TARGET NG? 3 “ENPEGy aa: 
REM FEEELEHE AEE HEHE FEE F ELE LETHE EEF HEHEHE HH ttt ttt t+tteett+tt+ttetet+et++e+essy 
| 
b++++ LOADING DATA Fre tttHEHHHHETE HHH HHH Ht ttt ttt tte tetttttttt+ttt+ettt+ 
| 

J=39 
RESTORE Ss 
FOR I=@ TO J 

REAO OtCl) 
NEXT J 
READ Te ,Patm,Rh,Tottm, Oia Efrea Nusske Otesckr come 


S3: DATA 3.651,3.48 ,2.85,2.04 190) (=3120 —1 25 (2 


DATA -3.68 ,-3.46 ,-2.391] ,-2.07 ,-1.04 2110, 1 2562 ee 
DATA 3.63 ,3.39,2.84 2.02 ,. 370 =. 180 (-1., 30 2 ee ee 
DATA -3.69,-3.46 ,-2.89 ,—-2.04 ,—. 330 ,- 150.1 22252-26722 05 
DATA 21.3,10098.6,49.8 30.59, .07534 S000 >| 93 20 (a7 
| 
| 
PEEEEE FE FHE FFE FE FEE HELE FEE EF ELLE LEE ELF E FEEL FEL E EET HEP Et He tt tet +ete settee 
| 
b++++ MODULE CONCLUSION +++tt tt ttt ttt FEF HHH FE HHH HHT HH Ft EF tte eH tttettttete+ 
| 
SUBEND 
LHEFEHEEFHE ELE ETE HLF E FETE EFF E HEF E EET FEE EEE EE EEE EF HET H HEHE t e+e ett+e tet e+tettt 
ISSSSSSSSSSSSSSSSSSSSSSSSSHSSSSSSSHSSSSSSSSSH$SSSSSHSSSSSSSSSSSSSSSSSSSSSSS$SS 
I 
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9960 

9972 

9980 

S950 

1@@20 
10010 
10020 
10030 
1004@ 
10050 
10050 
10072 
1008@ 
18030 
101@@ 
1011@ 
10120 
101352 
10140 
1015@ 
10160 
101708 
10180 
10130 
10200 
10210 
10220 
102350 
10240 
10250 
10260 
10270 
1@280 
10230 


ISSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSHSSHSSSSSSSSSSSSSSSSSSSHSSSSSSSSSSSSSHSSS 
SUB Sphere4(Dt(#),J,Tc ,Patm,Rh,Tottm,Dia,Efreq ,Nuspkr ,Diaspkr ,Gpre) 


REM FEFEFEEFEE FEF E EE FE FLEE FEF EFFEF FEE EFF FFE FEF EEE FFE L EF E+E HE HHH HHH + +44 
REM Tals MOBVEE CONTAINS THe DATA USED FOR DETERMINATION OF Er WITH 
REM TARGET N@. 4 EMPLOYED. 

REM FHFHEEEEFE FFE F FEE FEEL FEL HEFEE FFF F FEL EE ELF HEFL FEFELEFL FEE E FFE EH +++ + 44444 


{ 

1++4++ LOADING DATA t+tttEtEHH FTTH FFEFFEFE HEHEHE FH FFFE FFE FE FFE FETE FH +t +t te te tt+ttete+ 
| 
| 

J=39 

RESTORE $4 

FOR I=@ TO J 

ReaD Dt( i) 
NEXT I 
Rene CeratmM nt Igttmebta freq Nuspkr Diaspkr ,Gore 


oe OG 2 2Oe lO) leer Ome eOe=. 520221) .03 ,-1.64 .-2.07,-2.351 
Bi Vimeo Ole 2 eee ee ee. 440 20 1.02 e622 2.06 ,2.28 
Cone Once Os, 400" —, 540 751.04 — 1766 ,-2.10,-2.32 
Oita ees oee wie Oat. 400 1.340 ,1.04,1.64.2.08 2.28 


Omnia cie4, (00976045 6750054..05240 .5000.19,.0762 ,11@34 
| 


i 

PEEREEEEFFEE EFF F ELE EFF E LL EFE FFF FE FFE EE FHF F FE FE FFE FE HE FFE HF tt ett tt tee te tststet 
| 

l++++ MODULE CONCLUSION +++++++++ +44 +44 ttt ttt ttt t ete tettettett teste t++tt 
SUBEND 

PREEEEEE FEE FFE HF EEE EFL FH LHL FE FHF FFF E LFF FFAS FEAF EFF FF FFF E HEHE +H HFEF + tet ttt 
ISSSSSSSSSSSSSSSSSSSSSSSSHSSSSSSSSSSSSSHSSSSSSSSSSSSSSSSSSSSSSSSSSHSESSSSSS$SS 
| 


San 


10300 
1031 
1032 
10332 
10342 
1035 
1@360 
1937 
1@382 
1039@ 
104@@ 
1041@ 
1042@ 
10430 
10449 
12452 
10462 
10472 
10482 
104930 
19582 
10510 
1052@ 
10530 
10540 
1Q@55@ 
1Q@560 
12572 
1Q@580 
109592 
1Q@602 
1061 


| 
ISSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSESS 
SUB Trans30( J) Dias) 


REM FH EEE HEEFT HEE EE EEE HEHEHE EEE HEHEHE HEHEHE HEHE HHH Ht tte tt ttt tte tttettte+ et 
REM THIS MODULE CONTAINS DATA USED IN DETERMINATIONS CRS sie 

REM 3.0) ABP Ten, 

REM FAAEFE EEE AEE L AEE HEF E FALE AFF FEF EEE FEAL E H+ HEH H+ tt tet ttttettt+et+ttset 


| 
1++++ LOADING DATA +++ +++ttt+tt+ttHtFHHHFHH FFF FHF FHF ttt tt ttt tte tteetet+t+tetttes 
l 

J=39 

RESTORE 130 

POR —Om1o. 

READ Dt(I) 

NE) 

T3Q@: DATA .690,.680@,.590,.440,.320,.160,-.920 ,-. 230) 341 7 
DATA -.610,-.610,-.520 ,—.360,—.240 ,-.080 ,.110 (3205) 4a 
DATA .690,.680,.590,.430,.320,.150,-.040 ,-.240 |— 4707 ae 
DATA -.62@ ,-.600,—-.510,-.380 ,-.240 .—.070 ,. 1203) 3330 3527 


1 
EEEEEEEEE EEE EFE HEHE FEF EFA F EE EE FEAT HEE AEE THEE HEHEHE FE HH tt tt tet ttettttetetett 
1++++ MODULE CONCLUSION 4+4+4+4+4444+44+4444+4444 44444444444 744444777 47 
SUBEND 
PHASE EAA EAE HEHEHE HEHEHE AFH HEHEHE FE EEE FEFE HEHEHE EHH FETE FEF HEHEHE Ht ttt tt++etttt 
| 
ISSSSSSSSSSSESSSSSSSSSSSSSSSSSSSSSSSSSSESSSSSSSSSSSSSSSESSSSSSSSSSSSSSESSSSE 
| 
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18628 
18638 
10648 
1865@ 
18668 
18678 
18682 
10698 
1878 
1071@ 
18728 
18738 
18742 
12752 
18762 
10778 
18780 
1879@ 
18882 
10810 
1882@ 
18838 
1884@ 
18852 
1886@ 
18878 
1@88@ 
18898 
189@@ 
1831@ 
18928 
1095@ 


| 
ISS$S$SSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSFSSSSSSSSHSSFSS 


SUB Trans35( J ,0t (am) 


REM FEFEHEFEHEEE HFEF EEE HEHE FTE FFF HFEF HE FEET FEET EFF F THF F E+ HH Hee tet e tee tty 

REM Tis NOGUEE CONTAINS lala USeweIN DETERMINATION OF Er WITH 

REM Boy free eo. 

REM FEFETHEFE HEHEHE FE FEE LETHE EFF FEE F TET HE FE FEET EHH FE +++ HH Ft ttt eee teee+stest 

| 

P++++ LOADING DATA +t+ttt+ttrt++ttHHttt HHH H+t tt +++ +e ttt tt+t+ttt++ttt+et+stteet 
{ 

J—39 

RESTORE 135 


FOR I=@ TO J 


READ Dt(I) 
NEXT I 
iMorerOninmee so .57710..640..496 5.290 .080 —. 180 —-.400.-.610.-—. 730 
Vim? 707. 720) — 7 5004-,450 .—, 240 —,020,.240..460. .660,.780 
WRiuees SO v7 7 Oe O00 470 1.290 ..060.—.190.—.428 —-.60e .—-. 720 
Clim CeO eae. a 705-470 —2250.—-.010..246 480. .679 .79¢ 
| 
{ 
LEE EEE EEE HEHEHE FEE EEE EE HEHE EEE EEE EEE EEE EE HEE HEHEHE EET HEHE HHH H+ tte tte teetttet 
| 
l++++ MODULE CONCLUSION t+t+t+4+ttt+ttttttttttttt ttt ttttttttetttetteeeteettite 
{ 
SUBEND 
| 
PHEEEEEEEEEHEHHEE EEE EHF LE FHEE FEE EE EEE FEEL EH EFE FEE HEHE EH HEH He t+ tee teeetteetest 
| 
ISSSSSSSSSSSSSSSSSSSS$SSSSSSSSSSSSSSSSSSSSSHSSSHSSHSSSSSHSSSSSS$SSS$SSSSSSSSSss 


15:7 


10940 
1@950 
10960 
103970 
1098@ 
10990 
1188@ 
11010 
1182@ 
11832 
1104 
11052 
11@6@ 
1187@ 
11080 
110392 
1110 
11112 
11120 
11130 
11142 
1115 
1116@ 
11172 
1118 
11192 
112QQ 
Phz 10 
11220 
1123 
1124@ 
eZS0 
11250 
11270 
11280 
1129 
11300 
1131@ 
Pio28 
Hie ove 
11340 
11550 
1136 
11376 
FISee 


| 


ISS$SSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSESSSSSSSSSSSSSSSSSSSSSHSSSSSSSSSSSSSSSS 


SUB Noise3Q(J,Jt ,Dt(*) ,Tdt‘*), Jen, ,Rhn.Patmn cimeone 
REM FEET EEEEFEEHEE FEE EEE EEE SESE FEE EEE FETE FETE FET EE HEE EHH HEH Hee etteeee tet 


REM THIS MODULE CONTAINS DATA USED IN DETERMINATIONS OSE aiian 


REM 5 CUM AREE LED. 


REM LHEFLHHHE EEL TEE E EEE FLEE EEE FELT EEE EEE EEE H HEHEHE FEE EF HEHEHE t+ ete ete tet 


l+t+4++ LOADING DATA FFttttt ttt ttt tt HH tH Ft tHE He ttt ete tte te tte ttt tt+etete+ 


J=S9 

RESTORE N3Q 

FOR f=0°30 J 
READ Dt 7 ) 

NEXT I 

Jt=39 

RESTORE Tn3@ 

FOR T=05 708s 
READ Tdt(I) 

NEXT I 

READ Ten,Rhn ,Patmn,Efreqn 

N3Q: 

alee eee | | 


l 
J 


DATA 0 ,@,-.02,4.05 ,-.04,-.@6 ,-. 06 231 ere 
Lge 08 .~.07 (US 5. 04 ae 01 Ome 


DATA @,9,-.02 42.03 ,-.05.-.07,-: 05) 4.1 a eee 


“HAEA 212 s=e lee te 
DATA .01,0,0,-.03,-.04 ,-.97,-.09,-.11 
-.08 ,-.07,-.04 ,-.03,0,0,.01 


BATCH: sce 


=,08 ,-.07 ,-..@4=. Osea eee 


foe ee 


Tn3@: DATA .59...68 595.44) 23200 O02 eo ee eS 


DEA e<6i ,22 5 bo. S22 Gee cere Crean 
DATA .53,.59 ,Sae4ie ae 


pen. 4 OS 


ISqom@A wens, - 42a Sa 


DATA >. 6223.57.51 ,—-2 90 5.24 Ol ee oe enc 


DATA 20.351. /,101b.7, 0000 


PERE EEL EAL FE HEF H EAE EHEEE ELE HHH LE THEE SHEET EEE THEE HEE E HEHEHE FHP HHP HH Het + +ettt+ 


1+4++ MODULE CONCLUSION +++ 444 Ft ttt tt ttt HHH HHH FFF HHH ett +e ett eet tetett+e+ 


SUBEND 


PEALE ALHE EE EF HEHEHE HH HEHEHE HEHEHE LEH HHH HELA EFHEF HFT HEHE THE HH HEH +++ t++ 


ISSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSHSSSSSSSSSSSSSSSSSSSSSSHSSS 
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b1390 
114@@ 
11410 
11420 
1143@ 
11440 
11450 
11460 
11472 
1148@ 
11490 
11502 
tat O 
11520 
1155 
1154@ 
11552 
1156@ 
nS 70 
1158@ 
lo JO 
116@@ 
1161@ 
1162@ 
11632 
1164@ 
1165@ 
1166@ 
1167 
1168@ 
11639@ 
11700 
tae 10 
ier 2 OQ 
11730 
1174@ 
| ash 
11760 
11770 
1178@ 
ri730 
11802 
1181@ 


ISSSSSESSSESESSSSSSSSSSSSSSFSSSSSSSSSSSSSSSESSESSSSSSSSESSSSESSSSSSSSSESHFSSSSSS 


SUESNewseso(Jert Otter) Tdty* olen Rhmaeraton ,Efreqn) 
FEET E EE HE FE HFEF EEE THE EET HEE EEE FEET EEE EEE FEE HH + tet ete tgtetet+eetete+ 


REM 
REM 
REM 


REM 


THis 
Sea 


HODUeE COMPAINS DATA USED IN DETERMINATION OF Er WITH 
APEIEIER. 


FHEFHEE HEHEHE THE THFE TE HEAFEAFEFEFEE FETE E FEET ATEFE FETE HHH ett tt tte ett 


l++++ LOADING DATA tt+ttttttt Ft tT Ht ttHtFttEHt tt tee tet tettttetet+ttetettt 


| 
| 


U-59 
RESTORE N35 
FOR I=@ TO J 
READ Ot I> 
NEXT I 
Jt=39 
REowere InSo 
BOR 1-21 Out 
Renee att] >) 
NEXT I 
READ Tecn,Rhn ,Patmn ,Efrean 
VecteemiOGtie = 0Ge— 4g OO ue 0? 7065 — 604 22 0, . OZ 
Doi ooenOomncor.0o..04..05,.02 .0,.-.02 —.05 
Pe Oe eens - 0G —.05 ,-.05 -.04 —.02 ,.01,.03 
Palins .Oeeedor.O/ 05. .04..01 =. 0le=.04 .-.06 
pose) he (4. Cone aan ec Oo seocmeero.-~.06.>.51 .—.53 
DP Neoware Pee So. 764. —-.46,.—-.5,-<l eee 4...6 
Dime Sree eeeolbeend .. Gol oo 16.08 —-. 35 7—.54 
Ciena 145-700 -.45 —-,75 —. 06. 165,.41,.62 
OAthee leo oso tOlb.6b ,SU0e 


THEA FLHE EHH FEF EF E FETE FHF ELELE FEF EFF FEET LF EAL ETH +E TH+ Ht tet tee ett eteeetteet 


l++++ MODULE CONCLUSION +4+tt+tttttttttetttttettttttttttetetttettt+etetettst 


SUBEND 


| 


| 


PEAEHE FEET HEATH EET E FATE E FTTH HEHEHE HTH FEF ttt t+ teste te tettetttteetet+ttetttett 


[SSSSSSSSSSSSSSSSSSSSSSESSSHSSSSSESFSSESSSSSSSSESSSSSSSSESSSSSSSSHSSSSFESESSESS 
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APPENDIX F 


ACOUSTIC ARRAY CALIBRATION DATA 


This appendix contains data used to determine the calibration 
product E fs E.. Data units are contained in square brackets. 

For listings of the same type of data the units and uncertainty are 
explicitly given for the first datum. Remaining data in the listing have the 


Same units and uncertainty as the first datum. 


|. TWO WAY PROPAGATION PATH DATA 
Pre-amplifier gain G, and array impedance / were previously 
determined by Moxcey. 6,was determined to be 11,094 Z was determined 


tobe 15.8 [Q]. [Ref. 8] 
The remaining data are presented in the following tables. 
V is the voltage supplied by the HP3314A function generator to 


supply 


the pre-amplifier (see Figure 8). 
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Table 5a: Two Way Propagation Path Data for Sphere One, 
Collected 30 November 1987: Run *2: Data Re-measured 
Because First Run V_ (with 3.5 [V] Supply) Was Clipped. 


Diameter = 2a = .2546 + 0.0006 [m] 
20 2 O71 [eee] 


Rh =52.8 + 0.1 [%] 
P = 1016.8 + 0.2 [mb] 
No. pulses = 20 


Veuppty = 3-0 [V] 


t. = 30.31 [mS] - 400 [uS] = 29.91 + 0.1 [mS] 


t V, t Vi, 

2.25+0.0I[mS] 690+ 5{[mvV] 25 690 
TRIAS: 680 2.46 680 
APA 390 ne) 390 
2.28 440 2.48 430 
2.29 520 2.49 520 
250 160 250 150 
Fae oy = 10 2.0 | - 40 
257 5) EZ =-2240 
250 =A a5 -4270 
Syl =5 50) 254 -560 
O55 ay 0 Zea) -620 
256 -610 Tes -600 
Zoo S20 yaa =5 0 
2.358 =oe 2.58 =—o50 
2.39 -240 2.59 =240 
2.40 - BO 2.00 =7/ 0 
23) mae 2.61 120 
2.42 520 Zo7 550 
yea 490 Zo 300 
2.44 630 Paley 640 
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Table 5b: Two Way Propagation Path Data for sphere One Continued, 


Collected 30 November 1987; Run *2: Data Re-measured 


Because First Run V_ (with 3.5 [V] supply) Was Clipped. 


t 


533.68 +0.01[mS] 


53.69 
J5270 
oo 1 
SS7 2 
SSS 
55.74 
Bye ic' 
2.76 
S77, 
35.78 
535.79 
53.80 
35:0) 
50 O27 
53:85 
55.84 
BK Rele 
B)Siele, 
53°67 


V 
r 


954+ 0.03 [V] 


9.40 
Gee 
6.42 
5.88 
0.960 
SZ.07 
-4.90 
Se 
Bo 07 
-9.64 
~9.45 
SG.o2 
-6.41 
535 


-0.890 


2.14 
4.94 
20 
8.88 
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33.88 
33.89 
33.90 
33.91 
55.92 
59,93 
33.94 
SoS) 
33.96 
55:97 
33.98 
55.99 
34.00 
34.01 
534.02 
34.03 
34.04 
34.05 
34.06 
34.07 


9.60 
eee, 
G22 
6.26 
eo 
0.730 
-2.30 
ola 
~7,.40 
~8.96 
-9.64 
TELS 
=o.20 
Oz 
=J,0)I 
-0.650 
256 
Se 
133 
8) 2) | 


Table 6a: Two Way Propagation Path Data for Sphere Two, 
Collected 28 November 1987. 


Diameter = 2a = 10.13 + 0.01 [cm] 
ier 2eO! 2 Ont [Xe 


Rh = 49.8 + 0.1 [%] 
P = 1009.4 + 0.2 [mb] 
No. pulses = 20 


Veupply = 2-3 LVI 


t. = 30.95 [mS] - 400 [y5] = 30.55 + 0.1 [m5] 


t Ve t Ve 

Zoo O1lmSleesso + 5 mV} 2.46 8350 
fa Juli 770 2.4/7 70 
L206 640 2.48 630 
Fa FOE) 490 2.49 470 
2.50 290 FAD 19, 290 
7a8).| 80 ae) 60 
Pag BF 1 3h0 Vasa -190 
235 -400 ZS -420 
2.54 19) 18, 2.94 -600 
IS, a7 30) ante) 128 
250 -770 2.06 -770 
237 -720 Zoe =700 
Zoo =550 250 =o 
Yas, -430 2.99 -420 
2.40 -240 2.60 Low 
2.4 ae 78) Zo! =O) 
2.42 240 2.62 240 
2.43 460 Z.09 480 
2.44 660 2.64 670 
2.45 780 2.65 790 
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Table 6b: Two Way Propagation Path Data for Sphere Two Continued, 
Collected 28 November 1987. 


ti 


V 


r 
34.13+0.01[mS] 472+ 0.01 [V] 


34.14 
34.15 
34.16 
34.17 
34.18 
34.19 
34.20 
34.21 

34.22 
34.23 
34.24 
54.25 
34.26 
34.27 
34.28 
34.29 
34.30 
34.31 

34.32 


4.66 
4.16 
O25 
2.0] 


0.560 
-0.940 


ge Ore 
os 2 
-4.36 
-477 
-470 
-4.18 
25225) 
-1.99 


=(300 
0.960 


2.34 
oS) 
434 
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34.33 
34.34 
34.35 
34.36 
34.37 
34.38 
34.39 
34.40 
34.4] 
34.42 
34.43 
34.44 
34.45 
34.46 
34.47 
34.48 
34.49 
34.50 
34.51 
34.52 


4.74 
4.66 
4.13 
3.19 
1.95 
0.500 
= ioe 
-2.40 
gee) 2) 
-4 38 
-4.78 
-4.68 
-4.15 
=320 
-1.92 
-0.480 
1.02 
2.39 
3.54 
4.34 


Table 7: Two Way Propagation Path Data for Sphere Three, 
Collected 28 November 1987. 


Diameter = 2a = 7.634 + 0.003 [cm] 
T. = 21.340.) [°c] 


Rh = 49.8 + 0.1 [%] 
P = 1009.6 + 0.2 [mb] 
No. pulses = 20 


Veupaty 7 3-9 (VI 


ee 0199 tins) — 400nS)!= 3092 Ons) 


V, data is the same as sphere two's (see Table 6a). 


t Vv. t Ve 
33.98 +0.01[mS] 3.61 + 0.03 [VI 34.18 3.63 
33.99 3.40 34.19 3.39 
34.00 2.86 34.20 2.84 
34.01 2.04 34.21 2.02 
34.02 1.01 34.22 0.970 
34.03 -0.120 34.23 -0.180 
3404 -1.25 34.24 -1.30 
34.05 -2.25 34.25 -2.31 
34.06 -3.04 34.26 -3.08 
34.07 -3.54 34.27 -3.56 
34.08 -3.68 34.28 -3.69 
34.09 -3.46 34.29 -3.46 
34.10 ~2.91 34.30 -2.89 
34.11 ~2.07 3431 -2.04 
3412 ~1.04 34.32 -0.990 
34.13 0.110 3433 0.160 
34.14 1.25 34.34 1.29 
34.15 2.24 3435 2.26 
34.16 3.00 34.36 3.05 
34.17 3.48 3437 3.50 
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Table 8: Two Way Propagation Path Data for Sphere Four, 
Collected 28 November 1987. 


Diameter = 2a = 6.240 + 0.008 [cm] 
T.= 21.4201 [°c] 


Rh = 49.6 + 0.1 [%] 
P = 1009.9 + 0.2 [mb] 
No. pulses = 20 


Veuppty = 22 (VI) 


t. = 31.04 [mS] - 400 [5] = 30.64 + 0.1 [mS] 


V, data is the same as sphere two's (see Table 6a). 


t V_ t Vv, 
33.83 +0.01[mS] 2.28 + 0.03 [V] 34.03 2.29 
33.84 2.10 34.04 2.10 
33.85 1.68 34.05 1.68 
33.86 1.10 34,06 1.09 
33.87 0.420 34.07 0.400 
33.88 -0.320 34.08 -0.340 
33.89 -1.03 34,09 -1.04 
33,90 -1.64 34.10 -1.66 
33.91 -2.07 34.11 -2.10 
33.92 -2.31 3412 -2.32 
33.93 -2.34 34.13 -2.35 
33.94 -2.12 34.14 -2.12 
33.95 “1.71 34.15 -1.70 
33.96 -1.12 34.16 “1.11 
33.97 -0.440 34.17 -0.400 
33.98 0.320 34.18 0.340 
33.99 1.02 34.19 1.04 
34.00 1.62 34.20 1.64 
34.01 2.06 34.21 2.08 
34.02 2.28 34.22 2.28 
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Table 9: Noise Data for 3.0 [V] Supply Voltage; 
Collected 30 November 1987. 


T.= 20.820.1{(%c] Rh=51.720.1{%] P= 1016.7+02{mb] 


Veuppty = 3.0 [VI 


Ve data is the same as sphere one's (see Table 5a). 


No. pulses = 20 


t Vo t V 

33.73 20.01{mS] 0025 [mV] 34.03 -120 
33.74 00 34.04 -120 
S345 = 26 34.05 -110 
33.76 - 30 34.06 - 80 
33.77 - 40 34.07 =e 
33.78 - 60 34.08 - 40 
33.79 - 80 34.09 = HC 
33.80 -100 34.10 - 20 
33.81 =) ike 34.11 00 
33.82 -120 34,12 00 
33.83 -120 34.13 10 
33.84 -110 34.14 00 
33.85 -100 34.15 00 
33.86 - 80 34.16 - 30 
33.87 - 70 34.17 - 40 
33.88 st 34.18 - 70 
33.89 - 40 34.19 - 90 
33.90 - 10 34.20 -110 
33.91 00 34.21 -120 
33.92 00 34,22 -120 
33.93 00 34.23 -120 
33.94 00 34.24 -120 
33.95 ae 34.25 -110 
33.96 - HE 34.26 - 80 
33.97 - 50 34.27 - 70 
33.98 = 70 34.28 - 40 
33.99 = 150 34.29 = 36 
34.00 -100 34.30 00 
34.01 -110 3431 00 
34.02 -120 3432 10 
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Table 10a: Noise Data for 3.5 [V] Supply Voltage; 
Collected 50 November 1987. 


T= 21.3401 [%c] 


Ri = SS = Onl tel 
P = 1016.6 + 0.2 [mb] 


No. pulses = 20 

Viton =3.5[V] 

t Mh t Ve 
2.24 +0.01[{mS] 740 +5 [mV] 2.44 750 
2.29 850 2.45 830 
2.20 820 2.46 810 
Ze 720 2.4/7 700 
2.268 530 2.48 500 
2.29 380 2.49 360 
250 180 2.90 160 
225) = 160) Z| = G0 
2.52 =5/0 fl aye =Joe 
755 =350 255 Sa l0) 
254 -680 Faye -690 
Ls -770 Yieaye -760 
26 -750 2.96 =7a0 
Zot -640 3 =6a0 
Ze -460 PA ay 6 - 460 
2.59 =500 2.99 -280 
ZAG -100 2.60 - 80 
2.4] 140 Ze | 160 
Zaz 400 2.62 410 
2.43 600 Zao 620 
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Table !0b: Noise Data for 3.5 [V] Supply Voltage Continued; 
Collected 30 November 1987. 


t Ve t Me 

33.80 +0.01[mS] - 60+ 5 [mV] 34.00 - 70 
33.81 - 80 34.01 - 80 
33.82 - 80 34.02 - 90 
33.83 - 80 34.03 - 80 
33.84 - 70 34.04 - 80 
33.85 - 60 34.05 - 60 
33.86 - 40 34.06 - A0 
33.87 - 20 34.07 - 20 
33.88 00 34.08 10 
33.89 20 34.09 30 
33.90 50 34.10 40 
33.91 60 34.11 70 
33.92 60 34.12 60 
33.93 60 34.13 70 
33.94 40 34.14 50 
33.95 30 34.15 AO 
33.96 20 34.16 10 
33.97 00 34.17 = 10 
33.98 - 30 34.18 - 40 
33.99 - 50 3419 - 60 
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2. ONE WAY PROPAGATION PATH DATA 
One way propagation path data are presented in the following tables. 


T_ and Rh are estimated. The estimates were based on hygrothermographic 


comparisons and 28 and 30 November data. 


Table 11: One Way Propagation Path Data for 3.0 [V], 
5000 [Hz] Continuous Wave Supply from the HP3314A 
Function Generator; Collected 7 December 1987. 


(V__), = 0.246 + 0.001 [VI 


rms 
Be ZNO ie 
Rh = 50.0 [%] 
P= 1021.5 + 0.2 [mb] 
R=5.30[m] (calculated from two way propagation path data for 
spheres two, three and four) 


|_ (referenced to 107 '2 [W/m?2]) 


94.7 + 0.1 [dB] 
94.5 
94.7 
25.0 
94.8 
95.1 
94.2 
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Table 12: One Way Propagation Path Data for 3.5 [V], 
5000 [Hz] Continuous Wave Supply from the HP3314A 
Function Generator; Collected 7 December 1987. 


(V__), = 0.523 + 0.001 [V] 


rms 
Ho Za 008 bata 
Rh = 50.0 [%] 
P = 1021.5 + 0.2 [mb] 
R=5.30[m] (calculated from two way propagation path data for 
spheres two, three and four) 


|_ (referenced to 10°! [W/m?}) 


97.4 + 0.1 [dB] 
96.4 
95. | 
97.0 
96.2 
97.1 
96.4 
95.4 
96.3 
97.4 


Ley 
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